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

CGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAC  >  NC_000913/3269810‑3270101
                                                                                                                                                   |                                                                                                                                                 
cgcATCCGTCATGAATAAAGCGCACTCAACAATCTGCAACCAGCCGGAACGGTTTATTACGTAAATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaa                                                                                                                                                 <  5:444/149‑1 (MQ=255)
 gcATCCGGCATGAATAAAGCGGACTCAACATTCTGAAACCCGAAGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGATTTCTTCTGCGTATATTGTCGTTGACAGTCATTCAGCTTTCTTCCCCTCCATAAGTata                                                                                                                                                <  6:50658/149‑3 (MQ=21)
    tCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAc                                                                                                                                             <  1:252151/149‑1 (MQ=255)
            gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaa                                                                                                                                                <  8:318836/138‑1 (MQ=255)
            gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaa                                                                                                                                                >  7:318836/1‑138 (MQ=255)
            gAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgccg                                                                                                                                       >  8:249311/1‑147 (MQ=255)
             aaTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAg                                                                                                                                    >  6:333133/1‑149 (MQ=255)
              aTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGAACATCAGGTAAAAATGACCGATAAGACGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCCAAGc                                                                                                                                   >  1:61898/1‑149 (MQ=18)
               tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATTAAAAACCCTCCGAAGCg                                                                                                                                  >  3:256730/1‑149 (MQ=255)
                 aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                 <  5:333133/148‑1 (MQ=255)
                    cgcACTCAACAATCTGAAACCCGCCGGAACGGATTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGGGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttt                                                                                                                             <  8:85731/149‑1 (MQ=255)
                       aCTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGATATAGTATCCTGGACGGTAAGTCAAATTTATTCCTCGCCAAAAGTAAAAACCCGGCGCACCGGGttatt                                                                                                                           >  8:109519/1‑148 (MQ=18)
                        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTGGTTGCTCTTACTGCGCCACCAAAAGTAAAAACCCGACGCAGGGGGGTTTTAc                                                                                                                         >  7:65000/1‑149 (MQ=18)
                        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                         <  5:216490/149‑1 (MQ=255)
                        cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                         >  3:314840/1‑149 (MQ=255)
                         tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTAAATCAGGTAAAACTGAACGATAAGCCGATTTCTTTTGGGTAAAGTGTCTTGGAAAGTAATTAATCTTTCTTCCCCGACAAAATAAAAAAACCACCGAAGCGGGTTTTTACg                                                                                                                        >  6:349099/1‑149 (MQ=255)
                         tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTAAATCAGGTAAAACTGAACGATAAGCCGATATCTTTTGTATATATAGTAGTTGAAAGTAATTAATCGTTGTGACCCTAAAAAATAAAAAACCAACCGAAGCGGGTTTTTACg                                                                                                                        >  6:402936/1‑149 (MQ=255)
                          caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                        >  1:42292/1‑148 (MQ=255)
                          caacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                       >  7:14541/1‑148 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                       <  3:194810/148‑1 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                       >  3:119039/1‑148 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                       >  3:113681/1‑148 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                                      <  4:113681/149‑1 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                                      <  4:119039/149‑1 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTACAAAATAAAAAACCCGCCCAAGCGGGTTTTTACGt                                                                                                                       >  6:370899/1‑148 (MQ=255)
                           aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTAGTGTACAGAAATTCAGATTTCTCCACCAACAAAAGTAAAAACCCGCAGAAGAGTGTTTTTAGGt                                                                                                                       >  1:82501/1‑148 (MQ=17)
                              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTTATTAATCTTTCTGGACGGCAAAAAGTAAAAACCCGAGGAAGCGGGTTTTTCCGTtaaa                                                                                                                   >  7:166430/1‑149 (MQ=17)
                              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                   <  4:314840/149‑1 (MQ=255)
                              aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                   <  4:42813/149‑1 (MQ=255)
                               aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGCACAGTAATTCAACTATCTGCCCCACCAAAAGTAAAAAACCGCCCAAGCGGGTTTTTAAGTaaaa                                                                                                                   >  1:311633/1‑148 (MQ=17)
                               aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAACCGCTTTCTTTTGGGTATAGTGAAGAGGACACTAATTAAAACAACTCCCCCACCAAAAGAAAAAAACCCCCGAAGAGGGTTTTTACaaaaaa                                                                                                                   >  3:308837/1‑148 (MQ=12)
                                gatGAAACCGTCCGGAAAGGTTTATTACGTAAATCAGGTAAAACTGAACGATAAGCCTCTTTCTTTTGTTTATTGTGTCGTGGAGAGTCTTTCTTCTTTCTTCCCCTCTAAAAGTAAAAACCCGCCTAAGCGGGTTTTTACGTAAAAc                                                                                                                  <  4:256730/146‑1 (MQ=255)
                                    aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                              <  2:42292/148‑1 (MQ=255)
                                    aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTATTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                              <  6:98981/148‑1 (MQ=255)
                                       cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                          <  6:233436/149‑1 (MQ=255)
                                       cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                          <  7:249311/149‑1 (MQ=255)
                                        ccgccgTAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTGCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                         <  2:216791/149‑1 (MQ=255)
                                            cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGACGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                     <  6:303719/149‑1 (MQ=255)
                                             ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                 <  8:287829/120‑1 (MQ=255)
                                             ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                 >  7:287829/1‑120 (MQ=255)
                                               aaCGGTTTATTACGTACATCAGGTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                 <  8:14541/149‑1 (MQ=255)
                                                 cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                                                                                                <  7:146761/149‑1 (MQ=255)
                                                              aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                                                >  7:169407/1‑104 (MQ=18)
                                                              aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                                                <  8:169407/104‑1 (MQ=18)
                                                                ctcacgaaaaaCTGACCGATAAGCCTCTTTCATTTTCTGTTCGTCTCGTGGAAAGTCTTTCATCTTTCTGCCCCTCCAAAAGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                 <  8:344593/142‑1 (MQ=11)
                                                                aTCAGGTACAACTGACCGATAAGCCGCTTTCTTTTGTGTCTAGTGTCGTGGACAGTCATGCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                 <  7:109519/149‑1 (MQ=11)
                                                                aTCAGGTAAACCTGACCGATACGCCGCTTTCTATTGGGTCTAGTCTCGCGGACAGTCATTCATCTTTCTGCCACTCCAAAAGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                 <  8:391463/149‑1 (MQ=11)
                                                                aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAc                                                                                                                                             <  5:299330/89‑1 (MQ=40)
                                                                aTCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAc                                                                                                                                             >  6:299330/1‑89 (MQ=40)
                                                                 tctgttAAAACTGACCGATACGCCGCGTTCTTTTTGGTCTAGTGGCGGTGAGACTCATTCAGCTTTCTGCACCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                 <  7:412181/143‑1 (MQ=11)
                                                                              aCCCATCAGCCGCTCTCAGTTGCCTACACTCTCGCGGACACCCACGCATCTCTCAGCCCCTCCAAAAGCACAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                   <  8:65000/149‑1 (MQ=2)
                                                                                cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaaaa                                                                                                                                              >  6:413592/1‑72 (MQ=34)
                                                                                cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaaaa                                                                                                                                              <  5:413592/72‑1 (MQ=34)
                                                                                       ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                <  7:333357/127‑1 (MQ=35)
                                                                                        cGCTTTCATTTTTTTATTGTGTCGTCGACAGTCTTTATTCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcat                                                          <  3:217957/148‑1 (MQ=11)
                                                                                             tctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTTCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                    <  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)
                                                                                                                                                   |                                                                                                                                                 
CGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAC  >  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: