Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,140,862 A→T 6.7% L142H (CTT→CAT) ‡ hybD ← putative hydrogenase 2 maturation protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,140,8620AT 6.7% 229.4 / 5.8 75L142H (CTT→CAT) ‡hybDputative hydrogenase 2 maturation protease
Reads supporting (aligned to +/- strand):  ref base A (38/32);  new base T (2/3);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

TGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGATGC  >  NC_000913/3140714‑3141011
                                                                                                                                                    |                                                                                                                                                     
tGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCaac                                                                                                                                                      >  1:118182/1‑149 (MQ=255)
tGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAg                                                                                                                                                      <  2:77048/150‑1 (MQ=255)
 gCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCCGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgc                                                                                                                                                     >  2:810423/1‑150 (MQ=255)
  cATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgcg                                                                                                                                                    >  1:279975/1‑150 (MQ=255)
         cGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTc                                                                                                                                             <  2:788478/150‑1 (MQ=255)
         cGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTc                                                                                                                                             >  2:752542/1‑150 (MQ=255)
         cGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTc                                                                                                                                             >  2:757554/1‑150 (MQ=255)
           ccGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCGCTTGAGCAGGTTctg                                                                                                                                           <  1:1261166/150‑3 (MQ=255)
             gggCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCGCTTGAGCAGGTTctggc                                                                                                                                         <  1:401698/150‑5 (MQ=255)
             gggCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCGCTTGAGCAGGTTctggc                                                                                                                                         <  2:1003414/150‑5 (MQ=255)
              ggCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGCATCGCCTCCCGTGGGATGGCTTCCCCGCCAGATTCCCGCCGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCa                                                                                                                                        >  1:238660/1‑150 (MQ=255)
                cAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCAtt                                                                                                                                      >  1:1051792/1‑150 (MQ=255)
                cAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATAGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCAtt                                                                                                                                      >  2:359133/1‑150 (MQ=255)
                 aaTTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTg                                                                                                                                     <  1:118226/150‑1 (MQ=255)
                 aaTTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTg                                                                                                                                     <  2:720797/150‑1 (MQ=255)
                      cttcAAACGCTGTCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCa                                                                                                                                >  2:1177831/1‑150 (MQ=255)
                      cttcAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCa                                                                                                                                >  1:307585/1‑150 (MQ=255)
                      cttcAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCa                                                                                                                                >  1:1185382/1‑150 (MQ=255)
                       ttcAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCaa                                                                                                                               >  2:1054666/1‑150 (MQ=255)
                         cAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAgg                                                                                                                                                <  2:1181492/131‑1 (MQ=255)
                         cAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAgg                                                                                                                                                >  1:1181492/1‑131 (MQ=255)
                         cAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGGTCAATCATTGCTTCAAcc                                                                                                                             >  2:248981/1‑150 (MQ=255)
                          aaaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCg                                                                                                                            >  1:1351637/1‑150 (MQ=255)
                                 gCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGt                                                                                                                     >  1:516508/1‑150 (MQ=255)
                                     gTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAg                                                                                                                 <  1:1177831/150‑1 (MQ=255)
                                      tACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGc                                                                                                                <  2:1176264/150‑1 (MQ=255)
                                          tGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCg                                                                                                                      >  1:235234/1‑140 (MQ=255)
                                          tGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCg                                                                                                                      <  2:235234/140‑1 (MQ=255)
                                              ccTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGAtgtgt                                                                                                        >  1:111971/1‑150 (MQ=255)
                                                ttCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTgg                                                                                                      >  1:1224068/1‑150 (MQ=255)
                                                 tCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGc                                                                                                     <  2:513613/150‑1 (MQ=255)
                                                      gAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAg                                                                                                <  1:794756/150‑1 (MQ=255)
                                                           cTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGAtt                                                                                           >  1:555216/1‑150 (MQ=255)
                                                           cTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGAtt                                                                                           >  1:1360425/1‑150 (MQ=255)
                                                           cTGGAAACCTGATATATCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACCCAAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGAtt                                                                                           >  2:1403548/1‑150 (MQ=255)
                                                             ggAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTcc                                                                                         >  1:715638/1‑150 (MQ=255)
                                                             ggAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTcc                                                                                         >  2:948958/1‑150 (MQ=255)
                                                                aaCCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGTTTCAATCTTTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCtgg                                                                                      <  2:434191/150‑1 (MQ=255)
                                                                aaCCTGCTATCTCTCCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCggg                                                                                      <  2:961510/150‑1 (MQ=255)
                                                                   cTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGAta                                                                                   <  1:757554/150‑2 (MQ=255)
                                                                   cTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATc                                                                                   <  1:752542/150‑1 (MQ=255)
                                                                        aTCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGcc                                                                              >  2:997105/1‑150 (MQ=255)
                                                                         tctcTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCg                                                                             <  2:1343386/150‑1 (MQ=255)
                                                                         tctcTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCg                                                                             <  2:1379746/150‑1 (MQ=255)
                                                                                aGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAggg                                                                      <  2:97100/150‑1 (MQ=255)
                                                                                 gTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGt                                                                     >  1:1411765/1‑150 (MQ=255)
                                                                                  tCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTc                                                                    <  2:419895/150‑1 (MQ=255)
                                                                                   cATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCa                                                                   >  2:535206/1‑150 (MQ=255)
                                                                                   aaTGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCa                                                                   <  2:490774/149‑1 (MQ=255)
                                                                                     tGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGc                                                                 <  1:948958/150‑1 (MQ=255)
                                                                                       aaTCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCtt                                                               >  1:910795/1‑150 (MQ=255)
                                                                                            ccTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCg                                                          <  2:1185382/150‑1 (MQ=255)
                                                                                            ccTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCg                                                          <  1:810423/150‑1 (MQ=255)
                                                                                               cccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCgg                                                         >  1:872205/1‑148 (MQ=255)
                                                                                               cccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCgg                                                         <  2:872205/148‑1 (MQ=255)
                                                                                                  gTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACt                                                    <  1:639778/150‑1 (MQ=255)
                                                                                                     ggATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGc                                                 >  2:532109/1‑150 (MQ=255)
                                                                                                          gCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTg                                            >  2:122626/1‑150 (MQ=255)
                                                                                                                   ccAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAgg                                   >  1:554158/1‑150 (MQ=255)
                                                                                                                            cgcagcgcagCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAg                          <  2:1360425/150‑1 (MQ=255)
                                                                                                                              cagcgcagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGAc                        <  2:516508/150‑1 (MQ=255)
                                                                                                                                gcgcagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGt                      <  2:910795/150‑1 (MQ=255)
                                                                                                                                 cgcagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTc                     >  1:1176802/1‑150 (MQ=255)
                                                                                                                                 cgcagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTc                     >  2:706672/1‑150 (MQ=255)
                                                                                                                                   cagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCgg                   <  1:1403548/150‑1 (MQ=255)
                                                                                                                                   cagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCgg                   <  2:279975/150‑1 (MQ=255)
                                                                                                                                     gcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcc                 >  2:961415/1‑150 (MQ=255)
                                                                                                                                     gcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcc                 >  1:1205899/1‑150 (MQ=255)
                                                                                                                                     gattgAACCTGCGCTTGAGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcc                 >  1:1289635/5‑150 (MQ=255)
                                                                                                                                          aaCCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccaggcc            <  2:1051792/150‑1 (MQ=255)
                                                                                                                                          aaCCTGCGCTTGAGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccaggcc            >  1:251233/1‑150 (MQ=255)
                                                                                                                                           aCCTGCTCAAGCGCAGGGTGAATCATTGCTTCAAACCTCTTCGTTAAGCCCATGTGTGGCTCCAGCGATTCCGGGATAACGCCGCCCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAATCGCAGGGCCGACAACACGTCggccaggcca           >  2:773224/1‑150 (MQ=255)
                                                                                                                                            ccTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCaa          <  2:156293/150‑1 (MQ=255)
                                                                                                                                              tGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGc        <  1:473544/150‑1 (MQ=255)
                                                                                                                                                  cAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGAt    >  1:603270/1‑150 (MQ=255)
                                                                                                                                                  cAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGAt    >  2:971374/1‑150 (MQ=255)
                                                                                                                                                   aaGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGATg   <  1:997105/150‑1 (MQ=255)
                                                                                                                                                    aGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGATGc  >  1:225746/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGATGC  >  NC_000913/3140714‑3141011

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