Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,194 G→T 6.8% intergenic (+592/+22) 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,270,1940GT 6.8% 192.6 / 3.1 67intergenic (+592/+22)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base G (38/24);  new base T (3/2);  total (41/26)
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

CAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTG  >  NC_000913/3270063‑3270340
                                                                                                                                   |                                                                                                                                                  
cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGTGTAAAGCGTAGATGAAAGTCAATCATCATAAGACACCACAAGAAACAAAAACCGGAAGAAGAGGGTTTTTaaga                                                                                                                                    >  4:339213/1‑145 (MQ=2)
         aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGTTTTTgcttttggagggg                                                                                                                           >  1:159296‑M1/1‑135 (MQ=255)
              cgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGTTTTTgcttttggaggggcagaaa                                                                                                                     <  2:159296‑M1/149‑20 (MQ=255)
                     gCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAAGGCAATCATAGTTGGGCCCCGACAAAAACAATAACCCGCCGCAGGGGGTGGTTTCGGAAATGAGGTGGCACTGGCCGag                                                                                                              >  8:213113/1‑148 (MQ=2)
                        ggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATAGTTCTGCACCTACAAAAGCAAAAACCCGCAGCAGAGGGTTAATAAGTAAATCAGGGGAAAATGACCGATaat                                                                                                           >  4:189876/1‑148 (MQ=11)
                                     aacAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                                                                                                                               <  8:117631/98‑1 (MQ=14)
                                     aacAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATc                                                                                                                               <  1:128681/114‑1 (MQ=17)
                                       cAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCGTTCGGCCCCACCAAAAGCAAAAACCCGCCGCAGCCGGGTTTTACGTAAATCAGGTGAAAATGTCCGATAAGACGGGTTCTGtgggg                                                                                            >  1:370849/1‑145 (MQ=14)
                                              aaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                                                                 <  1:208055/105‑1 (MQ=32)
                                                   aCCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATATTTCTGCCCCACCAAAAGCAAAAACCCGCCGAAGACGGTTTTTACGTAAATCAGGTGAAAATGACCGATAAGACGGGTTCTGTCGTGGAAAGtcattcc                                                                                >  8:253769/1‑148 (MQ=17)
                                                       aTAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATc                                                                              >  7:210188/1‑147 (MQ=255)
                                                            ccGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGATACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                                                                         >  6:89246/1‑147 (MQ=255)
                                                             ctctttctttTGGGTATAGCTTCGGGGACAGTCACTCATGACTCTCCCGCTCCCAAAGCAACTACCCTGCTAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAg                                                                      <  5:110790/147‑1 (MQ=18)
                                                               ctttctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATAAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGc                                                                     >  7:18058/1‑148 (MQ=255)
                                                                 ttctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                                                                                                   <  2:169935/100‑1 (MQ=37)
                                                                 ttctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                                                                                                   >  1:169935/1‑100 (MQ=37)
                                                                   ctttTTGGTATAGGCTCGTGGACTGTCTTTGATCTTTGTGGCCCTCCAAACGGATATACCGGCCGAAGGGGGTTTTTACGTAATTCAGGTGAAACTGACCGATAAGGCGGGTTCTGTCGTGGACAGTGATTCATCTAGGCCAGCAATc                                                                 <  7:141502/148‑1 (MQ=18)
                                                                      ttGGGTATTGCGTCTTGTTCAGTCATGCATCTTTCTGCCAATCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGc                                                               <  3:363987/147‑1 (MQ=21)
                                                                        gggTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTACAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCAc                                                           >  8:99119/1‑149 (MQ=255)
                                                                         ggcatcgCGTCGGCGACATGCATTCATCTATATGGCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACt                                                          <  3:339213/143‑1 (MQ=21)
                                                                          gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTg                                                         >  7:181599/1‑149 (MQ=255)
                                                                              aGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCGCACTGGCt                                                      >  3:31442/1‑148 (MQ=255)
                                                                               gCGTCGTCGACAGTCATGCCTCCCTCTCCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCa                                                    <  5:380728/149‑1 (MQ=25)
                                                                                         cAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGaaa                                                                                                                       >  4:63177/1‑72 (MQ=34)
                                                                                         cAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGaaa                                                                                                                       <  3:63177/72‑1 (MQ=34)
                                                                                         cAGACATTCATCTTTCAGCCCCCCCCACAGCAAACCCCCGCCGAAGCGCGCTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTAc                                          <  3:1603/149‑1 (MQ=25)
                                                                                           gTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGCCCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTAccc                                        >  5:248813/1‑149 (MQ=255)
                                                                                              attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTgg                                                        >  4:325904/1‑130 (MQ=255)
                                                                                              attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTgg                                                        <  3:325904/130‑1 (MQ=255)
                                                                                              attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAagc                                                 >  5:30063/1‑137 (MQ=255)
                                                                                              attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAagc                                                 <  6:30063/137‑1 (MQ=255)
                                                                                              attcatCTTTCAGCCCCTCCAAACGCAAAAACCCGCGGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAACCTACCCgg                                      <  3:336834/148‑1 (MQ=255)
                                                                                               ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCCGCAATCGCTCACTGGCTCAAGAAGCCTACCCGGGt                                    >  2:88938/1‑149 (MQ=255)
                                                                                               ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCa                                                            >  7:293093/1‑125 (MQ=255)
                                                                                               ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCa                                                            <  8:293093/125‑1 (MQ=255)
                                                                                               ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                    >  8:287415/1‑149 (MQ=255)
                                                                                               ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                    >  1:125675/1‑149 (MQ=255)
                                                                                               ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                    >  2:94539/1‑149 (MQ=255)
                                                                                                  atCTTTCTGGCCCTTCAATAGGAATAACGCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACTGATAAGCGGTGTTCTGTGGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCa                                 <  2:315141/149‑1 (MQ=25)
                                                                                                  atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAATTCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                   >  5:7192/1‑147 (MQ=255)
                                                                                                     tttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                                                                                                       >  8:153194/1‑76 (MQ=35)
                                                                                                     tttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                                                                                                       <  7:153194/76‑1 (MQ=35)
                                                                                                     tttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTa                              >  4:114719/1‑149 (MQ=255)
                                                                                                       gcTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTa                              <  5:277322/146‑1 (MQ=255)
                                                                                                         tGCCCCTCCAAAAGCAAAAACCCGACGAACCGGGTTATTAAGTAAATCAGGTGAAACTGACCGATAACCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGAAATCGCTCACAGCCTCAACCAGCCTACCCGGGTTAAGTACtgg                          >  2:186694/1‑149 (MQ=25)
                                                                                                          gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTgg                                                        <  2:218480/118‑1 (MQ=255)
                                                                                                          gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTgg                                                        >  1:218480/1‑118 (MQ=255)
                                                                                                          gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGc                         >  2:12446/1‑149 (MQ=255)
                                                                                                           ccccTGCAACAGCATACACCCGCGGCGGCGGGTTTTTAGTTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGc                         <  3:146649/148‑1 (MQ=38)
                                                                                                            cccTCCAAAATAAGAAACCCGCCGTATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTtttgtgtatagtgtcgtggacagtcattcatctttctgcccctccaaaagcaaaaacccgcttcggcgggttttt                       >  7:306983‑M2/1‑74 (MQ=255)
                                                                                                            cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                                                                                                           <  7:275001/49‑1 (MQ=38)
                                                                                                               tCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTAc                    >  1:116351/1‑149 (MQ=255)
                                                                                                                ccAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTAcc                   >  4:131775/1‑149 (MQ=255)
                                                                                                                ccAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACAGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTAc                                          >  6:97453/1‑126 (MQ=255)
                                                                                                                ccAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACAGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTAc                                          <  5:97453/126‑1 (MQ=255)
                                                                                                                 caaacccAAAATCCCGCCGAAGCTGTGTTTTACGTATATCAGGGGACACCGCCCGATAAGCGGGTTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTAc                    <  1:408635/141‑1 (MQ=255)
                                                                                                                 cAAAAGCAAAAACCCGCCGAAGCGGGCTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCt                  >  2:6805/1‑149 (MQ=255)
                                                                                                                   aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCGCTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTa                >  3:50872/1‑149 (MQ=255)
                                                                                                                   aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTa                <  2:116351/149‑1 (MQ=255)
                                                                                                                    acgtaaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTAt               >  6:72283/5‑149 (MQ=38)
                                                                                                                    aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTAt               <  7:213113/149‑1 (MQ=255)
                                                                                                                    aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTAt               >  3:103192/1‑149 (MQ=255)
                                                                                                                      gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATg              >  7:327308/1‑148 (MQ=255)
                                                                                                                       taaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGGTCAGTACGGGCCGTACCTTATGaa            >  4:285723/2‑149 (MQ=38)
                                                                                                                        aaaaaCCCGCCGAAGCGGGTTTTTgcttttggaggggcagaaagatgaatgactgtccacgacactatacacaaaagaaagcggcttatcggtcagtttaacctgttttacgtaaaaacccgatacggcgggtttcttattttggaggg           >  7:306983‑M1/1‑24 (MQ=255)
                                                                                                                        aaaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAc           >  1:150206/1‑149 (MQ=255)
                                                                                                                             cccGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTa      >  3:62228/1‑149 (MQ=255)
                                                                                                                              ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAt     >  4:202468/1‑149 (MQ=255)
                                                                                                                                 ccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCGCTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTg  <  4:50872/149‑1 (MQ=255)
                                                                                                                                 ccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTg  <  8:327308/149‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                  
CAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTG  >  NC_000913/3270063‑3270340

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: