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

AATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAA  >  NC_000913/3269823‑3270110
                                                                                                                                                   |                                                                                                                                             
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                                                                                                                                  >  3:314840/1‑149 (MQ=255)
           cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                                  <  5:216490/149‑1 (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:42813/149‑1 (MQ=255)
                 aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                            <  4:314840/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                                                                                                                   <  7:249311/149‑1 (MQ=255)
                          cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                                   <  6:233436/149‑1 (MQ=255)
                           ccgccgTAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTGCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                                  <  2:216791/149‑1 (MQ=255)
                               cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCACTCCAAAAGTAAAAACCCGCCGACGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                              <  6:303719/149‑1 (MQ=255)
                                ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                          >  7:287829/1‑120 (MQ=255)
                                ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                          <  8:287829/120‑1 (MQ=255)
                                  aaCGGTTTATTACGTACATCAGGTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                          <  8:14541/149‑1 (MQ=255)
                                    cGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                                                                                                         <  7:146761/149‑1 (MQ=255)
                                                 aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                                                         <  8:169407/104‑1 (MQ=18)
                                                 aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                                                         >  7:169407/1‑104 (MQ=18)
                                                   ctcacgaaaaaCTGACCGATAAGCCTCTTTCATTTTCTGTTCGTCTCGTGGAAAGTCTTTCATCTTTCTGCCCCTCCAAAAGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                          <  8:344593/142‑1 (MQ=11)
                                                   aTCAGGTACAACTGACCGATAAGCCGCTTTCTTTTGTGTCTAGTGTCGTGGACAGTCATGCATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                          <  7:109519/149‑1 (MQ=11)
                                                   aTCAGGTAAACCTGACCGATACGCCGCTTTCTATTGGGTCTAGTCTCGCGGACAGTCATTCATCTTTCTGCCACTCCAAAAGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                          <  8:391463/149‑1 (MQ=11)
                                                    tctgttAAAACTGACCGATACGCCGCGTTCTTTTTGGTCTAGTGGCGGTGAGACTCATTCAGCTTTCTGCACCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                          <  7:412181/143‑1 (MQ=11)
                                                                 aCCCATCAGCCGCTCTCAGTTGCCTACACTCTCGCGGACACCCACGCATCTCTCAGCCCCTCCAAAAGCACAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTg                                                                            <  8:65000/149‑1 (MQ=2)
                                                                          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                                                     >  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)
                                                                                                                                                   |                                                                                                                                             
AATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAA  >  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: