Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,970 A→G 5.1% intergenic (+368/+246) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,9700AG 5.1% 257.5 / 5.1 99intergenic (+368/+246)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (41/52);  new base G (2/3);  total (44/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

ATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAA  >  NC_000913/3269824‑3270109
                                                                                                                                                  |                                                                                                                                           
aTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAg                                                                                                                                            <  7:273471/148‑1 (MQ=255)
aTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCCCTGTCTTTTGGGTATAGTGTCGTGGACAGACATTCATCATTCTGCCCCACCAAAAGAAAAAACCCGCCGAAg                                                                                                                                            >  5:117730/1‑148 (MQ=18)
 tAAAGCGCACTCAACAATCTGAAACCCTCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGc                                                                                                                                           <  1:39013/148‑1 (MQ=255)
 tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGc                                                                                                                                           <  5:312857/148‑1 (MQ=255)
  aaaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGTGAAAACCAGCCGAAggg                                                                                                                                          >  3:211822/1‑146 (MQ=255)
   aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                         <  5:6444/148‑1 (MQ=255)
     acgcACTCAACAAGCAGAAACCCGCCGGATCGGGTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTTGGTATAGTGTCGTTGACATTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                                                                                                                          <  6:206941/144‑1 (MQ=11)
      cgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGGAAAAACCCGGCGGATCggggg                                                                                                                                      >  7:715/1‑146 (MQ=255)
        cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtttt                                                                                                                                    <  4:389641/148‑1 (MQ=255)
         aCTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                                   <  6:324099/148‑1 (MQ=255)
           tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                                 >  3:342248/1‑148 (MQ=255)
           tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                                 >  6:442235/1‑148 (MQ=255)
            caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGttttt                                                                                                                                   >  2:288492/1‑145 (MQ=255)
            caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                                >  6:193097/1‑148 (MQ=255)
            caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                                >  1:444428/1‑148 (MQ=255)
            caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGTGTCGTGGACAGTCATTCATCATTAGGCACATCAAAAAGTAAAAACCCGACGAAGCGGGGATTTACg                                                                                                                                >  1:414795/1‑148 (MQ=18)
            caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCATTTGGGTATACTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                                >  1:309568/1‑148 (MQ=255)
             aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTa                                                                                                                                  <  1:288492/145‑1 (MQ=255)
              acaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTActtt                                                                                                                              <  2:206361/148‑4 (MQ=17)
              acaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                               >  5:179531/1‑147 (MQ=255)
              acaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGttatgcgga                                                                                                                               >  5:206941/1‑140 (MQ=255)
                aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                                                            <  5:442235/148‑1 (MQ=255)
                aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                                                            <  5:193097/148‑1 (MQ=255)
                      aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAACCCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                                                                                        <  3:135070/146‑1 (MQ=255)
                      aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                                                                                        <  3:141095/146‑1 (MQ=255)
                      aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                                                                                                        <  6:179531/146‑1 (MQ=255)
                       aaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTg                                                                                                                     <  6:329460/148‑1 (MQ=255)
                       aaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTg                                                                                                                     <  4:342248/148‑1 (MQ=255)
                       aaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTg                                                                                                                     <  2:444428/148‑1 (MQ=255)
                         cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                                                                                                   <  7:65490/148‑1 (MQ=255)
                          ccgccgTACCGGTTTATTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                                                     <  2:172096/129‑1 (MQ=11)
                             ccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGACGCGGGTTTTTACGTAAACCAGGTGAAAGTg                                                                                                               >  3:131068/1‑148 (MQ=14)
                             ccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTg                                                                                                               <  7:372508/148‑1 (MQ=255)
                              ctgAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                                                                                              <  5:314069/146‑1 (MQ=255)
                              cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGa                                                                                                              <  2:232352/148‑1 (MQ=255)
                               ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                             <  8:21297/148‑1 (MQ=255)
                                  aCGGTTTATTACGTACATCAGGTAAAACTTACCTATAATCCGCTATCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGa                                                                                                          <  7:95714/148‑1 (MQ=255)
                                  aCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                                                            >  1:156801/1‑130 (MQ=255)
                                  aCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                                                            <  2:156801/130‑1 (MQ=255)
                                    ggTTTATTAAGTAAATAATGTAAAACTGACCGATACGCCGCTTTCTTTTGGTTATAGTGTTGTCGATAGTCTTTGATCTTTCTTCCCCTCCAAAAGCAAAAACCCGCCGCATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATa                                                                                                        <  2:309568/148‑1 (MQ=11)
                                       ttattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGTTATAGTGTCGTGGAAAGTCATTAATCTTTATGCCCCTCCAAAATAAAAAACCACACGAATCGGGTTTTTACGTAAAACAGGTGAAACTGAcatcgaaa                                                                                                      >  6:283645/1‑140 (MQ=14)
                                             cGTACATCACTTAATAATGACCGATATTCCGCTTTCATTTGTGTATAGTGTCGTCGACTGTCTTTGTTCTTTCTGCCACTCCAAAATTATAAACCCGCCGTATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctt                                                                                                <  8:116710/147‑1 (MQ=2)
                                                 cATCAGGTAAAACTTCACGATAAGCCGCGTCCTTTTGGTTCTCTTCTCGCGGCCAGTCTATCCTCTTTCTCCCCCTCAAAAATAACAAACCCGCCGAACCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttcttt                                                                                           <  4:7886/148‑1 (MQ=2)
                                                         gaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGtttttgc                                                                                                                                 >  1:55713/2‑100 (MQ=11)
                                                         gaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGtttttgc                                                                                                                                 <  2:55713/101‑3 (MQ=11)
                                                         aaaaATGACCGATATTCAGGTTTCTTTTTAGGTTAGTGTCGGCGACTGCGTCTCATCTATATCGCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                                                                    <  5:412322/147‑1 (MQ=2)
                                                          aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGc                                                                                                                                           <  1:417146/91‑1 (MQ=33)
                                                          aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGc                                                                                                                                           >  2:417146/1‑91 (MQ=33)
                                                                agcGATATGCTGCTTTCTTTTTTTTTTAGTGTCGTGGAAAGTGTTTCTTATTTCTGCCACTCAAAAAGTACAAACCCGCCGTATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCg                                                                             <  8:328407/145‑1 (MQ=11)
                                                                  cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                                                                                                         <  7:141598/101‑1 (MQ=35)
                                                                  cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                                                                                                         >  8:141598/1‑101 (MQ=35)
                                                                         ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                                  >  5:268174/1‑85 (MQ=33)
                                                                             tttactttGGTTATATTGTCGTCGATATTCATTTATATTTATGGCACTCAAAAAGTTTAAACCCGCCGTATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcatt                                                                 <  4:211822/141‑1 (MQ=2)
                                                                                  tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                         >  6:203902/1‑101 (MQ=25)
                                                                                  tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                         <  5:203902/101‑1 (MQ=17)
                                                                                           aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgcgt                                                                               >  1:435459/1‑118 (MQ=17)
                                                                                             tgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTgg                                                                          >  8:127244/1‑121 (MQ=38)
                                                                                             tgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTgg                                                                          <  7:127244/121‑1 (MQ=38)
                                                                                               tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATTAAAAACCGGCCGAAGCGGGTTTTTAAGTAAAACATTTGAAACACACCGATAAGCCGCTTTCTTTTTGGTTTAGTGTCGTGGACAGTCATTCATCttt                                                          <  4:141344/135‑1 (MQ=11)
                                                                                                cGTGGACAGTCATTCATCTTTCTGCCCCTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTGATTCGTGTTTc                                                         >  3:141344/1‑135 (MQ=11)
                                                                                                  tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctt                                                                                            >  4:25352/1‑98 (MQ=37)
                                                                                                   ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                                                                                                   >  2:443829/1‑74 (MQ=39)
                                                                                                   ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                                                                                                   <  1:443829/74‑1 (MQ=39)
                                                                                                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctt                                                                                            <  2:217035/95‑1 (MQ=37)
                                                                                                     aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctt                                                                                            >  1:217035/1‑95 (MQ=37)
                                                                                                             tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaa                                                                                                                   >  7:274075/1‑64 (MQ=40)
                                                                                                             tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                          <  8:358083/121‑1 (MQ=38)
                                                                                                             tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                          >  7:358083/1‑121 (MQ=38)
                                                                                                              catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                          >  7:59044/1‑88 (MQ=38)
                                                                                                              catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                         >  4:144739/1‑89 (MQ=38)
                                                                                                              catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                         <  3:144739/89‑1 (MQ=38)
                                                                                                                            tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                          <  1:439977/106‑1 (MQ=40)
                                                                                                                            tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                          >  2:439977/1‑106 (MQ=40)
                                                                                                                            tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                  >  2:344885/1‑146 (MQ=37)
                                                                                                                             ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                           <  1:315095/72‑1 (MQ=39)
                                                                                                                             ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                           >  2:315095/1‑72 (MQ=255)
                                                                                                                             ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                   >  1:353327/1‑128 (MQ=38)
                                                                                                                             ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                   <  2:353327/128‑1 (MQ=38)
                                                                                                                             ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg               >  4:14556/1‑148 (MQ=37)
                                                                                                                             ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg               >  2:419625/1‑148 (MQ=37)
                                                                                                                             ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGcc                                                     >  1:231543/1‑110 (MQ=17)
                                                                                                                             ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGcc                                                     <  2:231543/110‑1 (MQ=17)
                                                                                                                             ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAATTAAAACAGCTGAAACTGACCGATAAGACGGTTTCTTTTGGTTATAGTGTCGTGGACAGGGATTAATAATTATTCACCTCCAAAAGCAAAAACCCGCCGAAGCGGGTGTTTGCg               >  4:217106/1‑148 (MQ=11)
                                                                                                                               aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                 <  7:20057/144‑1 (MQ=37)
                                                                                                                                caagctgaaaCCCGCCGTATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATACGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa            <  6:117730/141‑1 (MQ=11)
                                                                                                                                aaatcaaaaaCCCGCCGCACCGGGTTATTACGTAAAACAGGTCAAACTGACAGATAAGCCGCTTTCTTTTGGGTATAGTGTCGGGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa            <  5:194728/143‑1 (MQ=11)
                                                                                                                                aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                >  8:20057/1‑144 (MQ=37)
                                                                                                                                aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa            <  3:14556/148‑1 (MQ=37)
                                                                                                                                aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa             >  5:442423/1‑147 (MQ=37)
                                                                                                                                aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTc                                                                              >  3:75196/1‑82 (MQ=18)
                                                                                                                                 aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCg                                                                                                   >  6:215203/1‑60 (MQ=255)
                                                                                                                                   gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAcc         <  1:344885/148‑1 (MQ=38)
                                                                                                                                       aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaa                                              <  2:380696/107‑1 (MQ=35)
                                                                                                                                       aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaa                                              >  1:380696/1‑107 (MQ=35)
                                                                                                                                        atcccGCCGCACCGGGTTATTACGTACAACAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTg    <  4:131068/146‑1 (MQ=11)
                                                                                                                                         aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATc                                                             >  3:437171/1‑90 (MQ=38)
                                                                                                                                          cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc          <  5:144323/140‑1 (MQ=35)
                                                                                                                                          cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaa  <  6:442423/148‑1 (MQ=37)
                                                                                                                                          cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaa  <  1:419625/148‑1 (MQ=37)
                                                                                                                                                  |                                                                                                                                           
ATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAA  >  NC_000913/3269824‑3270109

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: