Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,153,792 T→G 5.9% intergenic (+13/‑224) sibB → / → mdtA sRNA antisense regulator of toxic IbsB protein/multidrug efflux system, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,153,7920TG 5.9% 393.6 / 11.4 129intergenic (+13/‑224)sibB/mdtAsRNA antisense regulator of toxic IbsB protein/multidrug efflux system, subunit A
Reads supporting (aligned to +/- strand):  ref base T (67/54);  new base G (3/5);  total (70/59)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 130 total aligned reads displayed.

GTGAGGGTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAATTCCT  >  NC_000913/2153644‑2153941
                                                                                                                                                    |                                                                                                                                                     
gTGAGGGTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTc                                                                                                                                                      >  2:1609569/1‑150 (MQ=255)
 tGAGGGTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTcc                                                                                                                                                     <  2:319141/150‑1 (MQ=255)
    gggTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGc                                                                                                                                                  <  2:1334297/150‑1 (MQ=255)
       tAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATTAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGa                                                                                                                                               >  1:149555/1‑150 (MQ=255)
       tAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGa                                                                                                                                               >  1:1116028/1‑150 (MQ=255)
          agCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGGGACGCTGaaac                                                                                                                                            <  2:804203/150‑2 (MQ=255)
           gCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGTCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGGGACGCTGaaacg                                                                                                                                           <  1:1063084/150‑3 (MQ=255)
            ccgggTTTCCCCCGCCCTGGTACTCTTTGTAAGCCGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGGGACGCTGaaacgg                                                                                                                                          <  1:1566306/148‑4 (MQ=37)
               ggTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAAtcgtcg                                                                                                                                       >  2:1168262/1‑150 (MQ=255)
                gTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCgg                                                                                                                                      <  2:553803/150‑1 (MQ=255)
                 tttCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGc                                                                                                                                     <  1:1010561/150‑1 (MQ=255)
                    cccccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTa                                                                                                                                  >  1:721163/1‑150 (MQ=255)
                     ccccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGCCCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTAc                                                                                                                                 >  2:1321256/1‑150 (MQ=255)
                         gcgcTGGTATTCTTAGTCAGCGCGGAAGCTTATGACTAAGAGCACCAGGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTAcctcc                                                                                                                             <  2:985061/147‑1 (MQ=25)
                              ggTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGCCCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCg                                                                                                                        >  2:115668/1‑150 (MQ=255)
                              ggTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACATCCTTTCg                                                                                                                        >  2:92259/1‑150 (MQ=255)
                                   tCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCTCCACGATGATGAGTAGCTTCATCATGACCCTTGCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATg                                                                                                                   <  1:330059/150‑1 (MQ=255)
                                     ttAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCa                                                                                                                 <  1:1665441/150‑1 (MQ=255)
                                         tAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCa                                                                                                             >  1:1287421/1‑150 (MQ=255)
                                         tAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCATTTCTCCATGCAAGCa                                                                                                             >  1:1165417/1‑150 (MQ=255)
                                               gggAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTAt                                                                                                       <  2:1539480/150‑1 (MQ=255)
                                               gggAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTAt                                                                                                       <  2:1116028/150‑1 (MQ=255)
                                               gggAAGCTTATGACTAAGAGCACCACGATGATCAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTAt                                                                                                       <  2:149555/150‑1 (MQ=255)
                                                     cTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTaa                                                                                                 >  1:408726/1‑150 (MQ=255)
                                                       tATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTTAcg                                                                                               >  1:132135/1‑150 (MQ=255)
                                                       tATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACATCCTTTCGCCATGCAAGCAGTCTATCGCTTAcg                                                                                               >  1:1312712/1‑150 (MQ=255)
                                                        aTGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATTGCCCCTTCCTCGGGCGGGGCTTTCCCGTTTCAGCGTACCGCTTAAAACGTCTGCTTCCCTCCTTTCGCCATGCAAGCAGTCTATAGCACAcgc                                                                                              >  1:131900/1‑150 (MQ=37)
                                                        aTGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAAcgc                                                                                              >  2:305210/1‑150 (MQ=255)
                                                        aTGACGAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAAcgc                                                                                              >  2:1311607/1‑150 (MQ=255)
                                                          gACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGt                                                                                            >  1:1702623/1‑150 (MQ=255)
                                                           aCTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCGATCCAAGAAGTCTATCGCTAACGCGTa                                                                                           >  2:1321957/1‑150 (MQ=255)
                                                              aaGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGAt                                                                                        <  2:1702623/150‑1 (MQ=255)
                                                                  gaaccacGATGAAGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATaaaa                                                                                    <  2:872904/148‑1 (MQ=255)
                                                                  gCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATaaaa                                                                                    <  2:726299/150‑1 (MQ=255)
                                                                  gCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATaaaa                                                                                    <  2:363398/150‑1 (MQ=255)
                                                                   caccacGATGATGAGGAGCGTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCTCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAAt                                                                                   <  1:288811/150‑1 (MQ=255)
                                                                      cacGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGt                                                                                >  2:1458359/1‑150 (MQ=255)
                                                                      cacGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGt                                                                                <  1:520493/150‑1 (MQ=255)
                                                                       acGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGtt                                                                               >  1:150111/1‑150 (MQ=255)
                                                                        cGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGttt                                                                              <  2:1332067/150‑1 (MQ=255)
                                                                        cGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTc                                                                             >  1:1006653/1‑150 (MQ=255)
                                                                          atgatgAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTcc                                                                            <  1:158693/150‑1 (MQ=255)
                                                                           tgatgaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCt                                                                           <  2:403104/150‑1 (MQ=255)
                                                                             atgaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCtgt                                                                         <  1:1574608/150‑1 (MQ=255)
                                                                             atgaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCtgt                                                                         <  1:1683748/150‑1 (MQ=255)
                                                                               gaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTgtgt                                                                       >  1:898719/1‑150 (MQ=255)
                                                                                aGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGtt                                                                      >  1:899453/1‑150 (MQ=255)
                                                                                 gTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGtta                                                                     <  1:605114/150‑1 (MQ=255)
                                                                                   aGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGttatt                                                                   >  1:351545/1‑150 (MQ=255)
                                                                                       tcatcaTGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTg                                                               >  2:1590764/1‑150 (MQ=255)
                                                                                         atcatGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGa                                                             <  1:73049/150‑1 (MQ=255)
                                                                                          tcatGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGAt                                                            >  2:1685451/1‑150 (MQ=255)
                                                                                          tcatGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGAt                                                            <  1:1250056/150‑1 (MQ=255)
                                                                                              gACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGt                                                        >  2:1267582/1‑150 (MQ=255)
                                                                                               aCCCTTTCCTTATTTCTGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTg                                                       <  2:1251510/150‑1 (MQ=255)
                                                                                               aCCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTg                                                       >  2:1552562/1‑150 (MQ=255)
                                                                                                cccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGc                                                      >  1:987189/1‑150 (MQ=255)
                                                                                                cccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGc                                                      >  1:12532/1‑150 (MQ=255)
                                                                                                cccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGc                                                      >  1:1326208/1‑150 (MQ=255)
                                                                                                     tCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCa                                                 <  1:1267582/150‑1 (MQ=255)
                                                                                                      ccTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCaa                                                >  1:1390010/1‑150 (MQ=255)
                                                                                                       cTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCaaa                                               >  2:211923/1‑150 (MQ=255)
                                                                                                         tatttatGGCCCCTTCCTCGTGAGGGGCGTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAAtg                                             <  2:288482/150‑1 (MQ=255)
                                                                                                          atttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAAtgt                                            >  2:238247/1‑150 (MQ=255)
                                                                                                          atttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAAtgt                                            >  2:1448477/1‑150 (MQ=255)
                                                                                                           tttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAAtgtg                                           <  1:542094/150‑1 (MQ=255)
                                                                                                           tttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAAtgtg                                           >  1:504907/1‑150 (MQ=255)
                                                                                                           tttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGAGTGCTCGCAAAtgtg                                           >  2:984125/1‑150 (MQ=255)
                                                                                                            ttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGc                                          <  2:987189/150‑1 (MQ=255)
                                                                                                             gaaagCCCCTCCCGAGGAAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAAGAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGcc                                         >  1:271210/5‑150 (MQ=21)
                                                                                                                ggCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTTCCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGt                                      >  1:237222/1‑150 (MQ=255)
                                                                                                                 gCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTc                                     >  2:1340108/1‑150 (MQ=255)
                                                                                                                   cccTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCAt                                   <  1:219722/150‑1 (MQ=255)
                                                                                                                    ccTTCCTCGGGACGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTCCCCGTCAtt                                  >  1:246178/1‑150 (MQ=255)
                                                                                                                     cTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAATTAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTc                                 >  2:267785/1‑150 (MQ=255)
                                                                                                                      ttCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCa                                <  2:23844/150‑1 (MQ=255)
                                                                                                                           cGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATTAAATTGTTTCCTGTGTTTTTACTGGATGCGTGGTCGCAAATGTGCCCGTCTTTCAGaggt                           >  1:1694254/1‑147 (MQ=21)
                                                                                                                             ggAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGAtt                         >  1:1035638/1‑150 (MQ=255)
                                                                                                                               aGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAATCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTcc                       <  1:115668/150‑1 (MQ=255)
                                                                                                                               aGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTcc                       <  2:919427/150‑1 (MQ=255)
                                                                                                                               aGGGGCTTTCCCGTTTCAGCGTCACGCTGAAATCGTCTGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTcc                       <  1:92259/150‑1 (MQ=255)
                                                                                                                                ggggCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCa                      >  2:474602/1‑150 (MQ=255)
                                                                                                                                  ggCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGa                    >  1:776574/1‑150 (MQ=255)
                                                                                                                                  ggCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGa                    >  1:983470/1‑150 (MQ=255)
                                                                                                                                  ggCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGa                    >  1:574673/1‑150 (MQ=255)
                                                                                                                                  ggCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGGTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGa                    >  1:437412/1‑150 (MQ=255)
                                                                                                                                   gCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGAc                   >  1:1619072/1‑150 (MQ=255)
                                                                                                                                   gCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGAc                   <  1:1685451/150‑1 (MQ=255)
                                                                                                                                        cccGTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCCTGAAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGAGTTATTAATGGATGCCTTCTCTCAAATGTGAACGTCATTCAGAAGATTACAGAAAgagt              >  2:91426/1‑150 (MQ=17)
                                                                                                                                         cgatttCAGCGGGACGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGtt             >  1:976440/4‑150 (MQ=255)
                                                                                                                                         ccGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAGCGCGTAGATAAAATAGTTTCCTGTGTTATTTCTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGCCAGTGtt             >  1:1643128/1‑150 (MQ=255)
                                                                                                                                         ccGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCAATTCAAGCTGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGtt             <  2:1165417/150‑1 (MQ=255)
                                                                                                                                          gatttCAGCGGGACGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGttt            >  1:1274687/3‑150 (MQ=255)
                                                                                                                                           gTTTCAGCGTCCCGCTGAAATCTTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAAt                                              <  1:896471/115‑1 (MQ=255)
                                                                                                                                           gTTTCAGCGTCCCGCTGAAATCATCGGCTTACCTCCTTTCGCCATGAAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTAATGGATGCGTGCTCGCAAAt                                              >  2:896471/1‑115 (MQ=255)
                                                                                                                                            tttCAGCGGGACGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCa          >  1:1371807/1‑150 (MQ=255)
                                                                                                                                             ttCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCATCATTCAGACGATTCCAGACAGTGTTTCAt         <  2:382596/150‑1 (MQ=255)
                                                                                                                                                cgCGTCCCGCTGAAATCTTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAAt      <  2:44749/149‑1 (MQ=255)
                                                                                                                                                 gCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAAtt     <  2:12532/150‑1 (MQ=255)
                                                                                                                                                    tCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTAATTACTGGATGCGTGCTCGCAGATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAATtcct  <  2:1568989/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGAGGGTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAATTCCT  >  NC_000913/2153644‑2153941

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: