Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,328,183 T→A 9.2% intergenic (+100/+56) yhbY → / ← greA RNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,328,1830TA 9.2% 467.2 / 27.9 139intergenic (+100/+56)yhbY/greARNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor
Reads supporting (aligned to +/- strand):  ref base T (64/62);  new base A (9/4);  total (73/66)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 140 total aligned reads displayed.

GGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGAT  >  NC_000913/3328035‑3328325
                                                                                                                                                    |                                                                                                                                               
 gTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATa                                                                                                                                                >  3:622842/1‑149 (MQ=255)
     tttATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATCAAAGGCCGCATAGCgg                                                                                                                                           <  4:738128/149‑1 (MQ=255)
     tttATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCgg                                                                                                                                           >  1:728277/1‑149 (MQ=255)
     tttATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCgg                                                                                                                                           >  1:401098/1‑149 (MQ=255)
       tATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGcc                                                                                                                                         >  3:403179/1‑149 (MQ=255)
        aTCGCCCAACTAAAGAATGTAAAATCTGGCTGCCACACTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGTGGGGTTTTCCGCAGGCAAGAGAGCAAAATGCCACGCTCTGCTCGTTGATAAAAGGCCGCATAGCGGCCt                                                                                                                                        >  4:1023816/1‑149 (MQ=255)
          cGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCttt                                                                                                                                      <  3:327300/149‑1 (MQ=255)
             ccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTc                                                                                                                                   <  3:846078/149‑1 (MQ=255)
             ccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTc                                                                                                                                   >  2:810083/1‑149 (MQ=255)
             ccAACTAAAGAACGGAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTc                                                                                                                                   >  2:671722/1‑149 (MQ=255)
               aaCTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCct                                                                                                                                 >  3:830189/1‑149 (MQ=255)
                 cTAAAGAACGTAAAATCTCGCTGCCACACTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCcttt                                                                                                                               <  1:87818/149‑1 (MQ=255)
                  tAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTATGCGGCCtttt                                                                                                                                     >  4:406368/1‑142 (MQ=255)
                   aaaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttct                                                                                                                             <  2:330760/149‑1 (MQ=255)
                   aaaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttct                                                                                                                             <  1:963694/149‑1 (MQ=255)
                    aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttctt                                                                                                                            >  2:180057/1‑149 (MQ=255)
                    aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttctt                                                                                                                            >  3:574186/1‑149 (MQ=255)
                    aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCcttt                                                                                                                               <  4:381466/146‑1 (MQ=255)
                    aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCcttt                                                                                                                               >  3:381466/1‑146 (MQ=255)
                    aaGAACGTAAAATCGCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttctt                                                                                                                            <  2:434221/149‑1 (MQ=255)
                      gAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttc                                                                                                                              >  1:845653/1‑145 (MQ=255)
                      gAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttc                                                                                                                              <  2:845653/145‑1 (MQ=255)
                      gAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTa                                                                                                                          <  1:948700/149‑1 (MQ=255)
                      gAACGTAAAATCTCGCTGCCACGCTAAGATTACCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTa                                                                                                                          <  1:788668/149‑1 (MQ=255)
                        aCGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACa                                                                                                                        <  4:67195/149‑1 (MQ=255)
                           tAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATa                                                                                                                     <  3:815399/149‑1 (MQ=255)
                           tAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATa                                                                                                                     >  2:969541/1‑149 (MQ=255)
                           tAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATa                                                                                                                     >  1:726248/1‑149 (MQ=255)
                           tAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATa                                                                                                                     <  2:401098/149‑1 (MQ=255)
                              aaTCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACAt                                                                                                                  >  1:325189/1‑149 (MQ=255)
                                         cACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGa                                                                                                       <  1:853179/149‑1 (MQ=255)
                                          aCGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAg                                                                                                      <  1:432334/149‑1 (MQ=255)
                                              tAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTAtt                                                                                                  >  3:48266/1‑149 (MQ=255)
                                                 gATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTggg                                                                                               <  3:356716/149‑1 (MQ=255)
                                                   ttATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTa                                                                                             >  2:585024/1‑149 (MQ=255)
                                                    tATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTaa                                                                                            <  1:475242/149‑1 (MQ=255)
                                                      tCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAAtt                                                                                          <  2:504841/149‑1 (MQ=255)
                                                        cTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCt                                                                                        >  4:329836/1‑149 (MQ=255)
                                                        cTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCt                                                                                        <  3:344279/149‑1 (MQ=255)
                                                          aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTa                                                                                      <  3:86965/149‑1 (MQ=255)
                                                             gTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAg                                                                                   <  2:450394/149‑1 (MQ=255)
                                                                  acacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAtt                                                                              >  1:1079770/1‑149 (MQ=255)
                                                                    acaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATc                                                                                                           <  3:315156/118‑1 (MQ=255)
                                                                    acaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTc                                                                             >  2:478154/1‑148 (MQ=255)
                                                                      aTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAc                                                                          >  3:137761/1‑149 (MQ=255)
                                                                      aTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAc                                                                          >  4:182290/1‑149 (MQ=255)
                                                                       tttGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAcc                                                                         >  3:304871/1‑149 (MQ=255)
                                                                       tttGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAcc                                                                         >  3:1080470/1‑149 (MQ=255)
                                                                       tttGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAcc                                                                         >  2:660061/1‑149 (MQ=255)
                                                                            tgtgTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAAt                                                                    <  4:448243/149‑1 (MQ=255)
                                                                            tgtgTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAAt                                                                    >  1:530313/1‑149 (MQ=255)
                                                                              tgtAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTa                                                                  <  2:883742/149‑1 (MQ=255)
                                                                              tgtAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTa                                                                  >  2:189386/1‑149 (MQ=255)
                                                                              tgcaaaaCGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAAGTa                                                                  >  3:718551/4‑149 (MQ=255)
                                                                                tAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACt                                                                <  4:137761/149‑1 (MQ=255)
                                                                                 aaaaCGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACtt                                                               >  2:7633/1‑149 (MQ=255)
                                                                                     cGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCaaa                                                           <  1:466067/149‑1 (MQ=255)
                                                                                      gAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAAt                                                          <  1:973956/149‑1 (MQ=255)
                                                                                      gAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAAt                                                          >  1:156672/1‑149 (MQ=255)
                                                                                        ggggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTc                                                        >  1:1045773/1‑149 (MQ=255)
                                                                                          ggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTa                                                      >  1:280775/1‑149 (MQ=255)
                                                                                            ttttCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACt                                                    >  1:949677/1‑149 (MQ=255)
                                                                                                 cGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTcgcc                                               <  3:188325/149‑1 (MQ=255)
                                                                                                          gagaGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGt                                                                                              >  1:979312/1‑93 (MQ=255)
                                                                                                            gagCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTCCGCCGCCCGGCGttt                                    <  4:796384/149‑1 (MQ=255)
                                                                                                                aaaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgat                                <  3:293022/149‑1 (MQ=255)
                                                                                                                aaaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgat                                <  4:518826/149‑1 (MQ=255)
                                                                                                                  aaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgatga                              <  3:103827/149‑1 (MQ=255)
                                                                                                                     gCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTgg                                                                                                <  2:8661/80‑1 (MQ=255)
                                                                                                                      ccACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACcac                          >  2:252083/1‑149 (MQ=255)
                                                                                                                      ccACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACcac                          <  4:875455/149‑1 (MQ=255)
                                                                                                                       cACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGt                                                                                 >  4:806084/1‑93 (MQ=255)
                                                                                                                       cACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGt                                                                                 <  3:806084/93‑1 (MQ=255)
                                                                                                                       cACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCaaa                                                           <  2:1080350/115‑1 (MQ=255)
                                                                                                                       cACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCaaa                                                           >  1:1080350/1‑115 (MQ=255)
                                                                                                                          gCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCacaaca                      <  1:189386/149‑1 (MQ=255)
                                                                                                                            tctGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcatc                    >  2:557829/1‑149 (MQ=255)
                                                                                                                             ctGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcatca                   <  3:234128/149‑1 (MQ=255)
                                                                                                                              tGTTCTTTGATAAAATGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcatcat                  >  2:976405/1‑149 (MQ=255)
                                                                                                                               gTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACatcatc                 <  1:660061/149‑1 (MQ=255)
                                                                                                                                  cGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttc              <  1:969541/149‑1 (MQ=255)
                                                                                                                                  cGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttc              >  2:626318/1‑149 (MQ=255)
                                                                                                                                    ttGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttctt            >  3:1006280/1‑149 (MQ=255)
                                                                                                                                     aggaaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttcttc           >  4:809473/5‑149 (MQ=255)
                                                                                                                                      ggaaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttcttct          >  4:326829/4‑149 (MQ=255)
                                                                                                                                      gATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttcttct          >  4:888843/1‑149 (MQ=255)
                                                                                                                                       gaaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttcttc           >  3:7530/3‑147 (MQ=255)
                                                                                                                                       gaaaaaGGCCGCTAAGCGGCCTTATTACTTTCTTTACAATACTTCAACATGTTGAGTATTGGGTAATTCTTACAGGTATTCGACCTTAATTACTTAATATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATcttcttctt         >  2:247434/3‑149 (MQ=255)
                                                                                                                                        tAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt        >  2:315104/1‑149 (MQ=255)
                                                                                                                                        aaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGtttt                                   >  2:959179/2‑122 (MQ=255)
                                                                                                                                         aaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTAc                                                                 <  2:464432/91‑1 (MQ=255)
                                                                                                                                         aaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTAc                                                                 >  1:464432/1‑91 (MQ=255)
                                                                                                                                         aaaaGGCCGCATAGCGGCCTTTTCCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCACCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTg       >  4:105204/1‑149 (MQ=255)
                                                                                                                                          aaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGc      >  3:979530/1‑149 (MQ=255)
                                                                                                                                          aaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGc      >  3:599607/1‑149 (MQ=255)
                                                                                                                                          aaaGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGc      <  1:180057/149‑1 (MQ=255)
                                                                                                                                          aaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCaaa                                                           <  4:321139/96‑1 (MQ=255)
                                                                                                                                          aaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCaaa                                                           >  3:321139/1‑96 (MQ=255)
                                                                                                                                           aaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGcc     >  3:514918/1‑149 (MQ=255)
                                                                                                                                              gCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGAt  <  4:48266/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGAT  >  NC_000913/3328035‑3328325

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: