Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,956 T→A 12.4% 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,9560TA12.4% 173.4 / 5.4 74intergenic (+354/+260)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (27/35);  new base A (7/3);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAA  >  NC_000913/3269809‑3270100
                                                                                                                                                   |                                                                                                                                                
ccgcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGt                                                                                                                                                  >  8:397589/1‑148 (MQ=255)
ccgcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTa                                                                                                                                                 <  4:259855/149‑1 (MQ=255)
 cgcATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTa                                                                                                                                                 >  7:82464/1‑148 (MQ=255)
         gCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACccgc                                                                                                                                        >  8:31624/1‑149 (MQ=255)
           aTGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgcc                                                                                                                                       >  7:368942/1‑148 (MQ=255)
              aaTAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAg                                                                                                                                   >  6:324684/1‑149 (MQ=255)
                tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCg                                                                                                                                 <  5:324684/149‑1 (MQ=255)
                  aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                                                                                               <  1:428672/149‑1 (MQ=255)
                  aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                <  3:195480/148‑1 (MQ=255)
                  aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                <  7:33475/148‑1 (MQ=255)
                       cACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                          <  8:29085/149‑1 (MQ=255)
                        aCTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaaa                                                                                                                                              >  3:126445/1‑128 (MQ=255)
                        aCTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTaaaa                                                                                                                                              <  4:126445/128‑1 (MQ=255)
                         cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                        <  7:123131/149‑1 (MQ=255)
                         cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTATTCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                          >  8:213957/1‑147 (MQ=255)
                          tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                <  4:8578/140‑1 (MQ=255)
                          tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                >  3:8578/1‑140 (MQ=255)
                          tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                       >  1:69502/1‑149 (MQ=255)
                            aacaaTGAGAAACCCGCCGCACCGGGTTATTACGTACAACAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                       <  1:145874/147‑1 (MQ=11)
                            aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                      <  4:176619/148‑1 (MQ=255)
                            aacaaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCCTCGAAAGAGAGTCATGTTTGTAAACCGACAAATCTAAAAACAAGCAGACGCGGGtttgaaagt                                                                                                                      >  4:212896/1‑142 (MQ=12)
                               aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                  <  8:82464/149‑1 (MQ=255)
                                 tCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATCTTTCTGCCCCTCCAAAAGAAAAAAACCGCAGAAGCGGGTTTTTACGTaaaa                                                                                                                  >  8:336539/1‑147 (MQ=255)
                                   tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGttt                                                                                                                            <  1:173829/135‑1 (MQ=255)
                                   tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGttt                                                                                                                            >  2:173829/1‑135 (MQ=255)
                                   tGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTCTCGTCGACAGTCTTTCATCGTTGAGCCCCTCAAAAAGTAAAAACCCGCCGAATGCGGTTTTTAAGTAAAACAgg                                                                                                              >  5:401901/1‑149 (MQ=17)
                                     aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                          <  2:245777/135‑1 (MQ=255)
                                     aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                          >  1:245777/1‑135 (MQ=255)
                                     aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                             <  7:397589/148‑1 (MQ=255)
                                        cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                         <  5:359304/149‑1 (MQ=255)
                                         ccgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                        <  4:354192/149‑1 (MQ=255)
                                          cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTAATTCATCTTTCTGCCCCTACAAAAGAAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGGGatcgg                                                                                                       >  5:55493/1‑145 (MQ=14)
                                             cgGAACTGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTTGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                    <  1:421407/149‑1 (MQ=255)
                                                aaCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGa                                                                                                 <  7:213957/149‑1 (MQ=255)
                                                 aCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                <  2:69502/149‑1 (MQ=255)
                                                    gTTTATTAGGTTCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATTGTGTCGTGGACAGTAATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGTTTTTTa                                                                                                                         <  5:66499/121‑1 (MQ=14)
                                                     tttATTACGTACATCAGGTTAATCTTACCGATAAGCCGCTTTCTTTTGTGTTTTGTGTCGTCGACAGTCTTTTATCTTTCTGCCCCTCCAAAATAAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGc                                                                                            <  2:4091/149‑1 (MQ=11)
                                                               aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaa                                                                                                                                              <  1:43864/89‑1 (MQ=25)
                                                               aCATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAaaaa                                                                                                                                              >  2:43864/1‑89 (MQ=255)
                                                                cctCATGTAAAAATGACCGATAAGCCGCTTTCTTTTGCGTTTAGTCTCGTGGCCAGTGCTTCCTCTTTCTGCCCCTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCtttt                                                                                 <  3:107419/147‑1 (MQ=11)
                                                                       tAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                        >  3:60899/1‑103 (MQ=14)
                                                                          aaatGACAGATACGCCGCGTTCTTTTGGGTATAGTTTGGTTGACAGTAATTGATCTTTGAGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgt                                                                        <  3:212896/145‑1 (MQ=11)
                                                                          aaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGaaa                                                                                                         >  7:39198/1‑115 (MQ=17)
                                                                            cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtttta                                                                                                                          <  2:26014/96‑2 (MQ=18)
                                                                            cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGtttta                                                                                                                          >  1:26014/1‑95 (MQ=18)
                                                                             tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTTGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                          <  8:244144/95‑1 (MQ=18)
                                                                             tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                          <  2:78977/95‑1 (MQ=33)
                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttc                                                                                     <  7:429308/119‑1 (MQ=35)
                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttc                                                                                     >  8:429308/1‑119 (MQ=35)
                                                                                                              tCGTGGACAGTCTTTCTAATTTCTCGCCCTCCAAAAGAATAAACCCGCCGAACCGGGTTATTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                   <  6:145552/149‑1 (MQ=11)
                                                                                                                        tgattcatATTTCTGCCCCTCCAAAATTACAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAcccg                         <  7:31624/147‑2 (MQ=11)
                                                                                                                        tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtc                                                           >  8:350700/1‑115 (MQ=39)
                                                                                                                        tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtc                                                           <  7:350700/115‑1 (MQ=39)
                                                                                                                        tcattcatCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                >  2:429505/1‑78 (MQ=255)
                                                                                                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtg                                                                       >  1:296004/1‑102 (MQ=38)
                                                                                                                         cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtg                                                                       <  2:296004/102‑1 (MQ=38)
                                                                                                                          attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                 >  8:324545/1‑91 (MQ=37)
                                                                                                                          attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                       >  7:25560/1‑117 (MQ=39)
                                                                                                                          attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattc                                                       <  8:25560/117‑1 (MQ=39)
                                                                                                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                                     >  7:193766/1‑51 (MQ=35)
                                                                                                                               tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcttt                                                                                      >  5:93167/1‑81 (MQ=38)
                                                                                                                               tCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcttt                                                                                      <  6:93167/81‑1 (MQ=255)
                                                                                                                                 tttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                               >  1:329279/1‑86 (MQ=38)
                                                                                                                                 tttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                               <  2:329279/86‑1 (MQ=38)
                                                                                                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                 <  1:240157/79‑1 (MQ=39)
                                                                                                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                 >  2:240157/1‑79 (MQ=39)
                                                                                                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                   <  6:192097/93‑1 (MQ=40)
                                                                                                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                   >  5:192097/1‑93 (MQ=40)
                                                                                                                                       ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                               >  8:213467/1‑80 (MQ=39)
                                                                                                                                        cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                               <  8:299017/79‑1 (MQ=39)
                                                                                                                                        cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                               >  7:299017/1‑79 (MQ=39)
                                                                                                                                        cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgc                                                                        >  5:189692/1‑85 (MQ=37)
                                                                                                                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                   <  7:390954/122‑1 (MQ=38)
                                                                                                                                         ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                   >  8:390954/1‑122 (MQ=38)
                                                                                                                                               aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAACCCGCTTGCGTATGGGAATAGTGGCGTGGCCCCTCCATCCTCTATCTGCCCCTACAACAGCAAACACCCGACAAATCGGGGGTGTACGtcac  >  6:181742/1‑146 (MQ=2)
                                                                                                                                                   tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTg                                                                  >  2:14070/1‑81 (MQ=14)
                                                                                                                                                   tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                                              <  8:224905/85‑1 (MQ=14)
                                                                                                                                                   tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                                              >  7:224905/1‑85 (MQ=14)
                                                                                                                                                   |                                                                                                                                                
CCGCATCCGGCATGAATAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAA  >  NC_000913/3269809‑3270100

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: