Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,677,840 C→A 6.3% G171W (GGG→TGG)  yhjG ← AsmA family protein YhjG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,677,8400CA 6.3% 313.3 / 7.8 96G171W (GGG→TGG) yhjGAsmA family protein YhjG
Reads supporting (aligned to +/- strand):  ref base C (46/43);  new base A (4/2);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

CTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGC  >  NC_000913/3677691‑3677985
                                                                                                                                                     |                                                                                                                                                  
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:1123575/1‑150 (MQ=255)
                tAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAata                                                                                                                                    >  1:635101/1‑150 (MQ=255)
                tAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAata                                                                                                                                    >  2:124145/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: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:1099886/1‑150 (MQ=255)
                                    aTCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCa                                                                                                                >  1:254851/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                                                                                  <  1:150652/150‑1 (MQ=255)
                                                                 cGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGt                                                                                  <  2:811328/150‑1 (MQ=255)
                                                                  gCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGtt                                                                                 <  1:124145/150‑1 (MQ=255)
                                                                           cAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGtt                                                                        <  1:1501515/150‑1 (MQ=255)
                                                                           cAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGtt                                                                        >  1:1097685/1‑150 (MQ=255)
                                                                            aGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTg                                                                       <  2:266891/150‑1 (MQ=255)
                                                                            aGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTg                                                                       <  2:1099886/150‑1 (MQ=255)
                                                                               ggCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCa                                                                    >  2:1032095/1‑150 (MQ=255)
                                                                                    gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                               >  2:871159/1‑150 (MQ=255)
                                                                                    gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                               >  1:1059484/1‑150 (MQ=255)
                                                                                    gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                               <  1:417537/150‑1 (MQ=255)
                                                                                    gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                               >  1:146439/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                                                         >  1:1377324/1‑150 (MQ=255)
                                                                                          gggATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGt                                                         >  2:917422/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                                                      >  2:138221/1‑150 (MQ=255)
                                                                                             aTCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCa                                                      <  1:966761/150‑1 (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:2813/1‑150 (MQ=255)
                                                                                                                                    atcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTc               >  1:1123521/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)
                                                                                                                                                 ggCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCTTTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTGCAGTTATTGTTCTTTTCAGAGAGGCGAATGAGGCGCGCGTCGGGCTTTTCGAGc  <  2:1241393/149‑1 (MQ=255)
                                                                                                                                                 cGCCCTTGATCGAAAAGAATATTATCCCGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGc  <  1:984055/150‑1 (MQ=255)
                                                                                                                                                 cGCCCTTGATCGAAAAGAATATTATCCAGCCGACACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGc  <  1:1114717/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTT‑GCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCGAGC  >  NC_000913/3677691‑3677985

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