Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 28,094 G→A V268M (GTG→ATG)  rihC → ribonucleoside hydrolase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091328,0940GA98.6% 519.6 / ‑4.6 140V268M (GTG→ATG) rihCribonucleoside hydrolase 3
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (75/62);  total (78/62)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 142 total aligned reads displayed.

GCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATC  >  NC_000913/27948‑28241
                                                                                                                                                  |                                                                                                                                                   
  cACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTg                                                                                                                                               <  2:1581030/151‑1 (MQ=255)
    cTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGat                                                                                                                                             >  2:163072/1‑151 (MQ=255)
     tACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGat                                                                                                                                             >  2:533429/1‑150 (MQ=255)
      aCCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGatat                                                                                                                                           >  4:390217/1‑151 (MQ=255)
      aCCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGat                                                                                                                                             <  1:455412/149‑1 (MQ=255)
        cGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCg                                                                                                                                         >  4:124277/1‑151 (MQ=255)
         gTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGa                                                                                                                                        <  3:517099/151‑1 (MQ=255)
         gTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGa                                                                                                                                        <  1:812013/151‑1 (MQ=255)
            gCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACgg                                                                                                                                     <  3:349824/151‑1 (MQ=255)
                cAGTATGCAAAGCGGCTTGCGAATGCACGATTTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCCGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGCTGGGTGCTATCGACGGtt                                                                                                                                   >  4:471012/1‑149 (MQ=255)
                 aGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGcc                                                                                                                                >  1:1479195/1‑151 (MQ=255)
                   tATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTg                                                                                                                              >  4:150611/1‑151 (MQ=255)
                    aTGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTgg                                                                                                                             >  2:1651572/1‑151 (MQ=255)
                      gCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTggg                                                                                                                            <  1:1014752/150‑1 (MQ=255)
                      gCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTggg                                                                                                                            >  1:1591727/1‑150 (MQ=255)
                      gCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGc                                                                                                                           <  1:643426/151‑1 (MQ=255)
                       cAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGc                                                                                                                           >  4:509464/1‑150 (MQ=255)
                              gCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAagccag                                                                                                                    >  1:828111/1‑150 (MQ=255)
                                ttGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAagccagcc                                                                                                                  >  2:971337/1‑150 (MQ=255)
                                ttGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAagccagcc                                                                                                                  >  1:569010/1‑150 (MQ=255)
                                    gAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGt                                                                                                             >  3:519996/1‑151 (MQ=255)
                                     aaTGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGt                                                                                                             >  2:1304104/1‑150 (MQ=255)
                                     aaTGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGt                                                                                                             >  1:451208/1‑150 (MQ=255)
                                     aaTGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTa                                                                                                            <  1:330966/151‑1 (MQ=255)
                                      aTGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCCGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTa                                                                                                            >  1:1506239/1‑150 (MQ=255)
                                              tctcTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGc                                                                                                    >  2:360554/1‑150 (MQ=255)
                                                tctGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCAtt                                                                                                 >  2:880068/1‑151 (MQ=255)
                                                 ctGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTg                                                                                                >  1:351398/1‑151 (MQ=255)
                                                   gcgcCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTgg                                                                                               <  1:435472/150‑1 (MQ=255)
                                                    cgcCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTgg                                                                                               <  1:1231283/149‑1 (MQ=255)
                                                      ccATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCt                                                                                           >  1:599880/1‑151 (MQ=255)
                                                             ctggctggTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGAtgtg                                                                                    >  2:1092301/1‑151 (MQ=255)
                                                               ggctggTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGaa                                                                                  <  1:556669/151‑1 (MQ=255)
                                                               ggctggTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGa                                                                                   <  1:331016/150‑1 (MQ=255)
                                                                 ctggTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGa                                                                                   <  1:839414/148‑1 (MQ=255)
                                                                  tggTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAgg                                                                               >  3:251697/1‑151 (MQ=255)
                                                                   ggTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAg                                                                                <  1:14655/149‑1 (MQ=255)
                                                                    gTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCt                                                                             <  1:1552640/151‑1 (MQ=255)
                                                                        gcCCGGACCTTTTCACTCTCAAACCCTGTTTTGTGGCAGGGGAAAATTAGGGGGATTTTACCCCAGGCAAGACGATGGTTGATATATACGGGTGCCTGGGCAAGCCAGGCCGTTTACAGGTGGTATTTGATCTGGATGTGAAAGGCtta                                                                           >  1:230549/1‑148 (MQ=255)
                                                                        gcCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCca                                                                         <  4:344150/151‑1 (MQ=255)
                                                                         cccGGACCTGTTCACTCTCAAACCCTGTTTTGTGGGAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTcc                                                                          >  1:592173/1‑149 (MQ=255)
                                                                         cccGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGAt                                                                                       >  3:404218/1‑136 (MQ=255)
                                                                         cccGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGAt                                                                                       <  4:404198/136‑1 (MQ=255)
                                                                          ccGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTCCCTCAGGCACGACGATGGTTGATATCGACTGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCcagc                                                                       <  1:1098654/151‑1 (MQ=255)
                                                                          ccGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCcagc                                                                       >  1:1651125/1‑151 (MQ=255)
                                                                          ccGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCcagc                                                                       >  3:400668/1‑151 (MQ=255)
                                                                          ccGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCcag                                                                        >  1:526411/1‑150 (MQ=255)
                                                                           cGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCcagca                                                                      <  2:1689276/151‑1 (MQ=255)
                                                                            ggACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCcagcag                                                                     <  4:415262/151‑1 (MQ=255)
                                                                            ggACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCcagc                                                                       >  4:566340/1‑149 (MQ=255)
                                                                              aCCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAgtg                                                                   >  2:140992/1‑151 (MQ=255)
                                                                               ccTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAgtgg                                                                  >  1:1247564/1‑151 (MQ=255)
                                                                                cTGTTCACTCTCAAACCCTGTTTTGTGGCAGTTGTACCTCAGGTCGAATTTTCCTCAGGCACGAGGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCGAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAgtgg                                                                  <  1:955480/150‑1 (MQ=255)
                                                                                  gTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAgtgggtg                                                               <  2:1014969/151‑1 (MQ=255)
                                                                                  gTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAgtgggt                                                                <  2:1090396/150‑1 (MQ=255)
                                                                                       ctctCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTg                                                           >  1:1556394/1‑150 (MQ=255)
                                                                                           cAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAg                                                         <  2:869668/148‑1 (MQ=255)
                                                                                           cAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAg                                                         <  1:1098218/148‑1 (MQ=255)
                                                                                           cAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGt                                                       <  1:536335/150‑1 (MQ=255)
                                                                                           cAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTg                                                      >  1:1191563/1‑151 (MQ=255)
                                                                                           cAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGGGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTg                                                      >  2:619522/1‑151 (MQ=255)
                                                                                           cAAACCCTGGTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAg                                                         <  2:885958/148‑1 (MQ=255)
                                                                                            aaaCCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTg                                                      <  1:464726/150‑1 (MQ=255)
                                                                                               cccTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCt                                                    >  2:1008215/1‑149 (MQ=255)
                                                                                                 cTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAAGGTACAGGGGGCATTGGATCTGGATGGGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTg                                                   >  1:182457/1‑148 (MQ=255)
                                                                                                  tGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCCGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGGTGTGAAAGGCTTCCCGCAGTGGGTGGCGGGGGTGCTGGct                                                >  1:287272/1‑150 (MQ=255)
                                                                                                   gTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGctct                                              >  2:158763/1‑151 (MQ=255)
                                                                                                   gTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGctct                                              <  2:101644/151‑1 (MQ=255)
                                                                                                     tttGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTgg                                            >  2:10382/1‑151 (MQ=255)
                                                                                                        gtgGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGcg                                          <  1:94772/150‑1 (MQ=255)
                                                                                                            cAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGcgtcgt                                      >  2:377702/1‑150 (MQ=255)
                                                                                                             aGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTaa                                    <  3:393871/151‑1 (MQ=255)
                                                                                                             aGTGGAAAATCAGCGCGAATTTACCTCAGGCACGACGATTGTTTATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTa                                     <  1:1638429/150‑1 (MQ=255)
                                                                                                                    aCTCAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTca                             <  2:1724498/151‑1 (MQ=255)
                                                                                                                       cAGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACAt                          <  2:852014/151‑1 (MQ=255)
                                                                                                                        aGGGCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGAGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATg                         >  3:271044/1‑151 (MQ=255)
                                                                                                                         gggCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGt                        <  2:915298/151‑1 (MQ=255)
                                                                                                                          ggCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGtt                       <  2:662266/151‑1 (MQ=255)
                                                                                                                          ggCGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGtt                       <  4:46472/151‑1 (MQ=255)
                                                                                                                            cGAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTAt                     >  2:1726999/1‑151 (MQ=255)
                                                                                                                             gAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTAtt                    <  1:1028881/151‑1 (MQ=255)
                                                                                                                             gAATTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTAtt                    >  2:803933/1‑151 (MQ=255)
                                                                                                                              aaTTTACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTg                   >  1:136200/1‑151 (MQ=255)
                                                                                                                                 ttACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGc                 >  1:191800/1‑150 (MQ=255)
                                                                                                                                 ttACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGAGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGc                 >  4:426888/1‑150 (MQ=255)
                                                                                                                                  tACTTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGCGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGATATTGGCa                >  1:1531341/1‑150 (MQ=255)
                                                                                                                                  tACCTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCAt               >  1:1120018/1‑151 (MQ=255)
                                                                                                                                    ccTCAGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATg              <  1:1015493/150‑1 (MQ=255)
                                                                                                                                     cTCAGGCACGACGATGGTTGATATCGAAGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGAACTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGTGGTGCTGGCTCTTGCGTCGTAACCTGTCACATGCTATTGGCATGc             >  1:114378/1‑150 (MQ=255)
                                                                                                                                        aGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtc         >  1:323616/1‑151 (MQ=255)
                                                                                                                                        aGGCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGGGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtc         >  1:706743/1‑151 (MQ=255)
                                                                                                                                         ggCACTACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCATCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGGTTCGAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACAGGTTATTGTCATGCAGtc         <  1:975982/150‑1 (MQ=255)
                                                                                                                                         ggCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtca        <  2:1117698/151‑1 (MQ=255)
                                                                                                                                          gCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtcat       <  1:1312451/151‑1 (MQ=255)
                                                                                                                                          gCACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtcat       >  2:851237/1‑151 (MQ=255)
                                                                                                                                           cACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGCAAGGCTTCAAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtcctt      >  1:1302069/1‑151 (MQ=255)
                                                                                                                                           cACGACGATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtcatt      >  2:612249/1‑151 (MQ=255)
                                                                                                                                             cgacgaTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtcattca    <  1:401492/151‑1 (MQ=255)
                                                                                                                                              gacgATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtcattcat   >  2:537315/1‑151 (MQ=255)
                                                                                                                                               acgATGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATc  <  2:467741/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATC  >  NC_000913/27948‑28241

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: