Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,196 A→G 7.1% intergenic (+594/+20) 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,1960AG 7.1% 191.2 / 5.3 70intergenic (+594/+20)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (39/26);  new base G (3/2);  total (42/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

CAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCT  >  NC_000913/3270063‑3270344
                                                                                                                                     |                                                                                                                                                    
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                                                            >  1:218480/1‑118 (MQ=255)
                                                                                                          gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTgg                                                            <  2:218480/118‑1 (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                   >  3:103192/1‑149 (MQ=255)
                                                                                                                    aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTAt                   <  7:213113/149‑1 (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)
                                                                                                                                     aGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCt  <  7:253769/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                    
CAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCT  >  NC_000913/3270063‑3270344

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: