Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,677,828 T→G 6.4% I175L (ATT→CTT)  yhjG ← AsmA family protein YhjG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,677,8280TG 6.4% 292.0 / 10.3 94I175L (ATT→CTT) yhjGAsmA family protein YhjG
Reads supporting (aligned to +/- strand):  ref base T (47/41);  new base G (4/2);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CGTCCCTGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTC  >  NC_000913/3677680‑3677972
                                                                                                                                                    |                                                                                                                                                 
cGTCCCTGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATACGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATg                                                                                                                                                  >  2:808476/1‑150 (MQ=255)
 gTCCCTGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATgg                                                                                                                                                 <  2:1546692/150‑1 (MQ=255)
 gTCCCTGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATgg                                                                                                                                                 <  1:970099/150‑1 (MQ=255)
      tGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATc                                                                                                                                            >  2:321493/1‑150 (MQ=255)
           cTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGccc                                                                                                                                       >  1:1089922/1‑150 (MQ=255)
           cTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGccc                                                                                                                                       >  1:516793/1‑150 (MQ=255)
              cAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTg                                                                                                                                    <  2:920719/150‑1 (MQ=255)
                 gCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATc                                                                                                                                 <  1:676262/150‑1 (MQ=255)
                 gCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATc                                                                                                                                 <  1:1473641/150‑1 (MQ=255)
                 gCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAAGGGCGGATCGCCATTGAtg                                                                                                                                 <  1:489606/150‑2 (MQ=255)
                  ccAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCg                                                                                                                                >  2:1419616/1‑150 (MQ=255)
                    aaaaaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGaa                                                                                                                              >  2:172458/1‑150 (MQ=255)
                      aaaCGTAATCGCCCACCTTTTATTTATACCGTTTACCTTTCGAGCCAGTAACTTCGCTGATCGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTGTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGaaaa                                                                                                                            >  2:1479185/1‑150 (MQ=255)
                       aaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAAGGGCGGATCGCCATTGAtgacaaag                                                                                                                           <  1:655946/150‑8 (MQ=255)
                           tAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAata                                                                                                                       >  2:124145/1‑150 (MQ=255)
                           tAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAata                                                                                                                       >  1:635101/1‑150 (MQ=255)
                           tAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAata                                                                                                                       >  2:1123575/1‑150 (MQ=255)
                            aaTCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAATGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAatat                                                                                                                      >  1:49453/1‑150 (MQ=255)
                                gCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGTTCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATc                                                                                                                  <  2:1610102/150‑1 (MQ=255)
                                gCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATc                                                                                                                  <  1:67922/150‑1 (MQ=255)
                                gCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATc                                                                                                                  >  1:1444280/1‑150 (MQ=255)
                                  ccACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCa                                                                                                                <  1:1164988/150‑1 (MQ=255)
                                   cACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAATGGATCAACAAAAATCTCCAGATCCGCTTTTCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAATAATATTATCCAg                                                                                                               <  2:689992/150‑1 (MQ=255)
                                     ccTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGcc                                                                                                             <  2:195447/150‑1 (MQ=255)
                                     ccTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGcc                                                                                                             <  2:1481415/150‑1 (MQ=255)
                                       ttttCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGa                                                                                                           <  1:1619439/150‑1 (MQ=255)
                                        tttCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGaa                                                                                                          >  2:447461/1‑150 (MQ=255)
                                        tttCTTTATCCGCTTTACATTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGaa                                                                                                          >  2:935726/1‑150 (MQ=255)
                                              tATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGAcc                                                                                                    <  2:1144065/150‑1 (MQ=255)
                                               aTCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCa                                                                                                   >  1:254851/1‑150 (MQ=255)
                                               aTCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCa                                                                                                   >  1:1099886/1‑150 (MQ=255)
                                               aTCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCa                                                                                                   <  2:839776/150‑1 (MQ=255)
                                                tCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCAt                                                                                                  >  1:1453262/1‑150 (MQ=255)
                                                  cGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGc                                                                                                <  2:1144450/150‑1 (MQ=255)
                                                   gCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGcc                                                                                               <  1:1541360/150‑1 (MQ=255)
                                                       tACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACg                                                                                           <  1:827967/150‑1 (MQ=255)
                                                       tACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACg                                                                                           >  2:62648/1‑150 (MQ=255)
                                                             tCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTg                                                                                    <  1:324611/150‑1 (MQ=255)
                                                              cGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGc                                                                                   >  1:1229775/1‑150 (MQ=255)
                                                                aTCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAACCGACCATGCCGACGGCTTTGCAt                                                                                 >  2:1171015/1‑150 (MQ=255)
                                                                 tCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCAtt                                                                                >  1:266891/1‑150 (MQ=255)
                                                                   cAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTcg                                                                              <  1:808476/150‑1 (MQ=255)
                                                                    aGTAACTTCGCTGAACGGCAGCGGCTTGCTTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTcgc                                                                             >  1:1613453/1‑150 (MQ=255)
                                                                    aGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTcgc                                                                             >  1:1495235/1‑150 (MQ=255)
                                                                     gTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCCCTTTGCTTACTTTGTCATCAATGGCGATACGCCGTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTcgcg                                                                            >  2:193225/1‑150 (MQ=255)
                                                                          ttCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCttt                                                                       <  2:126955/150‑1 (MQ=255)
                                                                           tCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTg                                                                      <  2:1177721/150‑1 (MQ=255)
                                                                            cGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGt                                                                     <  1:150652/150‑1 (MQ=255)
                                                                            cGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGt                                                                     <  2:811328/150‑1 (MQ=255)
                                                                             gCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGtt                                                                    <  1:124145/150‑1 (MQ=255)
                                                                                      cAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGtt                                                           <  1:1501515/150‑1 (MQ=255)
                                                                                      cAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGtt                                                           >  1:1097685/1‑150 (MQ=255)
                                                                                       aGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTg                                                          <  2:266891/150‑1 (MQ=255)
                                                                                       aGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTg                                                          <  2:1099886/150‑1 (MQ=255)
                                                                                          ggCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCa                                                       >  2:1032095/1‑150 (MQ=255)
                                                                                               gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                  >  1:146439/1‑150 (MQ=255)
                                                                                               gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                  >  2:871159/1‑150 (MQ=255)
                                                                                               gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                  <  1:417537/150‑1 (MQ=255)
                                                                                               gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                  >  1:1059484/1‑150 (MQ=255)
                                                                                                  tAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTaaa                                               >  1:1246928/1‑150 (MQ=255)
                                                                                                  tAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTaaa                                               >  2:1102219/1‑150 (MQ=255)
                                                                                                    aGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACg                                             <  1:1123575/150‑1 (MQ=255)
                                                                                                     tggATCAACAAAAATCTCCAGATACGCTTTGATTACTTTCTCATCAATGGCTATCCTCCCTTTATAGAAAAGAATATTATCCAGAAGAAACCACCATTACGACCGCTATGCATTCGACTAATTCTTATCAACGTTTCCAACATTAAACGt                                            >  2:966761/2‑150 (MQ=255)
                                                                                                     gggATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGt                                            >  1:1377324/1‑150 (MQ=255)
                                                                                                     gggATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGt                                            >  2:917422/1‑150 (MQ=255)
                                                                                                      ggATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGCCCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTc                                           >  2:1105184/1‑150 (MQ=255)
                                                                                                      ggATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTc                                           >  2:1006753/1‑150 (MQ=255)
                                                                                                       gATCAACAAAAATCTCCAGATACGCTTTGCTTACTTTGTCATCAATGGCGATCGGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGGCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTcc                                          <  2:1077898/150‑1 (MQ=255)
                                                                                                        aTCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCa                                         >  2:138221/1‑150 (MQ=255)
                                                                                                        aTCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCa                                         <  1:966761/150‑1 (MQ=255)
                                                                                                         tCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAg                                        <  1:138221/150‑1 (MQ=255)
                                                                                                          caacaaAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGt                                       <  2:77694/150‑1 (MQ=255)
                                                                                                           aacaaAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGtt                                      >  2:175165/1‑150 (MQ=255)
                                                                                                             caaAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGatat                                    >  1:1399706/1‑150 (MQ=255)
                                                                                                             caaAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTAt                                    >  1:1419725/1‑150 (MQ=255)
                                                                                                                   tgtcCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCt                              >  1:525286/3‑150 (MQ=255)
                                                                                                                    ctcCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCtt                             >  2:687094/1‑150 (MQ=255)
                                                                                                                    ctcCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCtt                             <  1:1419616/150‑1 (MQ=255)
                                                                                                                      ccAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCttt                            <  1:62648/150‑1 (MQ=255)
                                                                                                                           tCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCagag                      <  1:1572653/150‑1 (MQ=255)
                                                                                                                            ccGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCagaga                     >  1:1281784/1‑150 (MQ=255)
                                                                                                                            ccGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCagaga                     <  2:516793/150‑1 (MQ=255)
                                                                                                                              gCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAgg                   >  1:491951/1‑150 (MQ=255)
                                                                                                                                tttGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCg                 <  2:1453262/150‑1 (MQ=255)
                                                                                                                                     ttACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTGAGAGAGGCGAATCa            <  2:1399706/150‑1 (MQ=255)
                                                                                                                                      tACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAg           >  1:438946/1‑150 (MQ=255)
                                                                                                                                       aCTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAgg          >  2:39362/1‑150 (MQ=255)
                                                                                                                                         tttGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGgcg        <  2:1097685/150‑1 (MQ=255)
                                                                                                                                             tcatcaATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGcgcgcg    <  2:1229775/150‑1 (MQ=255)
                                                                                                                                             cgatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGcgcgcg    >  2:590916/3‑150 (MQ=255)
                                                                                                                                              gatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt   >  1:34338/2‑150 (MQ=255)
                                                                                                                                              catcaATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt   >  2:1372400/1‑150 (MQ=255)
                                                                                                                                               atcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTc  >  1:2813/1‑150 (MQ=255)
                                                                                                                                               atcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTc  >  1:1123521/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGTCCCTGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTC  >  NC_000913/3677680‑3677972

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: