Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,969 A→C 5.4% intergenic (+367/+247) 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,9690AC 5.4% 961.6 / 28.6 336intergenic (+367/+247)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (92/217);  new base C (2/20);  total (94/242)
Fisher's exact test for biased strand distribution p-value = 4.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

ATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCG  >  NC_000913/3269820‑3270117
                                                                                                                                                     |                                                                                                                                                    
 tGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGATCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGaa                                                                                                                                                     <  2:1128054/150‑1 (MQ=255)
 actatatAGCCCACTCCACAATCAGAAACCCGCCGGATCGGTTTATTACGTAAACCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                                                                                                                                     <  2:3671421/147‑1 (MQ=14)
    aTAAAGCGCACTCAACAATCTTAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCg                                                                                                                                                  <  1:5072349/150‑1 (MQ=255)
    aTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCg                                                                                                                                                  <  1:3582446/150‑1 (MQ=255)
     tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                                 <  2:5268149/150‑1 (MQ=255)
      aaaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACTTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                                                                                                                <  1:3963616/150‑1 (MQ=255)
      aaaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                                                                                                                <  2:4398523/150‑1 (MQ=255)
      aaaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                                                                                                                <  2:2711516/150‑1 (MQ=255)
      aaaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTGTTGGGTCTAGTGTCGTGGACCGTCCTTCATCGTTCTGCCCCGCCAAAAGCAAAAACCCGCCGCAGCggg                                                                                                                                                >  2:5098683/1‑150 (MQ=18)
         gcgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGTTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                                                                                                                             <  1:5605840/150‑1 (MQ=255)
          cgcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtttt                                                                                                                                            >  1:2822816/1‑150 (MQ=255)
            cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                                          <  1:5337920/150‑1 (MQ=255)
            cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCGGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTa                                                                                                                                          <  1:1992236/150‑1 (MQ=18)
             aCTCAACAATATTAAACCCGCCGTATCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                                         <  1:3279940/150‑1 (MQ=255)
                      tCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                                                                                                                <  1:3511211/150‑1 (MQ=255)
                       cTGACACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                                                                                                               <  2:575194/150‑1 (MQ=255)
                       cTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                                                                                                               >  1:71839/1‑150 (MQ=255)
                        tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGGAAAAACCCGCCGAAGCTGGATTGTAAGAAAAACAGGt                                                                                                                              >  1:3182443/1‑150 (MQ=17)
                        tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                              >  1:5538624/1‑150 (MQ=255)
                         gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTg                                                                                                                             >  2:190183/1‑150 (MQ=255)
                          aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGa                                                                                                                            <  1:3437470/150‑1 (MQ=255)
                            aCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                                          >  2:2509804/1‑150 (MQ=255)
                               cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCACCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTg                                                                                                                       >  2:616425/1‑150 (MQ=255)
                                 ccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGAc                                                                                                                     >  2:1267562/1‑150 (MQ=255)
                                  cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                                    <  1:4008850/150‑1 (MQ=255)
                                  cgGAACGGTATATTAGGTACATCTGGTAAATCTGAGCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                                    <  2:3144427/150‑1 (MQ=17)
                                    gatcGGTTTTTTACGTACATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGa                                                                                                                  <  2:1397335/147‑1 (MQ=11)
                                       cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                               <  2:1078303/150‑1 (MQ=255)
                                       cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                               <  2:3912654/150‑1 (MQ=255)
                                       cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                               <  1:4642667/150‑1 (MQ=255)
                                       cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                               <  2:4249209/150‑1 (MQ=255)
                                         gTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGc                                                                                                             <  1:3160625/150‑1 (MQ=255)
                                             attaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGct                                                                                                         >  1:427182/1‑150 (MQ=255)
                                             attaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGct                                                                                                         <  1:2794328/150‑1 (MQ=255)
                                                 cGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttct                                                                                                     <  2:5598077/150‑1 (MQ=255)
                                                    aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                                  <  2:4583842/150‑1 (MQ=40)
                                                    aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                                  <  2:4881684/150‑1 (MQ=40)
                                                    aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                                  <  2:5167667/150‑1 (MQ=40)
                                                    aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAAAGGTTAAAATTAACGCTAAGCCGggttctttg                                                                                                  >  1:1017237/1‑149 (MQ=11)
                                                     cATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                                 <  2:4974827/150‑1 (MQ=39)
                                                      aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                <  1:2058336/150‑1 (MQ=38)
                                                      aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                <  1:1556005/150‑1 (MQ=38)
                                                      aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                <  1:2558837/150‑1 (MQ=38)
                                                          ggTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                                                                            <  1:5789807/150‑1 (MQ=37)
                                                            tgaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAAATTTATTCCCCACCAAAAGAAAAAAACCGCCGAAGCGGGTTTTTAAGTAAAACAGGTGAAAAAGACCGATAATCCGCGTTCTTTTGGGTATAg                                                                                          >  1:1516971/3‑150 (MQ=11)
                                                              aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtg                                                                                        <  1:2434132/150‑1 (MQ=35)
                                                                 cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCg                                                                                     <  2:3972598/150‑1 (MQ=37)
                                                                  tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                                    <  2:3404069/150‑1 (MQ=37)
                                                                  tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                                    <  2:3198639/150‑1 (MQ=37)
                                                                  tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTAATTAATCTTTGTTACCCTACAAAATTAAAAACCCGCGGAAGCGGGTTTTTACGTAATACAGGTGAAAATTTACGATATGCGGGGTTGTTTTGGGTATAggttggt                                                                                    >  2:5204378/1‑144 (MQ=2)
                                                                      cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTAGAGTGTCGTGGAc                                                                                >  2:4777075/1‑150 (MQ=255)
                                                                      cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAAAGTCATTCATCTATCTGCCCCAACAAAATAAAAAACCAGCCGAATCGCCTTTTAACGTAAAACAGGTTAAACTTACCGAAAAGCCGCTTTCTTTTGTGTAAACTGTCTTgtcc                                                                                >  2:1612011/1‑147 (MQ=11)
                                                                       gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                                               <  2:5335525/150‑1 (MQ=37)
                                                                       gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                                               <  1:1022375/150‑1 (MQ=37)
                                                                             ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcatt                                                                         <  2:223478/150‑1 (MQ=37)
                                                                                  ttctttTGGGTATAGTGTCGTGGACAGTCATTAAACTTTCTGCCCCACCAAAAGTAAAAAACCGCCGAAGCGGGTTTTTAAGTAAAACAGGTGAAAATGACCGATAAGCCGCTTTCTTTTTGGTATAGTGTCGTGGACAGtcaatcatat                                                                    >  1:1830711/1‑148 (MQ=11)
                                                                                     ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTc                                                                 <  2:1954023/150‑1 (MQ=37)
                                                                                      tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCt                                                                <  1:3711256/150‑1 (MQ=37)
                                                                                      tttGCGTATAGTGTCGTGGACAGACATTCAACTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTATTACGTAAAACAGCTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCt                                                                <  1:5648780/150‑1 (MQ=11)
                                                                                       ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCATTCATCTTTCTg                                                               <  1:476168/150‑1 (MQ=37)
                                                                                        tGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                              <  1:1723821/150‑1 (MQ=37)
                                                                                        tGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                              <  1:4479612/150‑1 (MQ=37)
                                                                                               aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCa                                                       <  1:1628115/150‑1 (MQ=37)
                                                                                                gtgtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaa                                                      <  1:3859604/150‑1 (MQ=37)
                                                                                                  gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                    <  1:885053/150‑1 (MQ=37)
                                                                                                     gTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                                                               >  2:3771211/1‑120 (MQ=17)
                                                                                                      tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCACCAAAAGCaa                                                >  2:2689366/1‑150 (MQ=255)
                                                                                                       ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaa                                               <  1:3538585/150‑1 (MQ=37)
                                                                                                            gTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAccc                                          <  1:3898569/150‑1 (MQ=37)
                                                                                                             tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGTGTCGTGTAAAGTAATTCATATTTATGCCCCTCCAAAAGCaaaatcacg                                         >  2:151735/1‑145 (MQ=11)
                                                                                                                ttcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgccg                                      >  1:695014/1‑150 (MQ=37)
                                                                                                                 tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                                     >  2:4363693/1‑150 (MQ=37)
                                                                                                                  catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                    <  2:412574/150‑1 (MQ=37)
                                                                                                                  catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                    <  2:5041809/150‑1 (MQ=37)
                                                                                                                     cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                 <  1:4934868/150‑1 (MQ=37)
                                                                                                                      ttacTGCCACTCCAAAACTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                <  2:798248/147‑1 (MQ=255)
                                                                                                                         cTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                             <  2:2097443/150‑1 (MQ=37)
                                                                                                                         cTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                             <  2:2834684/150‑1 (MQ=37)
                                                                                                                          tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTGTAGAGTCTAGAGGAAAGTCAATCAAAATTAAGCACCAGCAATAGAAAAAACCAGCCGAGGCTGGttt                            >  2:754256/1‑150 (MQ=2)
                                                                                                                            ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGttttt                          <  1:4161254/150‑1 (MQ=255)
                                                                                                                             cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                         <  1:5258110/150‑1 (MQ=37)
                                                                                                                              ccTCCAAAAGTAAAACCCCGCCGAAGCGGTTTTTTACGTAAACCCGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttttgc                        <  2:2914269/150‑3 (MQ=255)
                                                                                                                              ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                        <  2:4196493/150‑1 (MQ=37)
                                                                                                                                 ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa                     <  1:2583356/150‑1 (MQ=37)
                                                                                                                                  cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                    >  1:1353217/1‑150 (MQ=37)
                                                                                                                                    caaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc                  <  1:3599643/149‑1 (MQ=38)
                                                                                                                                        gacaaaCCCGCCGCATCGGTTTTTTACGTACATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAgg              <  2:3527170/147‑1 (MQ=11)
                                                                                                                                         aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTAAAGTGTCGAGTCCAGTCATTCATATTTCTGACCCTACAAAAGACAATAACCGCCGAAGCGGGTGTTTACGTATACAAggg             >  2:2267057/1‑149 (MQ=11)
                                                                                                                                           acaCCCGCCGCACCGGTATACTACGTAAAACAGGTAAAACTGACCGAATCGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTACTGCCCCTACCCAAGCAAAAACCCGCAGAAGCTGGTTTTTTCGTAAACAAGGTGa           <  2:5438319/148‑1 (MQ=11)
                                                                                                                                           aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGa           >  1:3138046/1‑150 (MQ=37)
                                                                                                                                           aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTCAATCATCGTTCAGCCCCTACAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGGTAACCAGGTGa           >  1:4630919/1‑150 (MQ=11)
                                                                                                                                             aCCCGCCGATGAGGGTTTTTTCGTACAACAGGTGTAAATGACGGATAAGACGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTGAACTTTCCGCCCCTGAAATAGCTAAATGCCGCCGAAGGGGGTTTTTACGTAAACCTGGTGaaa         <  2:1622700/150‑1 (MQ=14)
                                                                                                                                             aCCCGCCGACGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa         <  2:3881021/150‑1 (MQ=37)
                                                                                                                                              cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc        <  2:3110388/150‑1 (MQ=37)
                                                                                                                                                  ccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAc    <  1:616425/150‑1 (MQ=37)
                                                                                                                                                  ccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAc    <  2:1978224/150‑1 (MQ=37)
                                                                                                                                                  ccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAc    >  2:2649485/1‑150 (MQ=37)
                                                                                                                                                   cgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc   <  2:2338608/150‑1 (MQ=37)
                                                                                                                                                   cgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGAcc   >  2:4062258/1‑150 (MQ=17)
                                                                                                                                                    gAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAAGCAGGTGAAACTGAcag  >  2:4828762/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCG  >  NC_000913/3269820‑3270117

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