Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,334 C→T 9.5% intergenic (‑477/+18) lysO ← / ← aqpZ L‑lysine exporter/water channel AqpZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,3340CT 9.5% 246.2 / 12.1 84intergenic (‑477/+18)lysO/aqpZL‑lysine exporter/water channel AqpZ
Reads supporting (aligned to +/- strand):  ref base C (37/39);  new base T (3/5);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGC  >  NC_000913/915189‑915482
                                                                                                                                                 |                                                                                                                                                    
ccGGGGGGGCTGGGGGGGGCGGCTATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAgg                                                                                                                                                  <  2:580196/150‑1 (MQ=255)
     ggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTc                                                                                                                                             >  2:214914/1‑150 (MQ=255)
     ggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTc                                                                                                                                             >  1:1565573/1‑150 (MQ=255)
        gctggggggggggGCTATTGCCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGGCGTAATGTCAGGCCtttttt                                                                                                                                          <  1:1040143/147‑5 (MQ=37)
        gctggggggggCGGCTTTGGGCCCCCCCGGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGGCGTAATGTCAGGCCtttttt                                                                                                                                          <  2:1718858/147‑5 (MQ=255)
        gccggggggggggGCTTTTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGGCGTAATGTCAGGCCtttttt                                                                                                                                          <  2:1941628/147‑5 (MQ=37)
        gCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACACCAGGCCTTCTGc                                                                                                                                          >  1:950802/1‑150 (MQ=255)
         gtggggggggggTCTATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                                                                                                                                         <  2:1676315/148‑1 (MQ=255)
         cAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                                                                                                                                         >  2:1714929/1‑150 (MQ=255)
           ggggggggggGCTTTGGGCCGCCCCGGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                                                                                                                                       <  2:464273/150‑1 (MQ=255)
           gggggCGGGGGGTTTGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                                                                                                                                       <  2:1927943/150‑1 (MQ=39)
           gggggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGGCGTAATGTCAGGCCtttttttaa                                                                                                                                       <  2:1724507/150‑8 (MQ=37)
           gggggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACACCAGGCCTTCTGCGtt                                                                                                                                       <  2:950802/150‑1 (MQ=255)
            ggggCGGGGGGTATTGGCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                                                                                                                                      <  2:1911922/150‑1 (MQ=39)
              ggcggcggcTTTTGGCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat                                                                                                                                    <  1:1268590/150‑1 (MQ=255)
              ggcggcggcGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat                                                                                                                                    <  1:1725381/150‑1 (MQ=255)
                   cggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc                                                                                                                               <  2:679204/150‑1 (MQ=255)
                    ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                                                                                                                              <  2:4420/150‑1 (MQ=255)
                        tttGGCCGCCCCGGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt                                                                                                                          <  2:1709680/149‑1 (MQ=255)
                           ggCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtttt                                                                                                                       >  1:78493/1‑150 (MQ=255)
                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc                                                                                                                     <  2:327296/150‑1 (MQ=255)
                                  cccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag                                                                                                                <  2:943117/150‑1 (MQ=255)
                                    cgccgcgcTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAggg                                                                                                              <  1:460346/147‑1 (MQ=255)
                                    cTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACCACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAggg                                                                                                              >  2:1802373/1‑150 (MQ=255)
                                       cgcgcTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTcc                                                                                                           <  1:1443630/150‑1 (MQ=255)
                                         cgcTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCgg                                                                                                         <  2:1147203/150‑1 (MQ=255)
                                         cgcTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCgg                                                                                                         >  2:1361874/1‑150 (MQ=255)
                                             ccTTGCGCCCGTGGCAATATGTTGCTTACCTCATGAAAGGAGCGCAACAAGATGATGAATCACCATATAACCCCATCTTAAAAAAAGGCCTGCCATTACGCCAGGCCTTCTCCGTTAATTAATCACGCTTTTGCAGCAGGGAACGGTaaa                                                                                                     >  2:165987/1‑150 (MQ=255)
                                             ccTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTaaa                                                                                                     >  1:620857/1‑150 (MQ=255)
                                              cTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAAt                                                                                                    >  2:1583367/1‑150 (MQ=255)
                                               ttGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATc                                                                                                   <  1:787549/150‑1 (MQ=255)
                                                    ccAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGacc                                                                                              >  1:1936485/1‑150 (MQ=255)
                                                        tGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCg                                                                                          >  1:1897672/1‑150 (MQ=255)
                                                          gCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGGATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGAt                                                                                        >  1:913219/1‑150 (MQ=255)
                                                          gCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGAt                                                                                        <  1:1099952/150‑1 (MQ=255)
                                                           cAATATGTTGCTTAGCTAATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAAAATCTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATAACGCTTTTCCAGAAGGGTCCGGTAAATCAGACCACCGATa                                                                                       >  2:245149/1‑150 (MQ=255)
                                                             atatGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAAt                                                                                     <  2:1565573/150‑1 (MQ=255)
                                                               atGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgc                                                                                   <  1:1620326/150‑1 (MQ=255)
                                                               atGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgc                                                                                   <  1:517643/150‑1 (MQ=255)
                                                                 gTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgccg                                                                                 >  2:703924/1‑150 (MQ=255)
                                                                   cccTTACCTCATCAACGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgccgcc                                                                               <  2:839664/148‑1 (MQ=255)
                                                                     cTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgacgccag                                                                             >  1:839664/1‑148 (MQ=255)
                                                                       tAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACa                                                                           >  2:388742/1‑150 (MQ=255)
                                                                         gCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAGAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAAt                                                                         >  1:983285/1‑150 (MQ=255)
                                                                         gCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAAt                                                                         <  1:214914/150‑1 (MQ=255)
                                                                          cTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGGCCGGTAAATCAGACAAGCGATAATGCCGCCGACAAtt                                                                        >  2:182903/1‑150 (MQ=255)
                                                                            cATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTgg                                                                      >  1:822033/1‑150 (MQ=255)
                                                                             aTGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGc                                                                     <  1:1583367/150‑1 (MQ=255)
                                                                             aTGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACAACCGATAATGCCGCCGACAATTGGc                                                                     >  1:1502213/1‑150 (MQ=255)
                                                                               gAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGcac                                                                   >  2:1864981/1‑150 (MQ=255)
                                                                               gAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGcac                                                                   >  2:1915057/1‑150 (MQ=255)
                                                                               gAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATCCCGCCGACAATTGGcac                                                                   >  2:1868869/1‑150 (MQ=255)
                                                                                   ggAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCaccacc                                                               >  2:184362/1‑150 (MQ=255)
                                                                                    gAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCAccc                                                              <  1:388742/150‑1 (MQ=255)
                                                                                      gcgcAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCag                                                            >  1:291173/1‑150 (MQ=255)
                                                                                       cgcAACAAGAAGATGAATCAACATATAACAACATCTTAAAAAAACGCCTGACAATACGCCAGGCCTTCTGCCATAATTAATCACGCTTTTCCAGCAGAGTCCGGAAAATAAGACAACCGATAATGCCGCAGACAAATGGCACAACACaga                                                           >  1:671162/1‑150 (MQ=255)
                                                                                         caacaaGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACa                                                                           >  1:1202984/1‑132 (MQ=255)
                                                                                         caacaaGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACa                                                                           <  2:1202984/132‑1 (MQ=255)
                                                                                         caacaaGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGGACCACCCagaag                                                         >  2:1981704/1‑150 (MQ=255)
                                                                                           acaaGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAAc                                                      <  1:1908686/150‑1 (MQ=255)
                                                                                              aGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCATCAGGGTCCGTTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACca                                                    <  1:245149/150‑1 (MQ=255)
                                                                                                  gatgaATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgt                                                >  1:1144282/1‑150 (MQ=255)
                                                                                                  gatgaATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgt                                                <  1:1714929/150‑1 (MQ=255)
                                                                                                    tgaATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgttg                                              <  1:703924/150‑1 (MQ=255)
                                                                                                    tgaATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgttg                                              <  2:1923040/150‑1 (MQ=255)
                                                                                                       aTCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTc                                           >  2:1262917/1‑150 (MQ=255)
                                                                                                       aTCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACAACCGATAAAGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTc                                           >  2:1257422/1‑150 (MQ=255)
                                                                                                            cATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATg                                      <  1:1077727/150‑1 (MQ=255)
                                                                                                               ataACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGccc                                   <  1:1027686/150‑1 (MQ=255)
                                                                                                                 aacaacATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAg                                 >  2:1022447/1‑150 (MQ=255)
                                                                                                                  acaacaTCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgc                                >  2:194972/1‑150 (MQ=255)
                                                                                                                      caTCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgccgcc                            <  2:1984505/150‑1 (MQ=255)
                                                                                                                           acgcagaaGGCCTGGCGTAATGTCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTg                         >  2:1787057/7‑148 (MQ=255)
                                                                                                                           acgcagaaGGCCTGGCGTAATGTCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTg                         <  1:1787057/142‑1 (MQ=255)
                                                                                                                               aaaaGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGAt                   <  2:620857/150‑1 (MQ=255)
                                                                                                                                aaaGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATa                  <  1:184362/150‑1 (MQ=255)
                                                                                                                                  aGGCCTGGCGTAATGTCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGc                >  1:882877/1‑150 (MQ=255)
                                                                                                                                  aGGCCTGGCGTAATGTCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGc                >  1:22658/1‑150 (MQ=255)
                                                                                                                                       tGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACcg           >  1:1214400/1‑150 (MQ=255)
                                                                                                                                        gACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACAACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACcgc          >  2:1449109/1‑150 (MQ=255)
                                                                                                                                          cATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCgg        <  1:1915057/150‑1 (MQ=255)
                                                                                                                                             tACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGtgc     <  1:470508/150‑1 (MQ=255)
                                                                                                                                             tACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGtgc     >  1:1372204/1‑150 (MQ=255)
                                                                                                                                                gCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGtgctgc  <  2:983285/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGC  >  NC_000913/915189‑915482

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