Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,967,887 G→A 25.1% A183V (GCC→GTC)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,967,8870GA25.1% 305.8 / 109.1 171A183V (GCC→GTC) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base G (58/70);  new base A (25/18);  total (83/88)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 178 total aligned reads displayed.

AGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCA  >  NC_000913/2967788‑2967986
                                                                                                   |                                                                                                   
aGCGCCCCGGTTAGGCGTTCGCGTGCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCAGTAAAGGTAACGTGGCATCCTGCCAGCCTACgg                                                                                                     <  1:969250/100‑1 (MQ=255)
aGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCgg                                                                                                     <  2:1420076/100‑1 (MQ=255)
 gcgcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGc                                                                                                    <  1:455400/100‑1 (MQ=255)
 gcgcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGAGGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGc                                                                                                    <  1:2221044/100‑1 (MQ=255)
  cgcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCa                                                                                                   >  1:4174203/1‑100 (MQ=255)
  cgcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCa                                                                                                   >  2:586771/1‑100 (MQ=255)
     cccGGTTAGTCGTTCGCGTTCCAGGGCCGGATACAGCGTTGATGCCTGATACACCTGTTCCATTGAAGGTAACGTGGCATCCTGCCAGCCTTCGGc                                                                                                    <  1:995745/96‑1 (MQ=255)
     cccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATc                                                                                                <  2:4318341/100‑1 (MQ=255)
      ccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCg                                                                                               <  2:1247753/100‑1 (MQ=255)
       cGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGc                                                                                              <  1:1391971/100‑1 (MQ=255)
        ggTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGcc                                                                                             <  1:3416220/100‑1 (MQ=255)
            aGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCacac                                                                                         >  2:1070198/1‑100 (MQ=255)
              tCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCa                                                                                       <  1:1644374/100‑1 (MQ=255)
               cGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCCTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAACCGCCACACCAg                                                                                      >  2:273196/1‑100 (MQ=255)
               cGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCa                                                                                       >  1:1989792/1‑99 (MQ=255)
               cGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAg                                                                                      >  1:3465007/1‑100 (MQ=255)
                    gcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtgct                                                                                 >  2:799489/1‑100 (MQ=255)
                      gTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGtgctgc                                                                               >  1:1017005/1‑100 (MQ=255)
                      gTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCAACACAAGGtgctgc                                                                               >  1:353190/1‑100 (MQ=255)
                       ttCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGcc                                                                              <  2:537421/100‑1 (MQ=255)
                         ccAGAGCCGGATCCAGCGTTGGTGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCgg                                                                            <  1:3891228/100‑1 (MQ=255)
                         ccAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCgg                                                                            >  1:1965168/1‑100 (MQ=255)
                         ccAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCgg                                                                            >  1:1766002/1‑100 (MQ=255)
                          cAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGCAACGTGGCATCCTGCCAGCCTTCGAACATCGCCACACCAGGtgctcc                                                                               >  1:1721967/1‑94 (MQ=255)
                           agagCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTa                                                                          >  2:3878391/1‑100 (MQ=255)
                             agCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAAt                                                                        >  1:908349/1‑100 (MQ=255)
                                cGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATgcg                                                                     <  1:179361/100‑1 (MQ=255)
                                 ggATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATgcg                                                                     <  1:886804/99‑1 (MQ=255)
                                  gATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGcgcg                                                                   <  2:3213296/100‑1 (MQ=255)
                                  gATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGcgcg                                                                   <  2:2306617/100‑1 (MQ=255)
                                  gATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGcgcg                                                                   <  2:1634658/100‑1 (MQ=255)
                                      cAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGAtt                                                               >  1:2395433/1‑100 (MQ=255)
                                      cAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGAtt                                                               >  1:1884187/1‑100 (MQ=255)
                                        gCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGGGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcg                                                             <  2:2097059/100‑1 (MQ=255)
                                         cGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgc                                                            >  1:4275812/1‑100 (MQ=255)
                                          gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                           <  2:3988562/100‑1 (MQ=255)
                                          gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                           <  1:1070196/100‑1 (MQ=255)
                                           ttGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGCATTcgc                                                            >  1:1280714/1‑98 (MQ=255)
                                              aTGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTg                                                       <  1:3330025/100‑1 (MQ=255)
                                              aTGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTg                                                       >  2:1782615/1‑100 (MQ=255)
                                              aTGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTg                                                       >  2:2043758/1‑100 (MQ=255)
                                               tGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTgg                                                      >  1:4005766/1‑100 (MQ=255)
                                               tGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTgg                                                      >  1:4281904/1‑100 (MQ=255)
                                                gCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGc                                                     <  2:1199628/100‑1 (MQ=255)
                                                gCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGc                                                     <  2:2316002/100‑1 (MQ=255)
                                                gCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGc                                                     >  2:2486264/1‑100 (MQ=255)
                                                   tGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGa                                                   <  1:2126553/99‑1 (MQ=255)
                                                   tGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGat                                                  >  1:3694176/1‑100 (MQ=255)
                                                      tACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATAtt                                               >  1:3749230/1‑100 (MQ=255)
                                                        cacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                             >  2:3251182/1‑100 (MQ=255)
                                                        cacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                             >  2:3553246/1‑100 (MQ=255)
                                                        cacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                             <  2:2246518/100‑1 (MQ=255)
                                                         acCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGcc                                            <  2:1138741/100‑1 (MQ=255)
                                                          tcttccgatctAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGccc                                           <  2:89239/90‑1 (MQ=255)
                                                          ccTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGccc                                           >  2:408710/1‑100 (MQ=255)
                                                          ccTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGccc                                           >  2:2379471/1‑100 (MQ=255)
                                                              ttCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAAc                                       <  1:2661745/100‑1 (MQ=255)
                                                                 cATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAg                                    >  1:2377513/1‑100 (MQ=255)
                                                                 cATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAg                                    <  1:2528029/100‑1 (MQ=255)
                                                                  aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                   >  1:697921/1‑100 (MQ=255)
                                                                  aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                   <  1:2226214/100‑1 (MQ=255)
                                                                  aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                   >  1:2257959/1‑100 (MQ=255)
                                                                   ttAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGc                                  >  2:1842106/1‑100 (MQ=255)
                                                                   ttAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGc                                  <  2:2986205/100‑1 (MQ=255)
                                                                    tAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCa                                 <  2:1332503/100‑1 (MQ=255)
                                                                     aaaGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAg                                >  2:1063867/1‑100 (MQ=255)
                                                                     aaaGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGATGCTGCCGGTAATGCGCGGAGTCGCGTCTGGCGATATTGCCCAAACAAGGCa                                 >  2:4165462/1‑99 (MQ=255)
                                                                     aaaGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCagatcggaaga                                >  1:89239/1‑90 (MQ=255)
                                                                       aGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTc                              <  1:1978708/100‑1 (MQ=255)
                                                                       aGGTAACGTGACATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTc                              >  2:1843775/1‑100 (MQ=255)
                                                                        ggTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGGTTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCa                             >  1:2638912/1‑100 (MQ=255)
                                                                        ggTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGCAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCa                             <  1:2439295/100‑1 (MQ=255)
                                                                          tAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAAc                           >  1:1676021/1‑100 (MQ=255)
                                                                           aaCGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACt                          >  2:646932/1‑100 (MQ=255)
                                                                           aaCGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACt                          <  2:3437703/100‑1 (MQ=255)
                                                                           aaCGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACt                          >  2:791506/1‑100 (MQ=255)
                                                                           aaCGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACt                          >  2:2432649/1‑100 (MQ=255)
                                                                            aCGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTg                         <  2:4206511/100‑1 (MQ=255)
                                                                            aCGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTg                         <  2:1435/100‑1 (MQ=255)
                                                                             gctggcaggatGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGc                        >  2:331039/11‑100 (MQ=255)
                                                                             cGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTg                         <  2:4274266/99‑1 (MQ=255)
                                                                             cGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGc                        <  2:4408737/100‑1 (MQ=255)
                                                                              gTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCg                       >  2:497307/1‑100 (MQ=255)
                                                                                ggCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGAc                     <  2:1592293/100‑1 (MQ=255)
                                                                                 gCAGCCTGCCAGCCTTCGGTAATCGCACTACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTGTGGCGATATTGCCCAAACAAGGCAGTCa                             <  1:2883872/91‑1 (MQ=255)
                                                                                  cATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTg                   <  2:2615116/100‑1 (MQ=255)
                                                                                   aTCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTga                  >  2:2384709/1‑100 (MQ=255)
                                                                                       tGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGGCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaa              >  1:1141957/1‑100 (MQ=255)
                                                                                        gCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaa              <  1:131568/99‑1 (MQ=255)
                                                                                         ccAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaaga            <  1:3786626/100‑1 (MQ=255)
                                                                                         ccAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaaga            <  2:787844/100‑1 (MQ=255)
                                                                                         ccAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaaga            >  2:945662/1‑100 (MQ=255)
                                                                                          cAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaagaa           <  1:2117908/100‑1 (MQ=255)
                                                                                           aGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAAt          >  2:1315211/1‑100 (MQ=255)
                                                                                             ccTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAg        >  2:3586984/1‑100 (MQ=255)
                                                                                              cTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGc       >  1:1444962/1‑100 (MQ=255)
                                                                                               ttCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCt      <  1:3940679/100‑1 (MQ=255)
                                                                                                tCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTg     >  2:1561878/1‑100 (MQ=255)
                                                                                                  ggCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAATAAGGCAGTCAACTTCGACTGAGAAAGAATAGCTg     >  1:626260/1‑98 (MQ=255)
                                                                                                  ggCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGcc   >  1:3025518/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCA  >  NC_000913/2967788‑2967986

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