Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,155,881 C→T 5.2% A294V (GCC→GTC)  argC → N‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,155,8810CT 5.2% 209.0 / 2.4 80A294V (GCC→GTC) argCN‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base T (2/2);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

GCCGCCTGAAATCGGGTGTGACCCAGGCGCAAGTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAATGATGAATCCATTAAT  >  NC_000913/4155734‑4156029
                                                                                                                                                   |                                                                                                                                                    
gccgccTCAAATCGGGTGTGACCCAGGCGCAACTCGCGCAAGTGTTACAGCAGGCGTATGCCCATACACCGATGGTGCGGCTGTATGACAACGGCGTTCCTGCGCAGAAAAATTTCGTTGGGCTGCCATTTTGCGATATCGGGATTGcc                                                                                                                                                     >  2:842666/1‑149 (MQ=255)
   gccTGAAATCGGCTGTGACCCAGGCGCAAGTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGtt                                                                                                                                                  <  1:383677/149‑1 (MQ=255)
         aaTCGGGTGTGACCCAGGCGCAAGTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGc                                                                                                                                            >  1:708819/1‑149 (MQ=255)
                ggggACCCAGGCGCAAGTCGCGCAAGTGATACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATc                                                                                                                                     <  2:480173/146‑1 (MQ=255)
                 gggACCCAGGCGCAAGTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATTACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCt                                                                                                                                    <  2:445868/147‑1 (MQ=255)
                    aCCCAGGCGCAAGTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGGCGTTCAGGGCGAGCATCTGat                                                                                                                                 >  1:265404/1‑149 (MQ=255)
                     cccAGGCGCAACTCGCGAAAGTGTAACAGCACGCGTAAGCCCATAAACCGCAGGGGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTGTGCCGTTCAGGGCGAGCATCTGatt                                                                                                                                <  1:702626/149‑1 (MQ=255)
                      ccAGGCGCAAGTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGTCGTTCAGGGCGAGCATCTGatta                                                                                                                               >  2:193099/1‑149 (MQ=255)
                      ccAGGCGCAAGTCGCGCAAGTGTTACAGCAGGAGTATGCCCATAAACCGCTGGTCCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGTCGTTCAGGGCGAGCATCTGatta                                                                                                                               >  2:646669/1‑149 (MQ=255)
                        tggCGGAAGTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACATAGGCGTTCCTGCGCTGAAAAATGTGGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGattatt                                                                                                                             <  1:99804/148‑1 (MQ=255)
                        cggCGCAAGGCGCGCCAGAGTACCACCAGTCGGATGCCCATACACCGCTGGTGCGGCGGTATGACACCGGCGTTCCGGCGCTGAATACTGGCGCTGGGCTGCCATTTTGCGAAATCGGGTTTGCCGTTCAGGGCGAGCATCTGattatt                                                                                                                             <  2:532867/148‑1 (MQ=255)
                         ggCGCAAGTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGATCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATtg                                                                                                                            >  1:712507/1‑149 (MQ=255)
                           cgcAAGTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATtgtg                                                                                                                          <  2:42273/149‑1 (MQ=255)
                            gcAAGTCGCGCAAGTGTTACACCCGGCGGATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCCTTCCGGCGCTGAAAAATGTCGTTGGTCTGCCATTTTGCTATATCGGGTTTGCCGTGCAGGGCGAGCATCTGATTATTGTgg                                                                                                                         >  2:788690/1‑149 (MQ=255)
                            gcAACTCTCTCAAGTGTTACATCAGGCGTATGCACATAAACCTCTGGTGCGTCTGTATGACAAAGGCGTTCAGTCGCTGAAAAATGTCGTTGGGCTTCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTgg                                                                                                                         <  2:859319/149‑1 (MQ=255)
                             cAAGTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCAGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGc                                                                                                                        >  1:372192/1‑149 (MQ=255)
                                gTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGAc                                                                                                                     >  1:309875/1‑149 (MQ=255)
                                   gcgcAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCATGGCGAGCATCTGATTATTGTGGCGACCGa                                                                                                                  >  2:160540/1‑149 (MQ=255)
                                   gcgcAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCGGAGTATTGTGGCGACCGa                                                                                                                  >  2:986857/1‑149 (MQ=255)
                                           gtTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACt                                                                                                          <  1:408000/149‑1 (MQ=255)
                                           gtTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACt                                                                                                          >  1:400325/1‑149 (MQ=255)
                                           gtTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACt                                                                                                          >  2:379331/1‑149 (MQ=255)
                                            ttACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCGGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACtt                                                                                                         >  1:949620/1‑149 (MQ=255)
                                               cagcagGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACt                                                                                                      <  1:190131/149‑1 (MQ=255)
                                                           gCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCgcggc                                                                                          >  1:696441/1‑149 (MQ=255)
                                                             ccATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTACCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGAGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGcggcgg                                                                                        <  2:375833/149‑1 (MQ=255)
                                                                  aaCCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACaa                                                                                   >  2:952783/1‑149 (MQ=255)
                                                                   aCCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAg                                                                                  >  1:644093/1‑149 (MQ=255)
                                                                   aCCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAg                                                                                  >  2:895067/1‑149 (MQ=255)
                                                                   aCCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCAAAGACAACTTACTGAAAGGCGCGGCGGCACAAg                                                                                  >  1:482392/1‑149 (MQ=255)
                                                                   aCCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGGCGTTGGGCTGCCATTTTGCGATATCGGGTTTGACGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAg                                                                                  >  2:248814/1‑149 (MQ=255)
                                                                   aCCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGAGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGGATCTGATTATTGTGGCGACCGAAGACAACTTAATGAAAGGCGCGGAGGCACAAg                                                                                  >  2:244879/1‑149 (MQ=255)
                                                                    ccGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACCAGc                                                                                 >  2:428438/1‑149 (MQ=255)
                                                                      ccTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTCGACGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCgg                                                                               <  1:353844/148‑1 (MQ=255)
                                                                                 tGTATGACAAAGGCGATCCGGCGCTGAAAAATGTCGTTGGGCTGCAATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATAGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAACCGGTACAGTGCGcc                                                                    >  2:715057/1‑149 (MQ=255)
                                                                                   tATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCaa                                                                  >  2:395643/1‑149 (MQ=255)
                                                                                    aTGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCTAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAat                                                                 >  1:163985/1‑149 (MQ=255)
                                                                                       aCAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGCGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATAtt                                                              >  1:762999/1‑149 (MQ=255)
                                                                                        cAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTc                                                             >  1:845386/1‑149 (MQ=255)
                                                                                        cAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAAAATTc                                                             >  1:595376/1‑149 (MQ=255)
                                                                                        cAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTAGCGATATCGGGTTTGCCGTTCAGCACGAGCATCTGATTAATGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTc                                                             >  1:246304/1‑149 (MQ=255)
                                                                                        cAAAGGCGTTCCGGCGCTGAAAAATCGCGTTGGGCTGCCATTTTGAGATCTCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTc                                                             <  2:106331/149‑1 (MQ=255)
                                                                                          aaGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGAATCTGATTATTGTGTCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCCCAAGCGGTACCGTGCTCCAATATTCGt                                                           >  2:993578/1‑149 (MQ=255)
                                                                                            ggCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGttt                                                         >  1:275982/1‑149 (MQ=255)
                                                                                             ctgTTCCTGCGATGAAAAATGTCGTTGGGCTGCCATTTTGCGATTTCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTc                                                        <  2:734153/147‑1 (MQ=255)
                                                                                               gTTCCGGCGCTGAAAAATGGCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCgg                                                      <  2:959764/149‑1 (MQ=255)
                                                                                                ttCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTc                                                        >  1:90065/1‑146 (MQ=255)
                                                                                                ttCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTc                                                        <  2:90065/146‑1 (MQ=255)
                                                                                                ttCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGACAATATTCGTTTCgag                                                     >  1:734153/1‑147 (MQ=255)
                                                                                                 accGGCGCTGAAAAATGGCGTTGGGCTGCCAGTTTGCGATATCGGGTGTGCCGTTCAGGGCGAGCAGCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCt                                                    <  2:868404/148‑1 (MQ=255)
                                                                                                    ggCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACTACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATg                                                 >  1:678320/1‑149 (MQ=255)
                                                                                                     gcgcTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGc                                                <  1:715864/149‑1 (MQ=255)
                                                                                                         tGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGaa                                            >  2:784472/1‑149 (MQ=255)
                                                                                                          gAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGaac                                           >  2:755227/1‑148 (MQ=255)
                                                                                                             aaaTGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGc                                        >  1:434278/1‑149 (MQ=255)
                                                                                                              aaTGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCa                                       <  1:825895/149‑1 (MQ=255)
                                                                                                                tGTCGTTGGGCTGCCATTTTGCGCTATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGt                                     <  1:132065/149‑1 (MQ=255)
                                                                                                                     ttGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCtt                                >  2:167141/1‑149 (MQ=255)
                                                                                                                     ttGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCtt                                <  2:696441/149‑1 (MQ=255)
                                                                                                                      tGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTa                               <  2:930865/149‑1 (MQ=255)
                                                                                                                        ggCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTAtt                             <  1:423630/149‑1 (MQ=255)
                                                                                                                           ggCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTaa                          <  2:220150/148‑1 (MQ=255)
                                                                                                                             ccATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAgg                        <  2:309875/149‑1 (MQ=255)
                                                                                                                              cATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAggg                       >  1:926724/1‑149 (MQ=255)
                                                                                                                              cATTTTGAGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCTGCACAAGCGGTAGAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAggg                       <  1:624616/149‑1 (MQ=255)
                                                                                                                                ttttGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTg                     >  2:291079/1‑149 (MQ=255)
                                                                                                                                 tttGCGAGTTAGGGTTTGCCGGTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGATGCGGCACAAGCGGTACAGTTCGCCAATATTCGTTTCGGCTATGCGGTAACGCAGTCTCTTATTTAAGGGTGc                    <  1:842666/149‑1 (MQ=255)
                                                                                                                                     cGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTAATTTAAGGGTGCAatg                >  1:167438/1‑149 (MQ=255)
                                                                                                                                      gATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAatga               <  2:400325/149‑1 (MQ=255)
                                                                                                                                      gATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTCCGCCAATATTCGTTTCGGCTATGCGGCAACGCAGTCGCTTATTTAAGGGTGCAatga               >  2:1017275/1‑149 (MQ=255)
                                                                                                                                       atatCGGGTTTGTCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAatgat              <  1:646669/149‑1 (MQ=255)
                                                                                                                                       atatCGGGTTTGTCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAatgat              <  1:193099/149‑1 (MQ=255)
                                                                                                                                       atatCGGGTTTGCCGTTCCGGGCGAGCATCTGATTATTGTGGCGACCAAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAatgat              <  2:482392/149‑1 (MQ=255)
                                                                                                                                       atatCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGACAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAatgat              <  2:15726/149‑1 (MQ=255)
                                                                                                                                       atatCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGTAACGCAGTCTCTTATTTAAGGGTGCTatgat              <  2:877211/149‑1 (MQ=255)
                                                                                                                                       atatCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAatgat              <  2:644093/149‑1 (MQ=255)
                                                                                                                                       atatCGGGTTTGCCGTTCAGGGCGAGAATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAatgat              <  2:753671/149‑1 (MQ=255)
                                                                                                                                             ggTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAGCGCAGTCTCTTATTTAAGGGTGCAATGATGAATcc        >  2:479014/1‑149 (MQ=255)
                                                                                                                                               tttGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAGGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAATGATGAATCCat      <  1:379331/149‑1 (MQ=255)
                                                                                                                                                 tGCCGTTCAGGGCGAGCATCTGATTATTGTTGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAATGATGAATCCatta    <  2:708819/149‑1 (MQ=255)
                                                                                                                                                  gCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAATGATGAATCCattaa   <  2:817544/149‑1 (MQ=255)
                                                                                                                                                   ccGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAATGATGAATCCattaat  <  1:60708/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCGCCTGAAATCGGGTGTGACCCAGGCGCAAGTCGCGCAAGTGTTACAGCAGGCGTATGCCCATAAACCGCTGGTGCGGCTGTATGACAAAGGCGTTCCGGCGCTGAAAAATGTCGTTGGGCTGCCATTTTGCGATATCGGGTTTGCCGTTCAGGGCGAGCATCTGATTATTGTGGCGACCGAAGACAACTTACTGAAAGGCGCGGCGGCACAAGCGGTACAGTGCGCCAATATTCGTTTCGGCTATGCGGAAACGCAGTCTCTTATTTAAGGGTGCAATGATGAATCCATTAAT  >  NC_000913/4155734‑4156029

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑