Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 321,039 G→A D70N (GAC→AAC)  ykgE → putative dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913321,0390GA98.4% 648.4 / ‑2.1 183D70N (GAC→AAC) ykgEputative dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (86/94);  total (88/95)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 192 total aligned reads displayed.

GGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCATCA  >  NC_000913/320949‑321128
                                                                                          |                                                                                         
  aTGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAAcaacga                                                                                       <  1:4595415/93‑1 (MQ=255)
  aTGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAAcaacga                                                                                       <  1:1577676/93‑1 (MQ=255)
   tgctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGAt                                                                                      >  1:3452748/1‑93 (MQ=255)
   tgctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGAt                                                                                      >  1:484095/1‑93 (MQ=255)
    gctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATc                                                                                     <  1:3244071/93‑1 (MQ=255)
    gctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATc                                                                                     >  2:3560086/1‑93 (MQ=255)
     ctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATcc                                                                                    <  2:3075249/93‑1 (MQ=255)
      tgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATccc                                                                                   <  2:109071/93‑1 (MQ=255)
        cGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCat                                                                                 <  2:2193927/93‑1 (MQ=255)
          gacAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCattt                                                                               >  1:1093629/3‑92 (MQ=255)
            cAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCattatt                                                                             <  2:2147771/93‑1 (MQ=255)
            cAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCattatt                                                                             <  1:4199355/93‑1 (MQ=255)
              gCCTGCGAGCAATAGGGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTc                                                                           >  2:2789715/1‑93 (MQ=255)
                 tGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCAcc                                                                        <  2:1841539/93‑1 (MQ=255)
                 tGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCAcc                                                                        <  2:1759481/93‑1 (MQ=255)
                 tGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCAcc                                                                        <  1:1999275/93‑1 (MQ=255)
                 tGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCAcc                                                                        <  2:858234/93‑1 (MQ=255)
                  gCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCg                                                                       >  2:4310343/1‑93 (MQ=255)
                  gCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCg                                                                       >  1:2990546/1‑93 (MQ=255)
                   cGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCgg                                                                      >  2:3036541/1‑93 (MQ=255)
                    gATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggc                                                                     >  2:1261641/1‑93 (MQ=255)
                     aTCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggct                                                                    <  2:1430898/93‑1 (MQ=255)
                     aTCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggct                                                                    <  2:4506870/93‑1 (MQ=255)
                     aTCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggct                                                                    <  1:2017256/93‑1 (MQ=255)
                       cAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCgggaaa                                                                  <  2:4005292/93‑5 (MQ=255)
                       cAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctgg                                                                  <  2:1826003/93‑1 (MQ=255)
                        aaTAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctggc                                                                 <  2:4334054/93‑1 (MQ=255)
                        aaTAGCGGTTATATCAAAGAAGCCATTCCAGCGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctggc                                                                 >  2:875943/1‑93 (MQ=255)
                         aTAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctggct                                                                >  2:608859/1‑93 (MQ=255)
                         aTAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctggct                                                                >  2:2983417/1‑93 (MQ=255)
                          tAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTCCACCGGCTGGctc                                                               >  2:3179047/1‑93 (MQ=255)
                           aGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGctct                                                              <  2:3154639/93‑1 (MQ=255)
                           aGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGctct                                                              <  1:2234478/93‑1 (MQ=255)
                                 tataTCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCAcc                                                        >  1:2076585/1‑93 (MQ=255)
                                  atatCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCt                                                       >  2:2676699/1‑93 (MQ=255)
                                   tatCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTa                                                      <  1:378576/93‑1 (MQ=255)
                                   tatCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTa                                                      <  1:2449634/93‑1 (MQ=255)
                                    atCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAt                                                     >  1:4053503/1‑93 (MQ=255)
                                    atCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAt                                                     <  1:260374/93‑1 (MQ=255)
                                       aaaGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAtcac                                                  >  1:3681859/1‑90 (MQ=255)
                                       aaaGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGcc                                                  >  2:2568938/1‑93 (MQ=255)
                                       aaaGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGcc                                                  >  1:1630202/1‑93 (MQ=255)
                                       aaaGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGcc                                                  <  2:2933028/93‑1 (MQ=255)
                                       aaaGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGcc                                                  >  1:90816/1‑93 (MQ=255)
                                          gaagCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTa                                               >  1:2339875/1‑93 (MQ=255)
                                           aagCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaa                                              >  2:2288281/1‑93 (MQ=255)
                                            agCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaaa                                             <  2:70191/93‑1 (MQ=255)
                                            agCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaaa                                             >  1:3872684/1‑93 (MQ=255)
                                               gATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAg                                          <  1:2135444/93‑1 (MQ=255)
                                                aTTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGt                                         >  1:1800161/1‑93 (MQ=255)
                                                aTTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGt                                         >  1:2743059/1‑93 (MQ=255)
                                                  tCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTa                                       <  2:1121003/93‑1 (MQ=255)
                                                  tCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGCTa                                       >  1:4155070/1‑93 (MQ=255)
                                                  tCCAGGGATGAAAAAGCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATTTCACCGGCTGGCTCTTGCGCCTATGCCGTAAAAAGTTa                                       <  2:535029/93‑1 (MQ=255)
                                                   ccAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAc                                      >  1:1341096/1‑93 (MQ=255)
                                                    cAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAcc                                     <  1:1473642/93‑1 (MQ=255)
                                                    cAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAcc                                     <  1:2184734/93‑1 (MQ=255)
                                                    cAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAcc                                     <  2:1395918/93‑1 (MQ=255)
                                                     aGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAccc                                    <  1:4686243/93‑1 (MQ=255)
                                                       ggATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATTCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGCTACCCGa                                  >  1:1838316/1‑93 (MQ=255)
                                                         aTGAAAAATCTGATCGCGGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCTCCTATGCCGTAAAAAGTTCCCCGACg                                <  1:893438/93‑1 (MQ=255)
                                                          tGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGt                               >  2:3626982/1‑93 (MQ=255)
                                                          tGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGt                               <  1:4411753/93‑1 (MQ=255)
                                                           gAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTa                              <  2:503093/93‑1 (MQ=255)
                                                           gAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTa                              <  2:4046609/93‑1 (MQ=255)
                                                           gAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGCCTCCTGCACCTATGCCGTACACACTTACCCGACGTa                              >  1:2622709/1‑93 (MQ=255)
                                                             aaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATc                            <  2:4139177/93‑1 (MQ=255)
                                                             aaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATc                            <  1:3436216/93‑1 (MQ=255)
                                                             aaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATc                            <  1:359430/93‑1 (MQ=255)
                                                             aaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATc                            >  2:2597206/1‑93 (MQ=255)
                                                              aaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCt                           <  2:4279868/93‑1 (MQ=255)
                                                              aaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCt                           <  1:4285808/93‑1 (MQ=255)
                                                              aaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCt                           <  2:4429581/93‑1 (MQ=255)
                                                               aaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTg                          <  1:3181106/93‑1 (MQ=255)
                                                               aaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTg                          <  2:1981706/93‑1 (MQ=255)
                                                                 tCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGc                        <  1:4650329/93‑1 (MQ=255)
                                                                   tGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCgg                      <  1:4546708/93‑1 (MQ=255)
                                                                    gATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGa                     <  1:736175/93‑1 (MQ=255)
                                                                     aTCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGAt                    <  1:1057513/93‑1 (MQ=255)
                                                                     aTCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGAt                    >  2:4067337/1‑93 (MQ=255)
                                                                      tCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACGCGACGGATCTGGCGGATg                   >  2:4179496/1‑93 (MQ=255)
                                                                       cgccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGa                  <  1:432810/93‑1 (MQ=255)
                                                                         ccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAAc                >  2:2037428/1‑93 (MQ=255)
                                                                          cgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAAcc               <  2:758590/93‑1 (MQ=255)
                                                                           gcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCt              >  2:2408877/1‑93 (MQ=255)
                                                                           gcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCt              <  2:3779387/93‑1 (MQ=255)
                                                                            cACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTg             <  1:549455/93‑1 (MQ=255)
                                                                             aCTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAGCCTGa            >  2:2425902/1‑93 (MQ=255)
                                                                              cTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACTCGACGTATCTGGCTGATGAACCTGaa           >  1:2826414/1‑93 (MQ=255)
                                                                               tGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAAt          >  1:4278192/1‑93 (MQ=255)
                                                                               tGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAAt          >  1:3716622/1‑93 (MQ=255)
                                                                                ggaggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATg         >  1:3094233/1‑93 (MQ=255)
                                                                                ggaggaTAACAACGATCCCATTAGTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATg         <  2:3381642/93‑1 (MQ=255)
                                                                                  aggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATggg       >  2:1328582/1‑93 (MQ=255)
                                                                                   ggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGc      <  1:821027/93‑1 (MQ=255)
                                                                                   ggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGc      >  2:4248467/1‑93 (MQ=255)
                                                                                    gaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCa     <  2:3744769/93‑1 (MQ=255)
                                                                                     aTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCAt    <  2:155848/93‑1 (MQ=255)
                                                                                     aTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCAt    >  2:3094976/1‑93 (MQ=255)
                                                                                       aacaacGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGCTGAACCTGAATGGGCATCa  >  2:1937276/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCATCA  >  NC_000913/320949‑321128

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: