Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,150,120 T→C 16.2% F153F (TTT→TTC yghA → putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,150,1200TC16.2% 191.4 / 27.8 88F153F (TTT→TTCyghAputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (32/41);  new base C (6/8);  total (39/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CGAAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCA  >  NC_000913/3150021‑3150219
                                                                                                   |                                                                                                   
cGAAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGttt                                                                                                     >  2:241136/1‑100 (MQ=255)
 gAAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGttt                                                                                                     >  1:1285857/1‑99 (MQ=255)
 gAAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTc                                                                                                    <  2:2041406/100‑1 (MQ=255)
  aaGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTACAGATATTAAAGACCTCACCAACAAACAGTTTCa                                                                                                   >  1:381989/1‑100 (MQ=255)
     cgcACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCaaaa                                                                                                >  1:1880038/1‑100 (MQ=255)
       cacaAGGCGTTAGGCGGGCTTGATATTATGGAGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGa                                                                                              <  2:1419224/100‑1 (MQ=255)
         caAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGAcc                                                                                            <  2:1288026/100‑1 (MQ=255)
          aaGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCt                                                                                           <  2:2014988/100‑1 (MQ=255)
               gTTAGGCGGGCTGGATATTATGTCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTGGcc                                                                                      <  2:166131/100‑1 (MQ=255)
               gTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGcc                                                                                      <  1:890087/100‑1 (MQ=255)
                   ggCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTa                                                                                  >  1:185189/1‑100 (MQ=255)
                   ggCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTa                                                                                  <  2:672980/100‑1 (MQ=255)
                      gggCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAAc                                                                                <  1:902253/99‑1 (MQ=255)
                      gggCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACg                                                                               <  1:772770/100‑1 (MQ=255)
                      gggCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACg                                                                               <  2:2124679/100‑1 (MQ=255)
                         cTGGATATTATGGCCCTGGTTGCCCTGCAACAGGTTTCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACg                                                                               <  2:858834/97‑1 (MQ=255)
                            gATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTcg                                                                         <  2:1254857/100‑1 (MQ=255)
                              tattatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTcg                                                                         <  1:633834/98‑1 (MQ=255)
                                 tatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTg                                                                    <  1:1880456/100‑1 (MQ=255)
                                 tatGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTACAAAAGCCCTTTGCCATTAAAGTTTTCGCGCt                                                                     >  2:1108202/1‑99 (MQ=255)
                                  atGGCGCTGGTCGCCGGAAAACAAGATGCCATTCAGGATATTGCAAACCTGACCAGCGAACAGTTTCACAAGACCTTTGCGATTAACGTTGTCGCGCTGt                                                                   >  2:1652539/1‑100 (MQ=255)
                                    ggCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTc                                                                 >  2:2011337/1‑100 (MQ=255)
                                    ggCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTc                                                                 <  1:365578/100‑1 (MQ=255)
                                       gcTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTgg                                                              <  2:306140/100‑1 (MQ=255)
                                        cTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTgcc                                                             <  1:1627280/100‑3 (MQ=255)
                                        cTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGc                                                             <  1:1798860/100‑1 (MQ=255)
                                        cTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGc                                                             <  2:1076163/100‑1 (MQ=255)
                                             cGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAAcc                                                        >  1:1913160/1‑100 (MQ=255)
                                             cGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAAcc                                                        >  2:1848273/1‑100 (MQ=255)
                                              gCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAAcc                                                        >  1:1241973/1‑99 (MQ=255)
                                               ccGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCa                                                      <  1:2391086/100‑1 (MQ=255)
                                               ccGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCa                                                      >  1:1467806/1‑100 (MQ=255)
                                                  ggAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGa                                                   >  1:1536122/1‑100 (MQ=255)
                                                  ggAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAAGa                                                   <  2:1101892/100‑1 (MQ=255)
                                                    aaaCAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAg                                                 >  1:1713332/1‑100 (MQ=255)
                                                     aaCAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGc                                                <  1:1266314/100‑1 (MQ=255)
                                                      aCAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCg                                               <  1:252603/100‑1 (MQ=255)
                                                       cAGGTTGCCATTCCGGATATTGCAGACATCACCAGCGAAAAGTTTCAAAAGACATTTGCCATTAAAGTTTTAGCGATGTTCTGGCTAACCTAGGAAGCGa                                              >  1:554247/1‑100 (MQ=255)
                                                        aGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGAt                                             <  1:2248739/100‑1 (MQ=255)
                                                         ggTTGCCATTCCAGATATTACAGACGTAACCAGCGAACAGTGTCAAAAGACCTTTGCCTTAAACGTTTACGCGCTGTTCTGACTAACCCAGGAACCg                                               >  2:1744475/1‑97 (MQ=255)
                                                          gTTGCCATTCCGTATATTGTAGATCTTATCAGCGAACAGTTTCAAAAGACCTTTGTCATTAACGTTTTCGTGCTGTTCTTGCTAACCCAGGAAGCGATcc                                           <  2:866164/100‑1 (MQ=255)
                                                          gTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATcc                                           <  1:554909/100‑1 (MQ=255)
                                                            tGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATccc                                          >  1:2003619/1‑99 (MQ=255)
                                                             gCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATccccc                                        <  2:2195974/100‑1 (MQ=255)
                                                             gCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATccccc                                        <  2:473438/100‑1 (MQ=255)
                                                             gCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATccccc                                        <  1:2011359/100‑1 (MQ=255)
                                                             gCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATacccc                                        >  1:135510/1‑100 (MQ=255)
                                                              ccATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCt                                       >  1:651130/1‑100 (MQ=255)
                                                              ccATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCt                                       >  1:2363752/1‑100 (MQ=255)
                                                              ccATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTCTTCGCGCTGTTCTGGCTAACCCAGCAAGCGATCCCCCt                                       <  2:1474736/100‑1 (MQ=255)
                                                               cATTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTg                                      <  2:1011717/100‑1 (MQ=255)
                                                                aTTCCGGTTATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTgn                                     <  1:1057109/100‑2 (MQ=255)
                                                                aTTCCGGATATTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGc                                     >  2:1775615/1‑100 (MQ=255)
                                                                 ttCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCAATCCCCCTACt                                    >  1:599627/1‑100 (MQ=255)
                                                                  tCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTa                                   <  1:2173440/100‑1 (MQ=255)
                                                                  tCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTAGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTa                                   >  2:392479/1‑100 (MQ=255)
                                                                    cGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAcc                                 >  1:2211419/1‑100 (MQ=255)
                                                                    cGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAcc                                 <  1:138817/100‑1 (MQ=255)
                                                                    cGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTAcc                                 <  1:2275822/100‑1 (MQ=255)
                                                                     ggATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCGCCTGCTACCg                                >  2:703385/1‑100 (MQ=255)
                                                                        tatTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGaaa                             <  2:1985432/100‑1 (MQ=255)
                                                                        tatTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGaaa                             >  1:1323476/1‑100 (MQ=255)
                                                                        tatTGCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGaaa                             >  1:2309834/1‑100 (MQ=255)
                                                                         atTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAg                            >  1:2059273/1‑100 (MQ=255)
                                                                          ttGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAgg                           <  1:700603/100‑1 (MQ=255)
                                                                           tGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGt                          >  1:2212025/1‑100 (MQ=255)
                                                                            gCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAgg                           >  2:1451851/1‑98 (MQ=255)
                                                                            gCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTg                         >  1:561580/1‑100 (MQ=255)
                                                                            gCAGACCTCACCAGCGAACAGTTTAAAAAGACCTTTGCCATCAACGTTTTCGCGCTGTTCTGGCTAACCCAGGACACGATCCCCCTGCTACCGAAAg                            >  2:1078892/1‑97 (MQ=255)
                                                                            gCAGACCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGCTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAgg                           >  2:602362/1‑98 (MQ=255)
                                                                              aGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCa                       <  2:2091871/100‑1 (MQ=255)
                                                                               gACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCaa                      >  2:1715482/1‑100 (MQ=255)
                                                                                aCCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                     <  2:2138631/100‑1 (MQ=255)
                                                                                aCCTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAg                     <  2:1783711/100‑1 (MQ=255)
                                                                                  cTCACCCGCGAACAGTTTCAAAAGCCCTTTGCCGCTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGt                    <  2:2019996/99‑1 (MQ=255)
                                                                                  cTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTa                   <  1:682523/100‑1 (MQ=255)
                                                                                  cTCACCAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTa                   >  2:1837504/1‑100 (MQ=255)
                                                                                   tCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTa                   <  1:926829/99‑1 (MQ=255)
                                                                                      ccAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCACCCTGCTACCGAAAGGTGCAAGTatcat               >  2:80543/1‑100 (MQ=255)
                                                                                      ccAGCGAACAGTTCCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatcat               >  2:395279/1‑100 (MQ=255)
                                                                                        agcgctcAGTTTCAAAAGAGTTGTGCCATTAACGTTTTCGCGCTTTTCTGGCTAACCCAGGAAGTGATCCCTCTGCTACCGAAAGGt                          <  2:2029802/81‑1 (MQ=255)
                                                                                        aGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatcatc              <  1:1401775/99‑1 (MQ=255)
                                                                                         ggcaaCAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTTTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTatc                 <  1:1605091/92‑1 (MQ=255)
                                                                                         gCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACAGAAAGGTGCAAGGatcatc              <  2:1285844/98‑1 (MQ=255)
                                                                                          cGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacc           <  1:1702911/100‑1 (MQ=255)
                                                                                          cGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacc           >  2:2207794/1‑100 (MQ=255)
                                                                                          cGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacc           >  2:1000137/1‑100 (MQ=255)
                                                                                           gAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacca          >  2:573645/1‑100 (MQ=255)
                                                                                           gAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATcacca          >  1:1920959/1‑100 (MQ=255)
                                                                                               aGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTc      <  1:399811/100‑1 (MQ=255)
                                                                                                 tttCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCg     >  1:211570/1‑99 (MQ=255)
                                                                                                  ttCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTc   <  1:613664/100‑1 (MQ=255)
                                                                                                  ttCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTc   <  2:561577/100‑1 (MQ=255)
                                                                                                  ttCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTc   <  1:457404/100‑1 (MQ=255)
                                                                                                  ttCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTc   <  1:1836985/100‑1 (MQ=255)
                                                                                                  ttCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTc   <  2:1301607/100‑1 (MQ=255)
                                                                                                  ttCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTc   <  1:2089525/100‑1 (MQ=255)
                                                                                                   tCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTc   <  2:662351/99‑1 (MQ=255)
                                                                                                   tCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCa  <  1:708469/100‑1 (MQ=255)
                                                                                                   tCAAAAGACCTTTGCCAGTAACGTTTGCGCGCTGTTCGGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCa  <  2:851244/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGAAGCGCACAAGGCGTTAGGCGGGCTGGATATTATGGCGCTGGTCGCCGGGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCA  >  NC_000913/3150021‑3150219

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