Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,932,094 C→G 16.1% intergenic (+35/+21) purT → / ← eda phosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,932,0940CG16.1% 243.1 / 39.2 95intergenic (+35/+21)purT/edaphosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase
Reads supporting (aligned to +/- strand):  ref base C (42/37);  new base G (4/11);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 9.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

GAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTG  >  NC_000913/1931996‑1932191
                                                                                                  |                                                                                                 
gagTGTTGTTGACGCCATTGAACGCTCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTGAACCCCGGCAAAAAAACGCTACAAAAATGCCCGATcc                                                                                                   <  1:2118208/99‑1 (MQ=255)
gagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCt                                                                                                  >  1:2478850/1‑100 (MQ=255)
 agTGTTGTTGACGCCATTGAACGCGCAAAGCACGCCGCCGGACAGATAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCAACAAAAATGCCCGATCCTc                                                                                                 <  1:1786946/100‑1 (MQ=255)
  gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCg                                                                                                >  2:28453/1‑100 (MQ=255)
     ttgttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATc                                                                                             >  2:1169622/1‑100 (MQ=255)
       gttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCgg                                                                                           >  2:1490823/1‑100 (MQ=255)
        ttgACGCCATTGAACGCGCGAAGCATGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCggg                                                                                          <  2:722334/100‑1 (MQ=255)
         tgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCgggn                                                                                         <  2:72445/100‑2 (MQ=255)
          gACGGCATGGAACGCGCGAAGCACTCCGCCGGACAGGTAAAAGTACAGGGTTAAAGCCGGCAAAAAAACGCTACAATAATGCCCGATCGAGGATCGGGCa                                                                                        <  1:1142515/100‑1 (MQ=255)
          gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                        <  2:730650/100‑1 (MQ=255)
          gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                        <  1:1039702/100‑1 (MQ=255)
          gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                        <  1:1719904/100‑1 (MQ=255)
          gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                                                        <  2:1183631/100‑1 (MQ=255)
          gACGCCATTGAACGCGCGAAGCACGCCGACGGACAGGTAAACGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                        <  1:115702/100‑1 (MQ=255)
           aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                        <  2:2193577/99‑1 (MQ=39)
           aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAt                                                                                       <  2:282192/100‑1 (MQ=255)
           aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGATAAACCAGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAt                                                                                       >  2:270782/1‑100 (MQ=255)
            cGCTATTGACCGCGCGCACCCCGCCGCCAGCTAGATAAAAGTACAGGCTTACACCCGGCAAAAAACCGCTACCAAAATGACCGATCCTCGATCGGGc                                                                                         <  1:1677965/97‑1 (MQ=255)
            cGCCATTGAACGCGGGAGGCACGCCGCCGGAAAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAATTGCCCGATCGAGGATCGGGCAt                                                                                       <  1:783538/99‑1 (MQ=255)
            cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAtt                                                                                      <  2:1187452/100‑1 (MQ=39)
             gCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAttt                                                                                     <  1:587840/100‑1 (MQ=255)
               cATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAttttt                                                                                   <  2:2289165/100‑2 (MQ=38)
               cATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTg                                                                                   <  1:1117690/100‑1 (MQ=255)
                atggAACGCGCGAAGCGCGCCGCCGGACAGGGAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAACAATGCCCGATCGAGGATCGGGCAtttttg                                                                                  <  1:1977102/97‑3 (MQ=37)
                aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAtttttg                                                                                  <  2:1780023/100‑3 (MQ=38)
                  tGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAAAGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGAc                                                                                 >  2:1509498/1‑99 (MQ=255)
                   gAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtt                                                                               >  2:2023881/1‑100 (MQ=255)
                    accgcgccTAGCACGCAGCCGCACAGGTAAACGTACACGCTAAAACAAGGCAACAAAACGCTACAAAAATGCCCTATCCTCGATCGGGCATTTTGACtt                                                                               <  1:1543587/97‑1 (MQ=255)
                    aaCGAGCAAAGTACGACGCAGATCAGGGAAAGCAACATGGTTAAAGCCGACAAAAATACGCTTCAAAAATGCCCGGTCCTCGATCGGGCATTTTGACttt                                                                              <  1:1801443/100‑1 (MQ=255)
                      cgcgcgAAGCACGCCGCCGGACAGGTACAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGCTCCTCGATCGGGCATTTTGACTTTTa                                                                            >  2:1281420/1‑100 (MQ=255)
                      cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtttt                                                                             <  1:2442588/99‑1 (MQ=255)
                      cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                            >  2:113885/1‑100 (MQ=255)
                      cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                            >  2:264048/1‑100 (MQ=255)
                      cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                            <  2:2380419/100‑1 (MQ=255)
                      cgcgcgAAGCACGCAGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAAAGCCCGATACTCGATCGGGCATTTAGACTTTTa                                                                            >  1:2297881/1‑100 (MQ=255)
                       gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTAc                                                                           >  1:2250764/1‑100 (MQ=255)
                       gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTAc                                                                           >  2:1257689/1‑100 (MQ=255)
                         gcgAAGCACGCCGGCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGAGACAAAAATGCCCGATCCTCGATCGGGCATTTTGGCTTTTACa                                                                          >  1:518219/1‑99 (MQ=255)
                          cgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGc                                                                        >  1:288554/1‑100 (MQ=255)
                          cgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGc                                                                        >  1:1723009/1‑100 (MQ=255)
                          cgAAGCACGCCGCAGGACAGGTAAAAGCACAGGGTTAAACACGGCAAAAAAACGCTACAACACTGCCCGCTCCTCGATCGGGCATTTTGACTTTTACAGc                                                                        <  1:1509489/100‑1 (MQ=255)
                            aaGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCtt                                                                      >  2:2237240/1‑100 (MQ=255)
                            aaGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCtt                                                                      >  2:2446728/1‑100 (MQ=255)
                              gCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAg                                                                    >  1:1265028/1‑100 (MQ=255)
                               cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                                                   >  2:2222412/1‑100 (MQ=255)
                               cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                                                   >  1:832174/1‑100 (MQ=255)
                                aCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgcg                                                                  >  2:2301230/1‑100 (MQ=255)
                                 cgccgcCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgcgc                                                                 >  1:755381/1‑100 (MQ=255)
                                     gccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCtt                                                              <  1:2287104/99‑1 (MQ=255)
                                     gccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTc                                                             <  1:2475731/100‑1 (MQ=255)
                                      ccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCt                                                            >  1:1850827/1‑100 (MQ=255)
                                       cgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTa                                                           <  1:2393080/100‑1 (MQ=255)
                                        ggACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                          <  2:2250779/100‑1 (MQ=255)
                                         gACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                                                         <  2:2332132/100‑1 (MQ=255)
                                          aCAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAg                                                        <  1:586697/100‑1 (MQ=255)
                                           cAGATACAAGTAGAGGATTGATCCCGGAGAGAAAACGCTACAAAAATTCCCGATCCTCGTTCGTGCATTTTGACTTTTACAGCTTa                                                                     <  1:2493353/86‑1 (MQ=255)
                                             ggTAAAAGTACAGTGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                     <  2:184712/100‑1 (MQ=255)
                                             ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                     >  1:2565473/1‑100 (MQ=255)
                                             ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                     >  2:890070/1‑100 (MQ=255)
                                             ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                     >  1:221866/1‑100 (MQ=255)
                                              gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                                                    <  1:1281892/100‑1 (MQ=255)
                                               tAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCa                                                   >  2:291892/1‑100 (MQ=255)
                                                aaaaGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAc                                                  >  1:1864527/1‑100 (MQ=255)
                                                   aGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgcg                                               <  1:291889/100‑1 (MQ=255)
                                                    gTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACgcgc                                              >  1:755673/1‑100 (MQ=255)
                                                     tACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGACCGGGCACTTTGACCATTACCGCTTAGCGCCTTCTACAGACTCACCCGcc                                             >  1:2269390/1‑100 (MQ=255)
                                                      aCAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCa                                            >  2:1901732/1‑100 (MQ=255)
                                                       cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg                                           >  2:978504/1‑100 (MQ=255)
                                                       cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg                                           >  2:906269/1‑100 (MQ=255)
                                                       cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg                                           >  2:753714/1‑100 (MQ=255)
                                                         gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                         >  1:251918/1‑100 (MQ=255)
                                                          ggTTAAACCCGGCAAAAAAACGCTACAAATATGCCCGATCCTCGATTGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt                                        <  1:2347273/100‑1 (MQ=255)
                                                           gTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa                                       <  2:86918/100‑1 (MQ=255)
                                                               aaCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATCTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTaa                                   <  2:2461663/100‑1 (MQ=255)
                                                                  ccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTATCGCCTTCTACAGCTTCACGCCCCAGCTTAGTAATg                                 >  1:1085414/1‑99 (MQ=255)
                                                                  ccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGc                                <  2:44914/100‑1 (MQ=255)
                                                                    ggCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCgg                              <  1:843243/100‑1 (MQ=255)
                                                                    ggCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCgg                              <  2:1723018/100‑1 (MQ=255)
                                                                     gCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTATAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGt                             <  2:1810601/100‑1 (MQ=255)
                                                                      cAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTc                            >  2:1024901/1‑100 (MQ=255)
                                                                       aaaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTc                            <  1:1807566/99‑1 (MQ=255)
                                                                           aaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAAt                       >  2:1480376/1‑100 (MQ=255)
                                                                            aaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATc                      <  2:668147/100‑1 (MQ=255)
                                                                                cTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGccc                  <  1:1173224/100‑1 (MQ=255)
                                                                                 tACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                 <  2:1500909/100‑1 (MQ=255)
                                                                                  aCAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                >  1:231390/1‑100 (MQ=255)
                                                                                   cAAAAATGCCCGATCCTCGATCGGGCATCTTCACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGc                <  1:514463/99‑1 (MQ=255)
                                                                                    caaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCtt              >  2:1034210/2‑100 (MQ=38)
                                                                                     aaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc             >  2:2012378/1‑100 (MQ=38)
                                                                                     aaaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc             >  1:2378197/1‑100 (MQ=255)
                                                                                       aaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCa           <  2:2226189/100‑1 (MQ=255)
                                                                                        aTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg          >  2:880161/1‑100 (MQ=39)
                                                                                        aTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg          <  1:1281415/100‑1 (MQ=255)
                                                                                         tGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgc         >  2:1886242/1‑100 (MQ=255)
                                                                                            ccGATCTTCGATCGGGCATTTTGACTTTTACAGCTTTGCGCCTTCTACAGCTTCACGAGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa      <  1:1042410/100‑1 (MQ=255)
                                                                                            ccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa      >  2:752588/1‑100 (MQ=255)
                                                                                                 ccTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTATAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCTCGCTTCCAGCGCATCTg  <  1:131802/99‑1 (MQ=255)
                                                                                                  |                                                                                                 
GAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTG  >  NC_000913/1931996‑1932191

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: