Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,969 A→C 7.3% 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 7.3% 185.1 / 2.4 74intergenic (+367/+247)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base C (3/3);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAA  >  NC_000913/3269823‑3270110
                                                                                                                                                   |                                                                                                                                             
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                                                                                                                                  >  3:314840/1‑149 (MQ=255)
           cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                                  <  5:216490/149‑1 (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:42813/149‑1 (MQ=255)
                 aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                            <  4:314840/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                                                                                                                   <  7:249311/149‑1 (MQ=255)
                          cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                                   <  6:233436/149‑1 (MQ=255)
                           ccgccgTAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTGC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                                  <  2:216791/149‑1 (MQ=255)
                               cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGACGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                              <  6:303719/149‑1 (MQ=255)
                                ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                          >  7:287829/1‑120 (MQ=255)
                                ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                          <  8:287829/120‑1 (MQ=255)
                                  aaCGGTTTATTACGTACATCAGGT‑AAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                          <  8:14541/149‑1 (MQ=255)
                                    cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                                                                                                         <  7:146761/149‑1 (MQ=255)
                                                 aCATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                                                         <  8:169407/104‑1 (MQ=18)
                                                 aCATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                                                         >  7:169407/1‑104 (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)
                                                    tctgttAAAACTGACCGATACGCCGCGTTC‑TTTTTGGTCTAGTGGCGGTGAGACTCATTCAGCTTTCTGCACCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                          <  7:412181/143‑1 (MQ=11)
                                                                 aCCCATCAGCCGCTCTC‑AGTTGCCTACACTCTCGCGGACACCCACGCATCTCTCAGCCCCTCCAAAAGCACAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                            <  8:65000/149‑1 (MQ=2)
                                                                          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                                                     >  3:175703/1‑115 (MQ=17)
                                                                                                                           cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGccc                                                     <  4:175703/115‑1 (MQ=17)
                                                                                                                           cccTCAACAAGTTTAAACCCGCAGTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGTTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCTTCGGCGGGttttt                   <  6:297877/149‑1 (MQ=11)
                                                                                                                            ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                  >  7:2503/1‑149 (MQ=37)
                                                                                                                            ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGa                                                                          >  7:232306/1‑93 (MQ=18)
                                                                                                                              tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                                 >  3:300145/1‑100 (MQ=40)
                                                                                                                              tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                                 <  4:300145/100‑1 (MQ=40)
                                                                                                                                   aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc           <  8:2503/149‑1 (MQ=38)
                                                                                                                                     gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                             <  3:287216/81‑1 (MQ=17)
                                                                                                                                     gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                             >  4:287216/1‑81 (MQ=17)
                                                                                                                                       aaaaaCCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                       >  1:302540/1‑101 (MQ=11)
                                                                                                                                       aaaaaCCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                       <  2:302540/101‑1 (MQ=11)
                                                                                                                                        aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt      <  8:37347/149‑1 (MQ=37)
                                                                                                                                            cccTCCGCACCGGGTTTTTACGTAAATAATTTCAAACTGACCGATAAGACGCTTTCTTTTGGGTATAGTGTCGTGGACATTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa  <  5:349099/149‑1 (MQ=11)
                                                                                                                                            cccGCCGCATCGGTTTTTTACGTAAATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa  <  5:370899/149‑1 (MQ=11)
                                                                                                                                             ccgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGt                                                                                       >  4:231429/1‑63 (MQ=35)
                                                                                                                                                 cgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGa                                                                                                           >  1:150275/1‑39 (MQ=40)
                                                                                                                                                   |                                                                                                                                             
AATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTC‑TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAA  >  NC_000913/3269823‑3270110

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