Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,956 T→C 8.6% 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,9560TC 8.6% 192.1 / 6.0 77intergenic (+354/+260)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (29/37);  new base C (2/4);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAC  >  NC_000913/3269810‑3270101
                                                                                                                                                   |                                                                                                                                                 
cgcATCCGTCATGAATAAAGCGCACTCAACAATCTGCAACCAGCCGGAACGGTTTATTACGTAAATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaa                                                                                                                                                 <  5:444/149‑1 (MQ=255)
 gcATCCGGCATGAATAAAGCGGACTCAACATTCTGAAACCCGAAGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGATTTC‑TTCTGCGTATATTGTCGTTGACAGTCATTCAGCTTTCTTCCCCTCCATAAGTata                                                                                                                                                <  6:50658/149‑3 (MQ=21)
    tCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAc                                                                                                                                             <  1:252151/149‑1 (MQ=255)
            gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaa                                                                                                                                                <  8:318836/138‑1 (MQ=255)
            gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaa                                                                                                                                                >  7:318836/1‑138 (MQ=255)
            gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgccg                                                                                                                                       >  8:249311/1‑147 (MQ=255)
             aaTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAg                                                                                                                                    >  6:333133/1‑149 (MQ=255)
              aTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGAACATCAGGTAAAAATGACCGATAAGACGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCCAAGc                                                                                                                                   >  1:61898/1‑149 (MQ=18)
               tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATTAAAAACCCTCCGAAGCg                                                                                                                                  >  3:256730/1‑149 (MQ=255)
                 aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                 <  5:333133/148‑1 (MQ=255)
                    cgcACTCAACAATCTGAAACCCGCCGGAACGGATTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGGGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                                                                                                             <  8:85731/149‑1 (MQ=255)
                       aCTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGATATAGTATCCTGGACGGTAAGTCAAATTTATTCCTCGCCAAAAGTAAAAACCCGGCGCACCGGGttatt                                                                                                                           >  8:109519/1‑148 (MQ=18)
                        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTGGTTGCTCTTACTGCGCCACCAAAAGTAAAAACCCGACGCAGGGGGGTTTTAc                                                                                                                         >  7:65000/1‑149 (MQ=18)
                        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                         <  5:216490/149‑1 (MQ=255)
                        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                         >  3:314840/1‑149 (MQ=255)
                         tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTAAATCAGGTAAAACTGAACGATAAGCCGATTTC‑TTTTGGGTAAAGTGTCTTGGAAAGTAATTAATCTTTCTTCCCCGACAAAATAAAAAAACCACCGAAGCGGGTTTTTACg                                                                                                                        >  6:349099/1‑149 (MQ=255)
                         tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTAAATCAGGTAAAACTGAACGATAAGCCGATATC‑TTTTGTATATATAGTAGTTGAAAGTAATTAATCGTTGTGACCCTAAAAAATAAAAAACCAACCGAAGCGGGTTTTTACg                                                                                                                        >  6:402936/1‑149 (MQ=255)
                          caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                        >  1:42292/1‑148 (MQ=255)
                          caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGT‑AAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                       >  7:14541/1‑148 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                       <  3:194810/148‑1 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                       >  3:119039/1‑148 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                       >  3:113681/1‑148 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                                      <  4:113681/149‑1 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                                      <  4:119039/149‑1 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTACAAAATAAAAAACCCGCCCAAGCGGGTTTTTACGt                                                                                                                       >  6:370899/1‑148 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTAGTGTACAGAAATTCAGATTTCTCCACCAACAAAAGTAAAAACCCGCAGAAGAGTGTTTTTAGGt                                                                                                                       >  1:82501/1‑148 (MQ=17)
                              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTTATTAATCTTTCTGGACGGCAAAAAGTAAAAACCCGAGGAAGCGGGTTTTTCCGTtaaa                                                                                                                   >  7:166430/1‑149 (MQ=17)
                              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                   <  4:314840/149‑1 (MQ=255)
                              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                   <  4:42813/149‑1 (MQ=255)
                               aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGCACAGTAATTCAACTATCTGCCCCACCAAAAGTAAAAAACCGCCCAAGCGGGTTTTTAAGTaaaa                                                                                                                   >  1:311633/1‑148 (MQ=17)
                               aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAACCGCTTTC‑TTTTGGGTATAGTGAAGAGGACACTAATTAAAACAACTCCCCCACCAAAAGAAAAAAACCCCCGAAGAGGGTTTTTACaaaaaa                                                                                                                   >  3:308837/1‑148 (MQ=12)
                                gatGAAACCGTCCGGAAAGGTTTATTACGTAAATCAGGTAAAACTGAACGATAAGCCTCTTTC‑TTTTGTTTATTGTGTCGTGGAGAGTCTTTCTTCTTTCTTCCCCTCTAAAAGTAAAAACCCGCCTAAGCGGGTTTTTACGTAAAAc                                                                                                                  <  4:256730/146‑1 (MQ=255)
                                    aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                              <  2:42292/148‑1 (MQ=255)
                                    aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TATTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                              <  6:98981/148‑1 (MQ=255)
                                       cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                          <  6:233436/149‑1 (MQ=255)
                                       cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                          <  7:249311/149‑1 (MQ=255)
                                        ccgccgTAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTGC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                         <  2:216791/149‑1 (MQ=255)
                                            cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGACGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                     <  6:303719/149‑1 (MQ=255)
                                             ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                 <  8:287829/120‑1 (MQ=255)
                                             ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                 >  7:287829/1‑120 (MQ=255)
                                               aaCGGTTTATTACGTACATCAGGT‑AAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                 <  8:14541/149‑1 (MQ=255)
                                                 cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                                                                                                <  7:146761/149‑1 (MQ=255)
                                                              aCATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                                                >  7:169407/1‑104 (MQ=18)
                                                              aCATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                                                <  8:169407/104‑1 (MQ=18)
                                                                ctcacgaaaaaCTGACCGATAAGCCTCTTTCATTTTCTGT‑TCGTCTCGTGGAAAGTCTTTCATCTTTCTGCCCCTCCAAAAGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                 <  8:344593/142‑1 (MQ=11)
                                                                aTCAGGTACAACTGACCGATAAGCCGCTTTC‑TTTTGTGTCTAGTGTCGTGGACAGTCATGCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                 <  7:109519/149‑1 (MQ=11)
                                                                aTCAGGTAAACCTGACCGATACGCCGCTTTC‑TATTGGGTCTAGTCTCGCGGACAGTCATTCATCTTTCTGCCACTCCAAAAGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                 <  8:391463/149‑1 (MQ=11)
                                                                aTCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAc                                                                                                                                             <  5:299330/89‑1 (MQ=40)
                                                                aTCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAc                                                                                                                                             >  6:299330/1‑89 (MQ=40)
                                                                 tctgttAAAACTGACCGATACGCCGCGTTC‑TTTTTGGTCTAGTGGCGGTGAGACTCATTCAGCTTTCTGCACCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                 <  7:412181/143‑1 (MQ=11)
                                                                              aCCCATCAGCCGCTCTC‑AGTTGCCTACACTCTCGCGGACACCCACGCATCTCTCAGCCCCTCCAAAAGCACAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                   <  8:65000/149‑1 (MQ=2)
                                                                                cGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaaaa                                                                                                                                              >  6:413592/1‑72 (MQ=34)
                                                                                cGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaaaa                                                                                                                                              <  5:413592/72‑1 (MQ=34)
                                                                                       ccGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                <  7:333357/127‑1 (MQ=35)
                                                                                        cGCTTTC‑ATTTTTTTATTGTGTCGTCGACAGTCTTTATTCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcat                                                          <  3:217957/148‑1 (MQ=11)
                                                                                             tc‑tttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTTCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                    <  5:239399/133‑1 (MQ=37)
                                                                                                ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                     >  7:364654/1‑98 (MQ=37)
                                                                                                ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                     <  8:364654/98‑1 (MQ=37)
                                                                                                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                           <  4:25891/66‑1 (MQ=35)
                                                                                                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                           >  3:25891/1‑66 (MQ=35)
                                                                                                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                                          >  8:114074/1‑115 (MQ=35)
                                                                                                          aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                                          <  7:114074/115‑1 (MQ=35)
                                                                                                                      aGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                         >  5:134394/1‑72 (MQ=39)
                                                                                                                                     tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACt                                                                                                        >  4:204107/1‑58 (MQ=255)
                                                                                                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                 >  6:310901/1‑64 (MQ=255)
                                                                                                                                        cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                   >  5:223620/1‑76 (MQ=39)
                                                                                                                                        cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt           >  7:37347/1‑148 (MQ=37)
                                                                                                                                        cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGccc                                            <  4:175703/115‑1 (MQ=17)
                                                                                                                                        cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGccc                                            >  3:175703/1‑115 (MQ=17)
                                                                                                                                        cccTCAACAAGTTTAAACCCGCAGTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGTTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCTTCGGCGGGttttt          <  6:297877/149‑1 (MQ=11)
                                                                                                                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa         >  7:2503/1‑149 (MQ=37)
                                                                                                                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGa                                                                 >  7:232306/1‑93 (MQ=18)
                                                                                                                                           tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                        <  4:300145/100‑1 (MQ=40)
                                                                                                                                           tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                        >  3:300145/1‑100 (MQ=40)
                                                                                                                                                aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc  <  8:2503/149‑1 (MQ=38)
                                                                                                                                                  gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                    >  4:287216/1‑81 (MQ=17)
                                                                                                                                                  gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                    <  3:287216/81‑1 (MQ=17)
                                                                                                                                                   |                                                                                                                                                 
CGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAC  >  NC_000913/3269810‑3270101

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: ‑