Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,424,646 C→A 6.7% noncoding (2138/2904 nt) rrlD ← 23S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,424,6460CA 6.7% 432.2 / 7.5 164noncoding (2138/2904 nt)rrlD23S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base C (78/71);  new base A (2/11);  total (80/84)
Fisher's exact test for biased strand distribution p-value = 1.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

CTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGACAG  >  NC_000913/3424505‑3424794
                                                                                                                                             |                                                                                                                                                     
ctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCt                                                                                                                                                 >  2:63231/1‑148 (MQ=11)
ctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATAAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCt                                                                                                                                                 >  2:115435/1‑148 (MQ=11)
   cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTaa                                                                                                                                               <  1:115435/147‑1 (MQ=32)
   cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTaa                                                                                                                                               <  3:140270/147‑1 (MQ=32)
   cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGc                                                                                                                                            <  1:59327/149‑1 (MQ=32)
     ttACTCTTTAGGAGGAGACCGCACCAGTCAAACTACCCACAAGACACTGTCCGCAACACGGATTACGGGTCAACTTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCAAACTTCTAAGCCt                                                                                                                                          >  3:237075/1‑149 (MQ=11)
        ctctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTcc                                                                                                                                        <  5:152770/148‑1 (MQ=32)
         tctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCa                                                                                                                                      >  2:58609/1‑149 (MQ=32)
           tttAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAcc                                                                                                                                    >  5:328295/1‑149 (MQ=32)
             tAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAc                                                                                                                                     <  5:258361/146‑1 (MQ=32)
             tAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAc                                                                                                                                     <  3:97455/146‑1 (MQ=32)
              aggaggAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTAt                                                                                                                                 >  4:91411/1‑149 (MQ=32)
                gaggagACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATcc                                                                                                                               >  1:128993/1‑149 (MQ=32)
                gaggagACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTACCACATATcc                                                                                                                               >  3:266845/1‑149 (MQ=17)
                 acgagaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGCTCAACGTTAGAACAACAAACATAAAAGGCTGGTATTTCAAGGTAGGCTCCATGCACACTGGCCTCCACACATCTAAGCCTCCCACCAATcc                                                                                                                               >  2:285272/3‑148 (MQ=11)
                    agaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTaca                                                                                                                           <  5:166185/149‑1 (MQ=32)
                    agaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTaca                                                                                                                           <  7:308970/149‑1 (MQ=32)
                      aCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCCGATTACGCGTCAACGTTAGAACAACAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCACACTGGCGTCCACACTTCTAAGCCTCCCAACTATCCTAcaca                                                                                                                         >  4:61920/1‑149 (MQ=11)
                           cccAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAg                                                                                                                    <  1:58609/149‑1 (MQ=32)
                            ccAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGGGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAgg                                                                                                                   <  1:250404/149‑1 (MQ=17)
                              aGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCt                                                                                                                              >  2:327426/1‑136 (MQ=32)
                               gTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCt                                                                                                                 >  8:81780/1‑148 (MQ=32)
                                tCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTc                                                                                                                >  7:277795/1‑148 (MQ=32)
                                tCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTc                                                                                                                >  6:143407/1‑148 (MQ=32)
                                   aaCTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATg                                                                                                            >  4:192063/1‑149 (MQ=32)
                                      tACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCCTGCAGCCTGGCGGCCACACTTCTAAGCCTCCCACCTATCCTACCCATCAAGGCTCAATGtt                                                                                                          >  3:267448/1‑148 (MQ=255)
                                      tACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGtt                                                                                                          >  7:237204/1‑148 (MQ=32)
                                      tACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGtt                                                                                                          >  7:74400/1‑148 (MQ=32)
                                       aCCCCCCCGACCCTGTCCGCAACCCGGATTACGGTTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGAATTTTTTTTTAATTTTTATGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCa                                                                                                        <  8:122500/149‑1 (MQ=2)
                                        cccACCAGACACTGTCCGCAACCCGGCTTACCGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAg                                                                                                       <  1:84351/149‑1 (MQ=17)
                                        cccACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCa                                                                                                        <  1:236981/148‑1 (MQ=32)
                                        cccACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAg                                                                                                       >  8:226037/1‑149 (MQ=32)
                                        cccACAAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTAAg                                                                                                       >  2:23921/1‑149 (MQ=17)
                                         ccaccaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCAAGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAgt                                                                                                      <  7:223939/149‑1 (MQ=17)
                                           accaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAgtgt                                                                                                    >  2:275321/1‑149 (MQ=32)
                                             caGACACTGTCCGCAACCCGGATTACGGGGCAACGATAGAACATAAAACATTAAAGGGTGGTATTTCAAGGTCGGCGCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCACAAGGTTCATTGTc                                                                                                   >  5:227600/1‑148 (MQ=11)
                                              aGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTc                                                                                                   >  6:8034/1‑147 (MQ=32)
                                                acacTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGc                                                                                               >  5:145369/1‑149 (MQ=32)
                                                  acTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAATGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCta                                                                                             <  1:98878/149‑1 (MQ=17)
                                                  acTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATTCATACTTTTTTTATAATTTTTATGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCta                                                                                             <  8:115443/149‑1 (MQ=11)
                                                   cTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCACACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCtat                                                                                            <  1:294316/149‑1 (MQ=17)
                                                      tCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCATAATTGTTTTTAAATTTTTATGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCtata                                                                                           <  2:260337/147‑1 (MQ=11)
                                                         gCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTaa                                                                                                                                               <  3:39471/93‑1 (MQ=33)
                                                           aaCCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCATCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAgg                                                                                    >  5:262085/1‑149 (MQ=17)
                                                            aCCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGt                                                                                   >  8:258658/1‑149 (MQ=32)
                                                             cccGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGtt                                                                                  >  1:266479/1‑149 (MQ=32)
                                                              ccGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGtt                                                                                  >  1:139830/1‑148 (MQ=32)
                                                              ccGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGtt                                                                                  >  8:127590/1‑148 (MQ=32)
                                                                 gATTACGGGTCAACGTTATAACATCAAACATTAAAGGGTGGTATTTCAAGTTCGGCTCCCTGCAGACTGGCGTCCAGACTTCTATGCCTCCCACCTATCGTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACg                                                                              <  1:285272/149‑1 (MQ=11)
                                                                          tCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTc                                                                     <  7:81780/149‑1 (MQ=32)
                                                                              cGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGCCTCAATGTTCAGAGTCAAGCTATAGTAAAGGTTCACGCGGACTTTACGTc                                                                 >  1:194708/1‑149 (MQ=11)
                                                                                ttAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTATTTACGTCtt                                                               >  3:19151/1‑149 (MQ=17)
                                                                                   gAAAATCAAACATTAAATGGTGGTATTTCAATGTCGGCTCCATGCAGATTTTTTTTTTAATTTTTATGCGTCCAACGATTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGc                                                             <  4:266845/148‑1 (MQ=2)
                                                                                     aCATCAAACATTAAAGGGGGGTATTTCACGGTCGGCTCCATGCAGACTGGCGGCCACACTTCTAAGCCTCCCACCTCTCCTACACCTCAAGGCTCAATGTTCATTGTCAAGCTATAGTAAAGGTTCACGGGGGCTTTCCGTCTTGCcgc                                                          <  7:82622/149‑1 (MQ=11)
                                                                                      cagcAAACATTAACGGGTGGTATTTCAAGGTCGTCTCCATGCAGATTTTTTTTTTAATTATAATGCGACCAACGATTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCcgcg                                                         <  4:176931/146‑1 (MQ=0)
                                                                                       aTCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCgg                                                        <  1:243872/149‑1 (MQ=32)
                                                                                       aTCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCgg                                                        <  5:136021/149‑1 (MQ=32)
                                                                                        tCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCggg                                                       >  4:182489/1‑149 (MQ=32)
                                                                                        tCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCggg                                                       >  4:102993/1‑149 (MQ=32)
                                                                                        tCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATAAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCAGCggg                                                       >  3:115068/1‑149 (MQ=17)
                                                                                         cAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCggg                                                       >  6:326732/1‑148 (MQ=32)
                                                                                         cAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTAACGGGGTATTTCCGTCTTGCCGCggg                                                       >  1:39670/1‑148 (MQ=17)
                                                                                           aaCATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTa                                                     >  7:132354/1‑148 (MQ=32)
                                                                                              aTTAAAGGGTTGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACt                                                 <  5:89233/149‑1 (MQ=17)
                                                                                              aTTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACt                                                 >  8:35700/1‑149 (MQ=32)
                                                                                               ttAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTg                                                >  6:134902/1‑149 (MQ=32)
                                                                                                 aaaGGGTGGTATTTCAAGGTCGGCTCCATGCAGAATGGCGTCAACACTTCTAATCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGc                                               <  6:227600/148‑1 (MQ=12)
                                                                                                  aaGGGTGGTATTTCAAGGTCGGCTCCATGCAGATTGGTTTTAAAACTTCTAAGCCTCCAACCTTTGCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCa                                              <  6:192187/148‑1 (MQ=11)
                                                                                                  aaGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCAt                                             >  8:280627/1‑149 (MQ=32)
                                                                                                     ggtgttATTTCAAGGTCGGCTCCATGCAGACTGGCGTACAAACTTATAAGCCTCCCACGTATTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCtt                                          <  2:76648/149‑1 (MQ=11)
                                                                                                       tggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACt                                                 <  2:335400/140‑1 (MQ=17)
                                                                                                        ggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTca                                        >  5:326454/1‑148 (MQ=32)
                                                                                                            tttCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCAc                                                                               <  3:149416/105‑1 (MQ=32)
                                                                                                            tttCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                                    >  8:73782/1‑148 (MQ=32)
                                                                                                            tttCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                                    >  2:155623/1‑148 (MQ=32)
                                                                                                            tttCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                                    >  4:136502/1‑148 (MQ=32)
                                                                                                            tttCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                                    >  7:192870/1‑148 (MQ=32)
                                                                                                              tccaGGTCGGCTCCATGCAGACTGTCGTTATCATTTCTAAGCCTCCAACGTTTCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGa                                  <  2:39670/145‑1 (MQ=11)
                                                                                                              tCAAGGTCGGCTCAATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTGAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCGCAGCGAg                                 >  3:59306/1‑149 (MQ=12)
                                                                                                                aaGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGt                                <  7:127590/148‑1 (MQ=32)
                                                                                                                aaGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCGCAATGTTCAGGGTCAAGCTATAGTAAAGCTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGtt                               >  2:209615/1‑149 (MQ=12)
                                                                                                                 aGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTc                              <  7:68669/149‑1 (MQ=32)
                                                                                                                    tCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCCCCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAAt                           >  7:104725/1‑149 (MQ=17)
                                                                                                                    tCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCTGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAAt                           <  6:328295/149‑1 (MQ=17)
                                                                                                                    tCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCttt                                                                      <  7:152938/106‑1 (MQ=32)
                                                                                                                    tCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCttt                                                                      >  8:152938/1‑106 (MQ=32)
                                                                                                                        cTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTATCCGTCTTGCCGCGGGTACACTGCATCTACACAGCGAGTTCAATTTCa                       >  6:144418/1‑149 (MQ=17)
                                                                                                                          ccATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTCCACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGGCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACt                     <  5:143407/149‑1 (MQ=17)
                                                                                                                           cATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCAc                      <  3:300336/147‑1 (MQ=32)
                                                                                                                           cATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACt                     >  6:213134/1‑148 (MQ=32)
                                                                                                                           cATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACt                     <  7:265347/148‑1 (MQ=32)
                                                                                                                            aTGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACt                     >  5:181325/1‑147 (MQ=32)
                                                                                                                                 gACTGGCGTGCACACTTCTTAGCCTCCCACCTTTCATACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGtctc              <  2:260853/149‑1 (MQ=255)
                                                                                                                                  aCTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCg             >  5:35566/1‑149 (MQ=32)
                                                                                                                                   cTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGTTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCgg            <  3:67511/149‑1 (MQ=17)
                                                                                                                                   cTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCgg            <  3:192063/149‑1 (MQ=32)
                                                                                                                                    tGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGc                                               <  1:221029/113‑1 (MQ=32)
                                                                                                                                     gtagtttacacTTCTAAGCGTCCAACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCggg           <  2:194708/141‑1 (MQ=255)
                                                                                                                                        gTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGga       <  7:42740/149‑1 (MQ=32)
                                                                                                                                            acatTTTTAAGCCTCCCACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGAc    <  8:192870/148‑1 (MQ=11)
                                                                                                                                             |                                                                                                                                                     
CTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGACAG  >  NC_000913/3424505‑3424794

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