Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,370,206 C→T 5.8% R63R (CGC→CGT 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,2060CT 5.8% 227.0 / 2.7 69R63R (CGC→CGTarnDundecaprenyl phosphate‑alpha‑L‑ara4FN deformylase
Reads supporting (aligned to +/- strand):  ref base C (37/28);  new base T (2/2);  total (39/30)
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.72e-01

GTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCT  >  NC_000913/2370058‑2370354
                                                                                                                                                    |                                                                                                                                                    
gTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTgcgc                                                                                                                                                      <  3:637945/149‑1 (MQ=255)
  ggCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTc                                                                                                                                                    >  2:655658/1‑149 (MQ=255)
     aCCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCaaa                                                                                                                                                 >  4:618619/1‑149 (MQ=255)
                      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                                                                                 <  3:654624/149‑1 (MQ=255)
                                                                     ttttCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGAtt                                                                                 <  1:575656/149‑1 (MQ=255)
                                                                      tttCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTg                                                                                <  3:365597/149‑1 (MQ=255)
                                                                      tttCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTg                                                                                <  3:444449/149‑1 (MQ=255)
                                                                          aGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAaga                                                                                     >  2:619826/1‑140 (MQ=255)
                                                                          aGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAaga                                                                                     <  1:619826/140‑1 (MQ=255)
                                                                            cGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATg                                                                          <  2:416277/149‑1 (MQ=255)
                                                                              tCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCCGCGTTTGAGCGCAGCATct                                                                                                                                       >  2:621019/1‑85 (MQ=37)
                                                                              tCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCCGCGTTTGAGCGCAGCATct                                                                                                                                       <  1:621019/86‑2 (MQ=37)
                                                                                 gCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAg                                                                                                           <  1:121171/111‑1 (MQ=255)
                                                                                 gCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAg                                                                                                           >  2:121171/1‑111 (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              >  2:571862/1‑149 (MQ=255)
                                                                                                                                        gAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTggg              >  2:669974/1‑149 (MQ=255)
                                                                                                                                        gAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTggg              <  1:571862/149‑1 (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  <  3:106146/149‑1 (MQ=255)
                                                                                                                                                    cTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCt  <  4:50128/149‑1 (MQ=255)
                                                                                                                                                    cTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCt  >  3:50128/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTACGGCCTGGCCAGGTAAAGAGATTGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCT  >  NC_000913/2370058‑2370354

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: