Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,370,215 G→A 5.5% A66A (GCG→GCA)  arnD → undecaprenyl phosphate‑alpha‑L‑ara4FN deformylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,370,2150GA 5.5% 242.0 / 3.5 73A66A (GCG→GCA) arnDundecaprenyl phosphate‑alpha‑L‑ara4FN deformylase
Reads supporting (aligned to +/- strand):  ref base G (38/31);  new base A (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

GCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCC  >  NC_000913/2370080‑2370363
                                                                                                                                       |                                                                                                                                                    
gcgTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATg                                                                                                                                         <  4:224865/149‑1 (MQ=255)
   tCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCt                                                                                                                                      >  3:436276/1‑149 (MQ=255)
          gAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGGCAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGatat                                                                                                                               >  2:504507/1‑149 (MQ=255)
          gAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGatat                                                                                                                               >  1:416277/1‑149 (MQ=255)
              tCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTa                                                                                                                           >  2:123133/1‑149 (MQ=255)
                ttGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACt                                                                                                                         >  4:365597/1‑149 (MQ=255)
                ttGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACCCGCCATCGCTTTATGGCTGGGATATTTTACt                                                                                                                         <  4:424356/149‑1 (MQ=255)
                 tGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTg                                                                                                                        <  1:655658/149‑1 (MQ=255)
                    gTAAGCATAATATTCAGGCCAGCATTTTTTTCAACGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCa                                                                                                                     >  1:196270/1‑149 (MQ=255)
                     tAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAg                                                                                                                    >  4:156760/1‑149 (MQ=255)
                      aaGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAgg                                                                                                                   <  4:293649/149‑1 (MQ=255)
                         cATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAc                                                                                                                <  3:156760/149‑1 (MQ=255)
                           taataTTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACgg                                                                                                              >  4:444449/1‑149 (MQ=255)
                           taataTTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACgg                                                                                                              >  1:344970/1‑149 (MQ=255)
                           taataTTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCCGCGTTTGAGCGCAGCATcttc                                                                                                                                              >  4:266410/1‑114 (MQ=255)
                           taataTTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCCGCGTTTGAGCGCAGCATcttc                                                                                                                                              <  3:266410/117‑4 (MQ=255)
                              tatTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCt                                                                                                           <  3:618619/149‑1 (MQ=255)
                                ttCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGGCAAGATGAGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTgg                                                                                                         >  1:402975/1‑149 (MQ=255)
                                ttCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTgg                                                                                                         >  2:443418/1‑149 (MQ=255)
                                       cAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTa                                                                                                  >  4:61197/1‑149 (MQ=255)
                                        aGCATTTTTTTCAGCGGCGGGCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTaa                                                                                                 >  4:628120/1‑149 (MQ=255)
                                         gCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTaaa                                                                                                <  4:356591/149‑1 (MQ=255)
                                         gCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTaaa                                                                                                <  1:555971/149‑1 (MQ=255)
                                          cATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGAACGGCCTGGCCAGGTAAag                                                                                               <  4:436276/149‑1 (MQ=255)
                                            tttttttCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAagag                                                                                             >  1:190272/1‑149 (MQ=255)
                                              tttttCAGCGTCGGCCCGGACAATATGGTCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGAt                                                                                           <  4:339997/149‑1 (MQ=255)
                                               ttttCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGAtt                                                                                          <  1:575656/149‑1 (MQ=255)
                                               ttttCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGAtt                                                                                          <  3:654624/149‑1 (MQ=255)
                                                tttCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTg                                                                                         <  3:365597/149‑1 (MQ=255)
                                                tttCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTg                                                                                         <  3:444449/149‑1 (MQ=255)
                                                    aGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAaga                                                                                              <  1:619826/140‑1 (MQ=255)
                                                    aGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAaga                                                                                              >  2:619826/1‑140 (MQ=255)
                                                      cGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATg                                                                                   <  2:416277/149‑1 (MQ=255)
                                                        tCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCCGCGTTTGAGCGCAGCATct                                                                                                                                                <  1:621019/86‑2 (MQ=37)
                                                        tCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCCGCGTTTGAGCGCAGCATct                                                                                                                                                >  2:621019/1‑85 (MQ=37)
                                                           gCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAg                                                                                                                    >  2:121171/1‑111 (MQ=255)
                                                           gCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAg                                                                                                                    <  1:121171/111‑1 (MQ=255)
                                                                       tGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCAtt                                                                  <  3:628120/149‑1 (MQ=255)
                                                                       tGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCAtt                                                                  <  2:196270/149‑1 (MQ=255)
                                                                        gggCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTc                                                                 >  2:554732/1‑149 (MQ=255)
                                                                        gggCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTc                                                                 <  3:508408/149‑1 (MQ=255)
                                                                        gggCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGGACGGCCTGGCCAGGTTAAGAGATTGGTCATGCAAATGCCGATATCATTc                                                                 >  2:178467/1‑149 (MQ=255)
                                                                           ccgccATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTg                                                              >  2:516194/1‑149 (MQ=255)
                                                                            cgccATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGa                                                             >  4:91137/1‑149 (MQ=255)
                                                                              ccATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAg                                                           >  4:213695/1‑149 (MQ=255)
                                                                                  ctctGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGc                                                       <  1:554732/149‑1 (MQ=255)
                                                                                        cGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACa                                                 <  1:123133/149‑1 (MQ=255)
                                                                                         gACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAg                                                                                                     <  3:172221/96‑1 (MQ=255)
                                                                                         gACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAg                                                                                                     >  4:172221/1‑96 (MQ=255)
                                                                                            tGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCAc                                             >  4:203274/1‑149 (MQ=255)
                                                                                             ggTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACg                                            <  1:443418/149‑1 (MQ=255)
                                                                                                 aaGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGt                                        >  1:355016/1‑149 (MQ=255)
                                                                                                      acaGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCgg                                                        >  1:617387/1‑128 (MQ=255)
                                                                                                      acaGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCgg                                                        <  2:617387/128‑1 (MQ=255)
                                                                                                      acaGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGcc                                   >  2:82769/1‑149 (MQ=255)
                                                                                                      acaGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGcc                                   >  1:626097/1‑149 (MQ=255)
                                                                                                      acaGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGcc                                   >  1:654487/1‑149 (MQ=255)
                                                                                                      acaGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGcc                                   >  1:7735/1‑149 (MQ=255)
                                                                                                       caGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCt                                  >  3:297716/1‑149 (MQ=255)
                                                                                                          tttttGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCa                               >  3:280166/1‑149 (MQ=255)
                                                                                                                  gAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTggg                       <  1:571862/149‑1 (MQ=255)
                                                                                                                  gAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTggg                       >  2:571862/1‑149 (MQ=255)
                                                                                                                  gAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTggg                       >  2:669974/1‑149 (MQ=255)
                                                                                                                   aaGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGa                      <  2:654487/149‑1 (MQ=255)
                                                                                                                     gATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATc                    >  4:415170/1‑149 (MQ=255)
                                                                                                                       tgctgcGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCAc                  <  3:415170/149‑1 (MQ=255)
                                                                                                                          tgcGCTCAAACGCGGCATCGCTTGATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCAt               >  2:531753/1‑149 (MQ=255)
                                                                                                                             gcTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTgg                        <  2:520931/137‑1 (MQ=255)
                                                                                                                             gcTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTgg                        >  1:520931/1‑137 (MQ=255)
                                                                                                                              cTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCt           <  4:50128/149‑1 (MQ=255)
                                                                                                                              cTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCt           >  3:50128/1‑149 (MQ=255)
                                                                                                                              cTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCt           <  3:106146/149‑1 (MQ=255)
                                                                                                                                 aaaCGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGc        >  1:57491/1‑149 (MQ=255)
                                                                                                                                       ggCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGccc  <  2:344970/149‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                    
GCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCC  >  NC_000913/2370080‑2370363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑