Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,677,835 C→T 6.3% R172R (CGG→CGA)  yhjG ← AsmA family protein YhjG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,677,8350CT 6.3% 293.1 / 6.3 95R172R (CGG→CGA) yhjGAsmA family protein YhjG
Reads supporting (aligned to +/- strand):  ref base C (48/40);  new base T (4/2);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

TGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCG  >  NC_000913/3677686‑3677982
                                                                                                                                                     |                                                                                                                                                    
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:1473641/150‑1 (MQ=255)
           gCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATc                                                                                                                                           <  1:676262/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                                                                                                                                 >  1:635101/1‑150 (MQ=255)
                     tAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAata                                                                                                                                 >  2:124145/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:1444280/1‑150 (MQ=255)
                          gCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATc                                                                                                                            <  1:67922/150‑1 (MQ=255)
                            ccACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCa                                                                                                                          <  1:1164988/150‑1 (MQ=255)
                             cACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAATGGATCAACAAAAATCTCCAGATCCGCTTTTCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAATAATATTATCCAg                                                                                                                         <  2:689992/150‑1 (MQ=255)
                               ccTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGcc                                                                                                                       <  2:1481415/150‑1 (MQ=255)
                               ccTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGcc                                                                                                                       <  2:195447/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                                                                                                     >  2:62648/1‑150 (MQ=255)
                                                 tACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACg                                                                                                     <  1:827967/150‑1 (MQ=255)
                                                       tCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑g                                                                                              <  1:324611/150‑1 (MQ=255)
                                                        cGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑Gc                                                                                             >  1:1229775/1‑150 (MQ=255)
                                                          aTCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAACCGACCATGCCGACGGCTTT‑GCAt                                                                                           >  2:1171015/1‑150 (MQ=255)
                                                           tCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCAtt                                                                                          >  1:266891/1‑150 (MQ=255)
                                                             cAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTcg                                                                                        <  1:808476/150‑1 (MQ=255)
                                                              aGTAACTTCGCTGAACGGCAGCGGCTTGCTTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTcgc                                                                                       >  1:1613453/1‑150 (MQ=255)
                                                              aGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTcgc                                                                                       >  1:1495235/1‑150 (MQ=255)
                                                               gTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCCCTTTGCTTACTTTGTCATCAATGGCGATACGCCGTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTcgcg                                                                                      >  2:193225/1‑150 (MQ=255)
                                                                    ttCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCttt                                                                                 <  2:126955/150‑1 (MQ=255)
                                                                     tCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTg                                                                                <  2:1177721/150‑1 (MQ=255)
                                                                      cGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGt                                                                               <  2:811328/150‑1 (MQ=255)
                                                                      cGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGt                                                                               <  1:150652/150‑1 (MQ=255)
                                                                       gCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGtt                                                                              <  1:124145/150‑1 (MQ=255)
                                                                                cAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGtt                                                                     >  1:1097685/1‑150 (MQ=255)
                                                                                cAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGtt                                                                     <  1:1501515/150‑1 (MQ=255)
                                                                                 aGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTg                                                                    <  2:1099886/150‑1 (MQ=255)
                                                                                 aGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTg                                                                    <  2:266891/150‑1 (MQ=255)
                                                                                    ggCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCa                                                                 >  2:1032095/1‑150 (MQ=255)
                                                                                         gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                            >  2:871159/1‑150 (MQ=255)
                                                                                         gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                            >  1:146439/1‑150 (MQ=255)
                                                                                         gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                            <  1:417537/150‑1 (MQ=255)
                                                                                         gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                            >  1:1059484/1‑150 (MQ=255)
                                                                                            tAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTaaa                                                         >  1:1246928/1‑150 (MQ=255)
                                                                                            tAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTaaa                                                         >  2:1102219/1‑150 (MQ=255)
                                                                                              aGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACg                                                       <  1:1123575/150‑1 (MQ=255)
                                                                                               tggATCAACAAAAATCTCCAGATACGCTTTGATTACTTTCTCATCAATGGCTATCCTCCCTTTATAGAAAAGAATATTATCCAGAAGAAACCACCATTACGACCGCTAT‑GCATTCGACTAATTCTTATCAACGTTTCCAACATTAAACGt                                                      >  2:966761/2‑150 (MQ=255)
                                                                                               gggATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGt                                                      >  2:917422/1‑150 (MQ=255)
                                                                                               gggATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGt                                                      >  1:1377324/1‑150 (MQ=255)
                                                                                                ggATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGCCCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTc                                                     >  2:1105184/1‑150 (MQ=255)
                                                                                                ggATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTc                                                     >  2:1006753/1‑150 (MQ=255)
                                                                                                 gATCAACAAAAATCTCCAGATACGCTTTGCTTACTTTGTCATCAATGGCGATCGGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGGCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTcc                                                    <  2:1077898/150‑1 (MQ=255)
                                                                                                  aTCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCa                                                   <  1:966761/150‑1 (MQ=255)
                                                                                                  aTCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCa                                                   >  2:138221/1‑150 (MQ=255)
                                                                                                   tCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAg                                                  <  1:138221/150‑1 (MQ=255)
                                                                                                    caacaaAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGt                                                 <  2:77694/150‑1 (MQ=255)
                                                                                                     aacaaAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGtt                                                >  2:175165/1‑150 (MQ=255)
                                                                                                       caaAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGatat                                              >  1:1399706/1‑150 (MQ=255)
                                                                                                       caaAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTAt                                              >  1:1419725/1‑150 (MQ=255)
                                                                                                             tgtcCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCt                                        >  1:525286/3‑150 (MQ=255)
                                                                                                              ctcCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCtt                                       <  1:1419616/150‑1 (MQ=255)
                                                                                                              ctcCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCtt                                       >  2:687094/1‑150 (MQ=255)
                                                                                                                ccAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCttt                                      <  1:62648/150‑1 (MQ=255)
                                                                                                                     tCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCagag                                <  1:1572653/150‑1 (MQ=255)
                                                                                                                      ccGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCagaga                               >  1:1281784/1‑150 (MQ=255)
                                                                                                                      ccGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCagaga                               <  2:516793/150‑1 (MQ=255)
                                                                                                                        gCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAgg                             >  1:491951/1‑150 (MQ=255)
                                                                                                                          tttGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCg                           <  2:1453262/150‑1 (MQ=255)
                                                                                                                               ttACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTGAGAGAGGCGAATCa                      <  2:1399706/150‑1 (MQ=255)
                                                                                                                                tACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAg                     >  1:438946/1‑150 (MQ=255)
                                                                                                                                 aCTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAgg                    >  2:39362/1‑150 (MQ=255)
                                                                                                                                   tttGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGgcg                  <  2:1097685/150‑1 (MQ=255)
                                                                                                                                       tcatcaATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGcgcgcg              <  2:1229775/150‑1 (MQ=255)
                                                                                                                                       cgatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGcgcgcg              >  2:590916/3‑150 (MQ=255)
                                                                                                                                        gatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt             >  1:34338/2‑150 (MQ=255)
                                                                                                                                        catcaATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt             >  2:1372400/1‑150 (MQ=255)
                                                                                                                                         atcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTc            >  1:1123521/1‑150 (MQ=255)
                                                                                                                                         atcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTc            >  1:2813/1‑150 (MQ=255)
                                                                                                                                               ggCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCtt      >  2:286004/1‑150 (MQ=255)
                                                                                                                                               ggCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCtt      >  1:346546/1‑150 (MQ=255)
                                                                                                                                               ggCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTCGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCtt      >  1:878562/1‑150 (MQ=255)
                                                                                                                                                  gATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGAATAAACGACCAGTTATTGTTCTTTTAAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTc   >  1:825581/1‑150 (MQ=255)
                                                                                                                                                   aTCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCg  <  1:935726/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCG  >  NC_000913/3677686‑3677982

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: ‑