Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,044,727 T→A 8.3% intergenic (+103/‑211) 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,7270TA 8.3% 242.6 / 11.6 97intergenic (+103/‑211)asnT/yeeJtRNA‑Asn/putative adhesin
Reads supporting (aligned to +/- strand):  ref base T (43/46);  new base A (6/2);  total (49/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 = 1.19e-01

GACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCA  >  NC_000913/2044578‑2044871
                                                                                                                                                     |                                                                                                                                                
gACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCttt                                                                                                                                                  >  2:443432/1‑150 (MQ=255)
   tGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATg                                                                                                                                               <  1:1367535/150‑1 (MQ=255)
   tGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATg                                                                                                                                               >  1:1513339/1‑150 (MQ=255)
    gTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATgg                                                                                                                                              >  1:561554/1‑150 (MQ=255)
    gTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATgg                                                                                                                                              <  1:1124501/150‑1 (MQ=255)
     ttAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGt                                                                                                                                             <  2:1535214/150‑1 (MQ=255)
       aaTCCGTATGTCACTGGTTCGAGTGCAGTCAGAGGATCCAAAATCTAAAAATTCGCTTTTTTAGCGGAATGTCACTGACCATAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTg                                                                                                                                           <  2:168998/150‑1 (MQ=39)
       aaTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTg                                                                                                                                           >  1:1693080/1‑150 (MQ=255)
       aaTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTg                                                                                                                                           <  2:157178/150‑1 (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                                                                                            >  1:450288/1‑150 (MQ=255)
                                                      aaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATaa                                                                                            <  2:775697/150‑1 (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                                                                                     >  1:1102237/1‑150 (MQ=255)
                                                             gCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGa                                                                                     >  2:716529/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                                                                         <  2:244175/150‑1 (MQ=255)
                                                                         cAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCagtag                                                                         <  1:155692/150‑1 (MQ=255)
                                                                         cAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCagtag                                                                         <  1:860767/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         >  2:918983/1‑150 (MQ=255)
                                                                                                                                         gcCCTCAACCATAAAGGTCGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATAtt         >  1:550120/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)
                                                                                                                                                aCCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCATCGATTATAAAACTAGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATct  >  2:292486/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCA  >  NC_000913/2044578‑2044871

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: