Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,044,729 T→A 8.2% intergenic (+105/‑209) 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,7290TA 8.2% 265.1 / 12.0 98intergenic (+105/‑209)asnT/yeeJtRNA‑Asn/putative adhesin
Reads supporting (aligned to +/- strand):  ref base T (43/47);  new base A (6/2);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

TGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCAGTATT  >  NC_000913/2044581‑2044876
                                                                                                                                                    |                                                                                                                                                   
tGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATg                                                                                                                                                    <  1:1367535/150‑1 (MQ=255)
tGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATg                                                                                                                                                    >  1:1513339/1‑150 (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:716529/1‑150 (MQ=255)
                                                          gCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGa                                                                                          >  1:1102237/1‑150 (MQ=255)
                                                          gCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGa                                                                                          <  2:1626738/150‑1 (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                                                                              <  2:244175/150‑1 (MQ=255)
                                                                      cAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCagtag                                                                              <  1:860767/150‑1 (MQ=255)
                                                                      cAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCagtag                                                                              <  1:155692/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:296993/150‑1 (MQ=255)
                                                                                                       tAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGcc                                             <  1:716529/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             >  2:157549/1‑150 (MQ=255)
                                                                                                                                       cccTCAACCATAAAGGTCGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTg             >  1:911001/1‑150 (MQ=255)
                                                                                                                                            gACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATc        <  1:644568/150‑1 (MQ=255)
                                                                                                                                             aCCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCATCGATTATAAAACTAGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATct       >  2:292486/1‑149 (MQ=255)
                                                                                                                                                  aaTGGTTGAGGGCGTTTTGCTATGCACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAACCTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCattatt  <  2:258718/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCAGTATT  >  NC_000913/2044581‑2044876

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: