Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,044,721 G→A 8.2% intergenic (+97/‑217) asnT → / → yeeJ tRNA‑Asn/putative adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,044,7210GA 8.2% 288.6 / 11.1 99intergenic (+97/‑217)asnT/yeeJtRNA‑Asn/putative adhesin
Reads supporting (aligned to +/- strand):  ref base G (43/46);  new base A (6/2);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-02

CGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATC  >  NC_000913/2044576‑2044870
                                                                                                                                                 |                                                                                                                                                     
cggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGCCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCt                                                                                                                                                   <  1:1692803/150‑1 (MQ=255)
 ggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCTACCtt                                                                                                                                                  >  2:937383/1‑150 (MQ=255)
 ggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGAGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATTGCGGGTTTTTAACATCTTAAGCGCCCTCGACCtt                                                                                                                                                  >  2:72735/1‑150 (MQ=255)
  gACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCttt                                                                                                                                                 >  2:443432/1‑150 (MQ=255)
     tGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATg                                                                                                                                              >  1:1513339/1‑150 (MQ=255)
     tGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATg                                                                                                                                              <  1:1367535/150‑1 (MQ=255)
      gTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATgg                                                                                                                                             <  1:1124501/150‑1 (MQ=255)
      gTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATgg                                                                                                                                             >  1:561554/1‑150 (MQ=255)
       ttAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGt                                                                                                                                            <  2:1535214/150‑1 (MQ=255)
         aaTCCGTATGTCACTGGTTCGAGTGCAGTCAGAGGATCCAAAATCTAAAAATTCGCTTTTTTAGCGGAATGTCACTGACCATAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTg                                                                                                                                          <  2:168998/150‑1 (MQ=39)
         aaTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTg                                                                                                                                          <  2:157178/150‑1 (MQ=255)
         aaTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTg                                                                                                                                          >  1:1693080/1‑150 (MQ=255)
         aaTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGAAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTg                                                                                                                                          >  1:1332517/1‑150 (MQ=255)
          aTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTATAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAAAATCTTAAGCGCACTCTACCTTTATGGTTGa                                                                                                                                         >  1:1130710/1‑150 (MQ=255)
           tcagTATGTCACTGGTTGGAGTCCAGTCAGAGGAGCCAAATTCTAAAATTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACTTTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAg                                                                                                                                        <  1:1274846/147‑1 (MQ=255)
           tCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGGTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAg                                                                                                                                        <  1:1291124/150‑1 (MQ=255)
             cGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCAACCATAAAGGTCGAggg                                                                                                                                      <  1:527340/150‑1 (MQ=38)
             cGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCATAAACGCCCTCAACCATAAAGGTCGAggg                                                                                                                                      <  2:1581683/150‑1 (MQ=38)
               tATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGCCCTTAGTTGACCATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGTGCg                                                                                                                                    >  1:1565487/1‑150 (MQ=255)
                  gTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGttt                                                                                                                                 >  2:1661723/1‑150 (MQ=255)
                    cACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCATAGTTGAACATTGATTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGTCGTTTTg                                                                                                                               >  2:1236862/1‑150 (MQ=255)
                     aCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGc                                                                                                                              <  2:9654/150‑1 (MQ=255)
                     aCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGc                                                                                                                              <  1:1227879/150‑1 (MQ=255)
                       tGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTa                                                                                                                            <  2:1315087/150‑1 (MQ=255)
                        ggTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTAt                                                                                                                           >  1:998363/1‑150 (MQ=255)
                        ggTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGATTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGACCTCGACCTTTATGGTTGAGGGCTTTTTGatat                                                                                                                           >  2:98637/1‑150 (MQ=255)
                         gTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATg                                                                                                                          >  1:1166944/1‑150 (MQ=255)
                             gAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGt                                                                                                                                   >  2:586703/1‑137 (MQ=255)
                             gAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGt                                                                                                                                   <  1:586703/137‑1 (MQ=255)
                                  cAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATc                                                                                                                 >  2:1046405/1‑150 (MQ=255)
                                   aGTCAGCGGAGCCAAATTGTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCa                                                                                                                <  2:836057/150‑1 (MQ=255)
                                    gTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAc                                                                                                               >  2:909656/1‑150 (MQ=255)
                                     tCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAcc                                                                                                              <  1:1318596/150‑1 (MQ=255)
                                       agagGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCAt                                                                                                            >  2:795011/1‑150 (MQ=255)
                                        gaggagCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCAtt                                                                                                           >  2:296993/1‑150 (MQ=255)
                                           gagCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTc                                                                                                        >  2:1697891/1‑150 (MQ=255)
                                           gagCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTc                                                                                                        >  1:1702075/1‑150 (MQ=255)
                                            agCCAAATTCTAAAAATTCGCTTGTTTAGCGCAATGTCACTGACCTTAGTTGAACTTTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTcc                                                                                                       <  1:25865/150‑1 (MQ=255)
                                            agCCAAAATCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTcc                                                                                                       <  1:641551/150‑1 (MQ=255)
                                               cAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCt                                                                                                    <  1:592835/150‑1 (MQ=255)
                                                  aTTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGa                                                                                                 >  2:644568/1‑150 (MQ=255)
                                                   ttCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGTCCTTAGTTAAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGAt                                                                                                <  2:1048996/150‑1 (MQ=255)
                                                    tCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGAtt                                                                                               >  2:939588/1‑150 (MQ=255)
                                                    gcTAAAAAAGCGAATTTTTAGCGCAATGTCACTGACCTTAGTTGATCATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCATCGAtt                                                                                               >  1:258718/2‑150 (MQ=255)
                                                      tAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATtat                                                                                             <  1:892203/150‑1 (MQ=255)
                                                        aaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAATCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATaa                                                                                           <  2:140753/150‑1 (MQ=255)
                                                        aaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATaa                                                                                           <  2:775697/150‑1 (MQ=255)
                                                        aaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATaa                                                                                           >  1:450288/1‑150 (MQ=255)
                                                           aTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAAc                                                                                        <  2:1488333/150‑1 (MQ=255)
                                                            ttCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACt                                                                                       <  1:1661723/150‑1 (MQ=255)
                                                             tCGCTTATTTAGCGCAATGTCACTGACCTTAGTTGACCATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCACTCGCCCTTTATGGTTGAGGGCGTGTTGCTATGAACGCCATCCCCATTTTCCCCTCGATTATAAAACtt                                                                                      >  2:290771/1‑150 (MQ=255)
                                                               gCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGa                                                                                    <  2:1626738/150‑1 (MQ=255)
                                                               gCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGa                                                                                    >  2:716529/1‑150 (MQ=255)
                                                               gCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGa                                                                                    >  1:1102237/1‑150 (MQ=255)
                                                                gttttttAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAg                                                                                   <  1:243641/149‑1 (MQ=255)
                                                                cTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAg                                                                                   <  2:211298/150‑1 (MQ=255)
                                                                cTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAg                                                                                   <  1:986414/150‑1 (MQ=255)
                                                                cTTTTTTAGCGAAATGTCACTGACCTTAGTTGAAAATTTTTTTTTAACGGATAGCGGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTGGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAg                                                                                   <  1:671122/150‑1 (MQ=255)
                                                                  tttttAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTATTAACTGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATTGTTGAGGTCTATTAGCTATGAACGCCATCACCATTTACCACTCGATTATAAAACTTGAgat                                                                                 >  1:225305/1‑148 (MQ=255)
                                                                    tttAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTAt                                                                               >  2:1083018/1‑150 (MQ=255)
                                                                    tttAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTGAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATTAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTAt                                                                               >  1:1055792/1‑150 (MQ=255)
                                                                      tAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGCCCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTc                                                                             >  1:69548/1‑150 (MQ=255)
                                                                       aGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCa                                                                            >  1:734250/1‑150 (MQ=255)
                                                                        gcgcAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACACCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCag                                                                           <  2:742676/150‑1 (MQ=255)
                                                                           cAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCagtag                                                                        <  1:860767/150‑1 (MQ=255)
                                                                           cAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCagtag                                                                        <  1:155692/150‑1 (MQ=255)
                                                                           cAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCagtag                                                                        <  2:244175/150‑1 (MQ=255)
                                                                           cAATGTCACTGACCTTAATTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCagtag                                                                        <  1:359446/150‑1 (MQ=255)
                                                                           aaaTGTCACTGACCTTATTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCagtag                                                                        <  2:1332517/149‑1 (MQ=255)
                                                                                tCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCATTAGTCTcc                                                                   >  1:302993/1‑150 (MQ=255)
                                                                                     gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCctct                                                              <  2:1693080/150‑1 (MQ=255)
                                                                                            gTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTTTGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACt                                                       <  1:1700865/150‑1 (MQ=255)
                                                                                            gTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACt                                                       <  2:1474253/150‑1 (MQ=255)
                                                                                               gAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTcac                                                    >  2:1488909/1‑150 (MQ=255)
                                                                                                aaCATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCCCCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAGCTcacc                                                   >  2:1725176/1‑149 (MQ=255)
                                                                                                   aTTGTTTTTTAACGGATAGAGGGTTTTGAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACAccc                                                <  2:267348/150‑1 (MQ=255)
                                                                                                        tttttAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAAc                                           >  1:1091497/1‑150 (MQ=255)
                                                                                                         ttttAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCGAAACt                                          >  2:700129/1‑150 (MQ=255)
                                                                                                         ttttAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACt                                          <  2:561554/150‑1 (MQ=255)
                                                                                                          tttAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTg                                         >  1:37462/1‑150 (MQ=255)
                                                                                                            tAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGcc                                       <  1:716529/150‑1 (MQ=255)
                                                                                                            tAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGcc                                       <  1:296993/150‑1 (MQ=255)
                                                                                                             aaCGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCt                                      <  1:909656/150‑1 (MQ=255)
                                                                                                                     gCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTTTTCAGTAGTCTCCCCTCTTGCAACTCACGCCCAAAAGTGCCTAACGaaaa                              >  1:253841/1‑150 (MQ=255)
                                                                                                                      cGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGTGTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATCCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAg                             <  1:882155/150‑1 (MQ=255)
                                                                                                                        ggTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCCACTCACACCCTAAACTGCCTTACGAATAGtt                           >  1:1556381/1‑150 (MQ=255)
                                                                                                                         gTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGtta                          <  1:922930/150‑1 (MQ=255)
                                                                                                                         gTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGATATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGtta                          <  1:263858/150‑1 (MQ=255)
                                                                                                                          tttttAACATCTTAAGCGCCCTCGCCCTTTATGGTTGAGGGCGTTTTGGTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCGTGCAACTCACCCCCAAAACTGCCTAACGAAAAGttat                         >  1:153048/1‑150 (MQ=255)
                                                                                                                                   agcaaaacgcCCTCAACCATAAAGGTCGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAACACTGCCTAACGAAAAGTTATTAATTTTCa                >  2:476880/8‑150 (MQ=255)
                                                                                                                                    cTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATAAATTTTCaa               >  1:161561/1‑150 (MQ=255)
                                                                                                                                       aaGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGACCGCCATCACCATTTTCCCCTCGATTATGAATCTTGAGTTATTCAGTAGTCTCCCGTCTTGCAACTCACACCCAAAACTGCCTAACGAATAGTTATTAATTTTCAATCa            >  1:257326/1‑150 (MQ=255)
                                                                                                                                       aaGCGCCCTCGAACTTTATTGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCa            <  2:225305/150‑1 (MQ=255)
                                                                                                                                           gcCCTCAACCATAAAGGTCGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATAtt        >  1:550120/1‑150 (MQ=255)
                                                                                                                                           gcCCTCAACCATAAAGGTCGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATAtt        >  2:918983/1‑150 (MQ=255)
                                                                                                                                           gcCCTCAACCATAAAGGTCGAGGGCGCTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATAtt        >  2:1661104/1‑150 (MQ=255)
                                                                                                                                            cccTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTg       >  1:388778/1‑150 (MQ=255)
                                                                                                                                            cccTCAACCATAAAGGTCGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTg       >  1:911001/1‑150 (MQ=255)
                                                                                                                                            cccTCAACCATAAAGGTCGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTg       >  2:157549/1‑150 (MQ=255)
                                                                                                                                                 gACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATc  <  1:644568/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATC  >  NC_000913/2044576‑2044870

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