Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,967 C→T 5.0% intergenic (+365/+249) 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,9670CT 5.0% 281.9 / 4.3 99intergenic (+365/+249)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (42/52);  new base T (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 = 9.34e-01

Only 100 of 104 total aligned reads displayed.

GAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAA  >  NC_000913/3269822‑3270109
                                                                                                                                                 |                                                                                                                                              
gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgccg                                                                                                                                               >  7:21297/1‑147 (MQ=255)
gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgccg                                                                                                                                               >  1:47953/1‑147 (MQ=255)
  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)
    aaaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgccg                                                                                                                                               <  7:294340/143‑1 (MQ=255)
    aaaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgccg                                                                                                                                               >  8:294340/1‑143 (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)
              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)
                                               cGTACATCACTTAATAATGACCGATATTCCGCTTTCATTTGTGTATAGTGTCGTCGACTGTCTTTGTTCTTTCTGCCACTCCAAAATTATAAACCCGCCGTATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctt                                                                                                <  8:116710/147‑1 (MQ=2)
                                                   cATCAGGTAAAACTTCACGATAAGCCGCGTCCTTTTGGTTCTCTTCTCGCGGCCAGTCTATCCTCTTTCTCCCCCTCAAAAATAACAAACCCGCCGAACCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttcttt                                                                                           <  4:7886/148‑1 (MQ=2)
                                                           gaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGtttttgc                                                                                                                                 <  2:55713/101‑3 (MQ=11)
                                                           gaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGtttttgc                                                                                                                                 >  1:55713/2‑100 (MQ=11)
                                                           aaaaATGACCGATATTCAGGTTTCTTTTTAGGTTAGTGTCGGCGACTGCGTCTCATCTATATCGCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                                                                    <  5:412322/147‑1 (MQ=2)
                                                            aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGc                                                                                                                                           >  2:417146/1‑91 (MQ=33)
                                                            aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGc                                                                                                                                           <  1:417146/91‑1 (MQ=33)
                                                                  agcGATATGCTGCTTTCTTTTTTTTTTAGTGTCGTGGAAAGTGTTTCTTATTTCTGCCACTCAAAAAGTACAAACCCGCCGTATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCg                                                                             <  8:328407/145‑1 (MQ=11)
                                                                    cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                                                                                                         >  8:141598/1‑101 (MQ=35)
                                                                    cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                                                                                                         <  7:141598/101‑1 (MQ=35)
                                                                           ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                                  >  5:268174/1‑85 (MQ=33)
                                                                                    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                                                          >  7:358083/1‑121 (MQ=38)
                                                                                                               tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCttt                                                          <  8:358083/121‑1 (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)
                                                                                                                                                 |                                                                                                                                              
GAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAA  >  NC_000913/3269822‑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: