Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,492,667 C→T V229M (GTG→ATG)  idnO ← 5‑keto‑D‑gluconate‑5‑reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,492,6670CT98.4% 241.2 / NA 63V229M (GTG→ATG) idnO5‑keto‑D‑gluconate‑5‑reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/28);  minor base A (0/1);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATT  >  NC_000913/4492519‑4492814
                                                                                                                                                    |                                                                                                                                                   
aTGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCg                                                                                                                                                   <  1:533608/151‑1 (MQ=255)
 tGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCgc                                                                                                                                                  >  1:872901/1‑151 (MQ=255)
  ggCATAGTCCTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAGCAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCgca                                                                                                                                                 <  2:784213/151‑1 (MQ=255)
           tttATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATc                                                                                                                                        >  2:367293/1‑151 (MQ=255)
           tttATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATc                                                                                                                                        >  2:592559/1‑151 (MQ=255)
            ttATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATc                                                                                                                                        <  1:283626/150‑1 (MQ=255)
                 tctTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTcc                                                                                                                                  >  1:880949/1‑151 (MQ=255)
                  ctTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCt                                                                                                                                 >  2:205825/1‑151 (MQ=255)
                    tattaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGc                                                                                                                               >  2:44154/1‑151 (MQ=255)
                     attaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGc                                                                                                                               <  1:322752/150‑1 (MQ=255)
                     attaGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGc                                                                                                                               <  2:780121/150‑1 (MQ=255)
                        aGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGAt                                                                                                                           <  1:644191/151‑1 (MQ=255)
                        aGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGAt                                                                                                                           >  2:611935/1‑151 (MQ=255)
                        aGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGAt                                                                                                                           <  1:156093/151‑1 (MQ=255)
                        aGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGAt                                                                                                                           >  1:11446/1‑151 (MQ=255)
                                  aaGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCg                                                                                                                 >  1:166311/1‑151 (MQ=255)
                                     cagcGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTgc                                                                                                              >  1:131526/1‑151 (MQ=255)
                                        cgcTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTgcggc                                                                                                           <  1:688355/151‑1 (MQ=255)
                                        cgcTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGGGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTgcggc                                                                                                           >  2:185832/1‑151 (MQ=255)
                                             tCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCggggg                                                                                                      >  2:325187/1‑151 (MQ=255)
                                                cgcTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTcc                                                                                                   >  2:1003521/1‑151 (MQ=255)
                                                  cTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTcc                                                                                                   <  1:812179/149‑1 (MQ=255)
                                                   tGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGtt                                                                                                >  2:1013157/1‑151 (MQ=255)
                                                     cACGTAGGGGGGCTGTTAAACAGCCGCTAACATGCCGCCATCAACAGACAACAGGTGGCCGTTTACGAAATCAGTGGCTTTTGAAGAAAGGAACATCGCATCACCAGTCAGTTCCTGCCGATCTCCCCAGCATGCGGCGGGGGTCCGttt                                                                                               <  2:869205/150‑1 (MQ=255)
                                                           ggggggCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAAc                                                                                        <  1:936284/151‑1 (MQ=255)
                                                            gggggCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAAcc                                                                                       <  2:373946/151‑1 (MQ=255)
                                                              gggCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAAcc                                                                                       >  1:462489/1‑149 (MQ=255)
                                                              gggCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAg                                                                                     <  2:10117/151‑1 (MQ=255)
                                                              gggCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAg                                                                                     <  2:859286/151‑1 (MQ=255)
                                                                     tAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGt                                                                                >  1:143936/1‑149 (MQ=255)
                                                                      aaaCAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCgg                                                                                 >  1:717443/1‑147 (MQ=255)
                                                                      aaaCAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGa                                                                              >  1:845977/1‑150 (MQ=255)
                                                                         cAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAgg                                                                           >  2:901533/1‑150 (MQ=255)
                                                                                 aaCATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTc                                                                    <  2:377703/149‑1 (MQ=255)
                                                                                 aaCATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGGGAAGGCTTCGTcc                                                                   >  2:170005/1‑150 (MQ=255)
                                                                                  aCATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTc                                                                 <  1:50090/151‑1 (MQ=255)
                                                                                   cATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAATGCTTCGTCCTCa                                                                <  1:1027446/151‑1 (MQ=255)
                                                                                      gccgccATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAAcc                                                             <  1:260219/151‑1 (MQ=255)
                                                                                      gccgccATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAAc                                                              <  2:92292/150‑1 (MQ=255)
                                                                                        cgccATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCCTCGGCCTCAAcc                                                             >  2:784730/1‑149 (MQ=255)
                                                                                         gccATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGt                                                          <  1:660923/151‑1 (MQ=255)
                                                                                               aacaaacaACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTa                                                    <  1:548105/151‑1 (MQ=255)
                                                                                                acaaacaaCAGGTGGCCGGTTACGAAATCAGAGGCTTTTGAAGAAAGGAACAACGCAGCACCAATCAGTTCCTGGGGATCTCCCCAGCGTGCGGCGGGGGTTCGGTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAg                                                   <  2:1034856/151‑1 (MQ=255)
                                                                                                  aaacaaCAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTc                                                 <  1:175783/151‑1 (MQ=255)
                                                                                                  aaacaaCAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAATGCTTTAGTc                                                 >  1:252259/1‑151 (MQ=255)
                                                                                                           gTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCCTCCTCCTCAACCAGTGCTTTAATCATTTTTGtt                                        >  2:781722/1‑151 (MQ=255)
                                                                                                             ggCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGGGAAGGCCTCGTCCTCAACCAGTGCTTTTGTCATTTCTGtttt                                      >  2:292029/1‑151 (MQ=255)
                                                                                                                       gAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGccc                            <  2:384752/151‑1 (MQ=255)
                                                                                                                       gAAATCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGccc                            >  2:725700/1‑151 (MQ=255)
                                                                                                                           tCAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGGCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCgg                          >  2:389281/1‑149 (MQ=255)
                                                                                                                            cAGAGGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGgcg                        <  1:617285/150‑1 (MQ=255)
                                                                                                                               agGCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCaa                     >  1:348593/1‑150 (MQ=255)
                                                                                                                                 gCTTTTGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATAc                  <  1:1002810/151‑1 (MQ=255)
                                                                                                                                   ttttGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGGGCAATAcc                 >  1:159115/1‑150 (MQ=255)
                                                                                                                                   ttttGAAGAAAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGGGCTTTAGTCATTTCTGTTTTTAAATAGCCCGGCGCAATAcc                 >  1:861190/1‑150 (MQ=255)
                                                                                                                                        aagaaAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGAc           <  2:720299/151‑1 (MQ=255)
                                                                                                                                        aagaaAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGAc           >  1:936477/1‑151 (MQ=255)
                                                                                                                                        aagaaAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCCGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGAc           >  1:239836/1‑151 (MQ=255)
                                                                                                                                         agaaAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGATTGCACAACCAGGAGGTGAAGGCCTCGGCCTCAAACAGTGATTTAGACATTTTTGTTTTGAAAAAGACCGGGGCAATACTGTTGaac          >  2:233140/1‑149 (MQ=255)
                                                                                                                                          gaaAGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGAcca         <  2:940519/151‑2 (MQ=255)
                                                                                                                                           aaaGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCCTCGGCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTg        >  2:5233/1‑151 (MQ=255)
                                                                                                                                            aaGGAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGGGGAGGCCCCGACCTCAAACAGGGCCTTAGTCCTTTCTGTTTTTCAACAGCCCCGCGCAATACCGtttgc           >  2:1034516/1‑144 (MQ=255)
                                                                                                                                              ggAACATCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAat     <  2:994715/151‑1 (MQ=255)
                                                                                                                                                 acatCGCAGCCCCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGCCCTGAatat   <  2:551924/150‑1 (MQ=255)
                                                                                                                                                 acatCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATAtt  >  1:639909/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGTAGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGGCGGGGGTCCGTTTGCACAACCAGGCGGTGAAGGCTTCGTCCTCAACCAGTGCTTTAGTCATTTCTGTTTTGAAATAGCCCGGCGCAATACCGTTGACCTGAATATT  >  NC_000913/4492519‑4492814

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