Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,954 T→C 8.5% intergenic (+354/+260) 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,9540TC 8.5% 198.7 / 4.3 86intergenic (+354/+260)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (34/40);  new base C (3/5);  total (37/49)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAG  >  NC_000913/3269808‑3270102
                                                                                                                                                  |                                                                                                                                                    
cgcATCCGGGATGAATACAGCGCACTCCACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaa                                                                                                                                                    <  2:915420/149‑1 (MQ=255)
cgcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGCTCAGCCGCTTTCTTTTGGGTATAGTGTAGTGGACCGTCATTCCTCTTTCTGCCCCTCAAAAAGTaa                                                                                                                                                    >  1:613550/1‑149 (MQ=25)
  cATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaaa                                                                                                                                                  <  1:1263467/149‑1 (MQ=255)
   atcgcgCATGACTAAAGCGCACTCAACACTCTGAATCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTTCCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAGGTaaata                                                                                                                                                 <  2:635451/144‑3 (MQ=21)
      cGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAccc                                                                                                                                              <  2:736881/149‑1 (MQ=255)
      cGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAccc                                                                                                                                              <  1:926330/149‑1 (MQ=255)
      cGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAccc                                                                                                                                              >  1:349596/1‑149 (MQ=255)
              attaaGCGCACTCAACCATCCGGAGCCCGCCGGAGCGGTTTATTACTCACAGCGGCTAAAACTGACCGATAAGCCGCTTTCTTTGGGGCATGGTGTCGTGCACAGTCATTCATCTTCCTCCCCCTCCAAAAGTAAAAACCCGCCGAAGc                                                                                                                                      <  2:808951/146‑1 (MQ=17)
                aaaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTCCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                    <  2:245384/149‑1 (MQ=255)
                  aGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTTGAGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgggc                                                                                                                                  >  1:647594/1‑148 (MQ=255)
                    cgcACTCAACAATCTGGAACCCGCCGGAAGGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCAGCCGAAGCGGGttt                                                                                                                                >  1:849751/1‑149 (MQ=255)
                      cACTCAACAATCTGCAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                              <  1:577433/149‑1 (MQ=255)
                         tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATAAGGTAAAACTGACCGATACGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCCTATTTCTGCCCCTCCAAAAGTAAAAATCTGCCGAAGCGGGTTTTTagg                                                                                                                           >  2:1166422/1‑147 (MQ=18)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACGGTCATTCATCTTTCTGCCACTCCAAAAGTAAAAACCCCCCGAAGCGGGgttttcggga                                                                                                                         >  1:425111/1‑144 (MQ=18)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                                         >  2:191250/1‑149 (MQ=255)
                            acaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTActttt                                                                                                                        <  1:429934/149‑5 (MQ=18)
                              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                      <  2:855894/149‑1 (MQ=255)
                              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                      <  1:687906/149‑1 (MQ=255)
                              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                      <  1:1047625/149‑1 (MQ=255)
                                tCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAGAACCCGCTGAAGCGGGTTTTTACGTAAAAcc                                                                                                                    >  2:551693/1‑148 (MQ=255)
                                tCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTGCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACa                                                                                                                    >  2:551745/1‑149 (MQ=255)
                                 ctgagccccGCCGGAACTCTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGTGGACAGTCATCCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                                                                                   <  2:1142322/143‑1 (MQ=17)
                                 cTGAACCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAg                                                                                                                   <  1:522708/149‑1 (MQ=255)
                                   gAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                 <  1:684443/149‑1 (MQ=255)
                                     aaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAAGTAAAACAGGTGa                                                                                                               <  1:331301/149‑1 (MQ=255)
                                            gggAAGGGGTTGTTACGTACATCAGGTAAACCGGGCCGATAGGCCGCTTTCTTTTGGGTATAGTGCCGTGGACAGTCATTCATCTTTCTGCACCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAGGTAAAACAGGTGAACCTGag                                                                                                        <  2:990479/148‑2 (MQ=14)
                                                  ggTTTATTACGTACATCGGGTAAAACTGACCGATTAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAGCAGGTGAAACTGACCGATaa                                                                                                  <  2:59575/149‑1 (MQ=255)
                                                      tattaCGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAGTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGATACTGACCGATAAGACg                                                                                              >  2:910162/1‑149 (MQ=255)
                                                       attaCGTACATCAGGTAAACCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGc                                                                                             <  2:892630/149‑1 (MQ=12)
                                                          aCGCACATCAGGTAAAACCGAACCATCCGCCGCTTTCCTTGCGGTAGAGTGTCGTGGACAGGCCTGCGTCTTTCCGCCCCTCCCAAACTACATCCCCGCCCAAGCGGGTTGTTACGTAGACCAGGTGAAACTGACCGATAAGCCGcttt                                                                                          <  2:314393/149‑1 (MQ=2)
                                                             tACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttctt                                                                                       <  1:153412/149‑1 (MQ=255)
                                                              aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                    >  2:112769/1‑103 (MQ=25)
                                                               ccgcAGGAAAGACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                     <  2:1328481/146‑1 (MQ=255)
                                                               cagcAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                     <  1:1034552/146‑1 (MQ=38)
                                                               cATCAGGTAAATCTGACCGGTAAGCCTCTTTCTTTTGGGTGTATCTTCGTTGACAGTCACTCATCTTTCTGCCCTGTCAAAAGGAAAAGCCCGCCGCATCGGGTTTTTACGTAAAGCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                     <  2:390527/149‑1 (MQ=11)
                                                               cATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                     <  1:1160425/149‑1 (MQ=39)
                                                                aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                    <  2:1339580/149‑1 (MQ=39)
                                                                 ccAGGTAAGACGGACCGAAAAGTCGCTTTCGTTTGGGTATAGTGTCGGGGACAGTCATTAATCTTTCGTCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                   <  2:117005/148‑1 (MQ=11)
                                                                     gTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                                                               <  2:420391/149‑1 (MQ=37)
                                                                     gTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGtat                                                                               <  1:478392/149‑1 (MQ=38)
                                                                      tAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTgtgcgat                                                                              >  2:1318761/1‑143 (MQ=255)
                                                                      tAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                   >  1:609152/1‑144 (MQ=11)
                                                                       aacaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGGCGTGGACAGTCATTCATCTTTCTGCCCCTCCCAAAGTAAAACCCCGCCGAAGCGTGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                                             <  1:551745/146‑1 (MQ=11)
                                                                       aaaaCTGACCGATAAGCCGCTTTCTTTTGGGTATATTGTCGTGCACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAATCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                                             <  1:551693/149‑1 (MQ=255)
                                                                            tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCg                                                                        <  1:961296/149‑1 (MQ=38)
                                                                             gagcGATATGCCGCTTTCTGTTGGGGATAGTGTCGTGGACAGCCATGTAACTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                       <  2:610544/146‑1 (MQ=11)
                                                                              gccGATAGGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                      <  2:116739/148‑1 (MQ=17)
                                                                                  aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGTGTCGTGGACa                                                                  >  1:446947/1‑149 (MQ=38)
                                                                                   tAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGAATAGTGTCGTggcaag                                                                 >  1:1140447/1‑145 (MQ=12)
                                                                                     aGCCGCTTTCTTTTGGGTATTCTGTCGTGGACAGTCGTCCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtc                                                               <  2:349596/149‑1 (MQ=11)
                                                                                     aGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtc                                                               <  1:765038/149‑1 (MQ=38)
                                                                                       ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                                                                                                  >  1:1226783/1‑96 (MQ=37)
                                                                                        cGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCCACTCATAGATCGCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcatt                                                            <  1:1326110/149‑1 (MQ=11)
                                                                                         tctttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCAGCACCTGCCGGAATAAGACCCCGCCGAAGTGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                           <  1:1266279/148‑1 (MQ=11)
                                                                                           tttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattcat                                                         <  1:70221/149‑1 (MQ=38)
                                                                                                    ggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGATCGATAAGCCGCTTTCTTTTGGGTATAGCATCGTGAACCGTCATACATCtttatccc                                                >  2:1268216/1‑144 (MQ=11)
                                                                                                        taGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACGGACCGATAAGCCGCTGTCATTTGGGTATAGCGGCGTGGAAAGTCAGTCAGCTTACTGCCCCTc                                            >  2:679624/1‑149 (MQ=11)
                                                                                                           tgtCGTGGACAGTCACTCTACTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCag                                         >  1:993782/1‑148 (MQ=255)
                                                                                                             tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                    >  1:433360/1‑56 (MQ=38)
                                                                                                             tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTGCGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCAAGTCATTCATCTTTCTGCCCCTCCcaaa                                       >  1:184775/1‑149 (MQ=11)
                                                                                                             tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                       >  1:1194080/1‑149 (MQ=38)
                                                                                                             tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                       <  1:293816/149‑1 (MQ=38)
                                                                                                             tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                       >  1:433567/1‑149 (MQ=38)
                                                                                                             acGTGGACAGTCATTTATCTTTCGGCCCCTCCAAAAGTAAAAACCCGTCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                       <  2:90815/148‑1 (MQ=255)
                                                                                                              cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                      >  1:428999/1‑149 (MQ=38)
                                                                                                                   aCAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                             >  2:524774/1‑57 (MQ=38)
                                                                                                                     aGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAc                               >  1:623754/1‑149 (MQ=38)
                                                                                                                       tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                                         >  2:499063/1‑57 (MQ=38)
                                                                                                                        ggtatgccTTTCGGCCCCTCCAAAAGGAAAAACCCGCCGGAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                            <  2:1115612/142‑1 (MQ=11)
                                                                                                                        cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                            <  2:623754/149‑1 (MQ=17)
                                                                                                                         attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc                           <  1:389782/149‑1 (MQ=38)
                                                                                                                             atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                       >  2:121137/1‑133 (MQ=17)
                                                                                                                                tttCTGCCCCCCCAAAAGCAGACACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                    <  2:901486/149‑1 (MQ=255)
                                                                                                                                 ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                   <  1:515773/149‑1 (MQ=38)
                                                                                                                                   cagCCCATCCAACACGAAAAACCCGCCGAAGCGGGTTTTCACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGt                 <  2:451819/147‑1 (MQ=11)
                                                                                                                                        ccTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa            <  1:1112225/149‑1 (MQ=38)
                                                                                                                                         cTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCg                                                                                                      >  2:952774/1‑58 (MQ=255)
                                                                                                                                         cTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCCAGTCACTCAGCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGGTTTTAc           >  2:95852/1‑149 (MQ=11)
                                                                                                                                          tCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTACAACAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg          <  2:691159/149‑1 (MQ=255)
                                                                                                                                           cccaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGGATAGTGTCGTGGACGGTCAGTCATTTTTTTGCTCCTCCACAGGGCATAGCCGGCCGCAGCGGGTGTTTACGt         >  2:1082639/4‑149 (MQ=11)
                                                                                                                                            cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGGTTTTACGTa        >  2:1299477/1‑149 (MQ=255)
                                                                                                                                            cAAAAGTAAAAACCCGCCGAAGCGGGGTTTTACGTAAAACAGGGGAAACTGACCGGTAAGCCGGGTTCTATCGGGGATAGTGACGGGGACAGTCATGCATCCTCCCGCTCACCCAATAGCAAAAACCTGCTGAAGTGGGTTCTTACggg        >  2:1299562/1‑147 (MQ=255)
                                                                                                                                             aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGATTTGACGTaa       >  1:205645/1‑149 (MQ=255)
                                                                                                                                             aaaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa       >  1:894706/1‑149 (MQ=38)
                                                                                                                                              aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGACAGTCATCTTTCTGCCCATCCAAAAGCAAAAACCCGCCGAAGCGGGCTTTTACGGaaa      >  2:1260817/1‑149 (MQ=11)
                                                                                                                                                aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                                                               >  2:253283/1‑74 (MQ=40)
                                                                                                                                                  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg  <  1:717468/149‑1 (MQ=38)
                                                                                                                                                  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGTGGATAGTGTCGTGGACAGTCCTTCATCTTTCTGCCCCCCCAAAAGCAAAACCCGGCCGCAGCGGGGTTTTACGCAAATCaa  >  2:715079/1‑148 (MQ=11)
                                                                                                                                                  |                                                                                                                                                    
CGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAG  >  NC_000913/3269808‑3270102

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