Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,967,887 G→A 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,8870GA96.8% 311.6 / ‑1.1 93A183V (GCC→GTC) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (50/40);  total (52/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTG  >  NC_000913/2967791‑2967983
                                                                                                |                                                                                                
gcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACaa                                                                                               >  1:388455/1‑100 (MQ=255)
 ccccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTCAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAAt                                                                                              >  2:2371515/1‑100 (MQ=255)
 ccccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAAt                                                                                              >  2:2101243/1‑100 (MQ=255)
   ccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCg                                                                                            <  1:2511197/100‑1 (MQ=255)
   ccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCg                                                                                            >  2:1231751/1‑100 (MQ=255)
    cGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGc                                                                                           <  2:2130176/100‑1 (MQ=255)
     ggTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGAGAATCGcc                                                                                          <  1:1453267/100‑1 (MQ=255)
       tCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCcac                                                                                        <  1:2544499/100‑1 (MQ=255)
          gTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACAcc                                                                                     <  1:1569751/100‑1 (MQ=255)
               tCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtg                                                                                >  2:2682298/1‑100 (MQ=255)
                cgcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtgc                                                                               <  1:2317588/100‑1 (MQ=255)
                     tCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCg                                                                          >  2:2188388/1‑100 (MQ=255)
                         gagCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTaa                                                                      <  1:655617/100‑1 (MQ=255)
                          agCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAAt                                                                     <  2:1965332/100‑1 (MQ=255)
                           gCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAAt                                                                     >  2:257480/1‑99 (MQ=255)
                           gCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATg                                                                    >  2:605495/1‑100 (MQ=255)
                             cGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATgcg                                                                  >  1:608319/1‑100 (MQ=255)
                              ggATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATgcgc                                                                 >  2:133180/1‑100 (MQ=255)
                              ggATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATgcgc                                                                 >  2:2317038/1‑100 (MQ=255)
                              ggATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATgcgc                                                                 >  2:1843827/1‑100 (MQ=255)
                                aTCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCgg                                                               >  2:196215/1‑100 (MQ=255)
                                   cAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGAtt                                                            >  1:1369926/1‑100 (MQ=255)
                                   cAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGAtt                                                            >  2:2149853/1‑100 (MQ=255)
                                    aGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCGACACCAGGAGCTGCCGGTAATGCGCGGATTc                                                           >  2:1861717/1‑100 (MQ=255)
                                    aGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTc                                                           >  2:661301/1‑100 (MQ=255)
                                       gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                        <  2:2265831/100‑1 (MQ=255)
                                          gATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCt                                                     >  2:1880903/1‑100 (MQ=255)
                                            tGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTgg                                                   >  1:930494/1‑100 (MQ=255)
                                            tGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTgg                                                   >  1:881256/1‑100 (MQ=255)
                                            tGCCTGATACAACTGGTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGCCATGCGCGGCTTCCAGTCtcg                                                   >  2:1178483/1‑98 (MQ=255)
                                             gCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGc                                                  <  1:528205/100‑1 (MQ=255)
                                              gcTGATACACCTGTTCCATTAAAGGAAACGTGGCATACTGCCAGCCTTCGACAATCGCCACACCACGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGc                                                  >  1:669409/2‑99 (MQ=255)
                                              ccTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCg                                                 <  1:694453/100‑1 (MQ=255)
                                              ccTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCg                                                 <  1:1483969/100‑1 (MQ=255)
                                               cTGATACTCCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGa                                                <  1:1366713/100‑1 (MQ=255)
                                               cTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGa                                                <  1:593632/100‑1 (MQ=255)
                                               cTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGa                                                >  2:542488/1‑100 (MQ=255)
                                                 gATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGata                                              >  1:477097/1‑100 (MQ=255)
                                                 gATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGata                                              >  2:1642796/1‑100 (MQ=255)
                                                    acacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATAtt                                            >  1:2011738/1‑99 (MQ=255)
                                                     cacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                          >  1:374038/1‑100 (MQ=255)
                                                     cacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCCCACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATAtt                                            >  1:99686/1‑98 (MQ=255)
                                                     cacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                          >  2:1241281/1‑100 (MQ=255)
                                                     cacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                          >  2:1241234/1‑100 (MQ=255)
                                                     cacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                          >  1:566300/1‑100 (MQ=255)
                                                     cacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                          >  1:2210160/1‑100 (MQ=255)
                                                       ccTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGccc                                        >  2:2245646/1‑100 (MQ=255)
                                                        cTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCa                                       <  2:1369917/100‑1 (MQ=255)
                                                        cTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCa                                       <  1:1970765/100‑1 (MQ=255)
                                                         tGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCaa                                      >  2:230956/1‑100 (MQ=255)
                                                          gTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCaaa                                     >  2:1163010/1‑100 (MQ=255)
                                                          gTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCaaa                                     >  2:2171124/1‑100 (MQ=255)
                                                            tCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACa                                   >  2:2072930/1‑100 (MQ=255)
                                                            tCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACa                                   >  1:2125638/1‑100 (MQ=255)
                                                             ccATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACaa                                  <  2:377020/100‑1 (MQ=255)
                                                             ccATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACaa                                  >  2:2059828/1‑100 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTTCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                <  2:1347203/100‑1 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                <  1:139255/100‑1 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                <  2:1570419/100‑1 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                <  2:1625736/100‑1 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                <  1:2054142/100‑1 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                <  2:2000833/100‑1 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                >  1:2270887/1‑100 (MQ=255)
                                                                ttAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGc                               <  1:2630559/100‑1 (MQ=255)
                                                                 tAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCa                              <  1:378959/100‑1 (MQ=255)
                                                                  aaaGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAg                             <  1:1179420/100‑1 (MQ=255)
                                                                  aaaGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAg                             <  1:1231758/100‑1 (MQ=255)
                                                                  aaaGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAg                             <  1:2049839/100‑1 (MQ=255)
                                                                  aaaGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAg                             >  1:86653/1‑100 (MQ=255)
                                                                   aaGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGCCTGGCGATATTGCCCAAACAAGGCAGt                            <  1:2374510/100‑1 (MQ=255)
                                                                    aGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTc                           >  1:655458/1‑100 (MQ=255)
                                                                      gTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCaa                         >  1:1251191/1‑100 (MQ=255)
                                                                         aCGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTg                      <  1:1107020/100‑1 (MQ=255)
                                                                         aCGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTg                      <  1:940559/100‑1 (MQ=255)
                                                                          cGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTg                      <  2:2206651/99‑1 (MQ=255)
                                                                          cGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGc                     <  2:279393/100‑1 (MQ=255)
                                                                          cGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGc                     >  2:124705/1‑100 (MQ=255)
                                                                           gTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCg                    >  2:1655720/1‑100 (MQ=255)
                                                                              gCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACt                 >  1:1279301/1‑100 (MQ=255)
                                                                               cATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTg                <  2:152018/100‑1 (MQ=255)
                                                                                 tCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTgag              <  1:1843842/100‑1 (MQ=255)
                                                                                 tCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTgag              >  1:1665029/1‑100 (MQ=255)
                                                                                  ccTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTgaga             >  2:725016/1‑100 (MQ=255)
                                                                                  ccTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTgaga             <  1:634994/100‑1 (MQ=255)
                                                                                    tGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaa           <  2:608314/100‑1 (MQ=255)
                                                                                        aGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAAt       >  1:738245/1‑100 (MQ=255)
                                                                                         gCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATa      <  2:2098707/100‑1 (MQ=255)
                                                                                         gCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATa      <  1:1656344/100‑1 (MQ=255)
                                                                                         gCCTTCGACAATCGCCACACCAGGTGCTGCCGGAAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTCAGAAAGAATa      >  2:1808018/1‑100 (MQ=255)
                                                                                          ccTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAg     <  2:1726499/100‑1 (MQ=255)
                                                                                          ccTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAg     >  2:362913/1‑100 (MQ=255)
                                                                                           cTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGc    <  2:86651/100‑1 (MQ=255)
                                                                                             tCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTg  >  1:742058/1‑100 (MQ=255)
                                                                                                |                                                                                                
GCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTG  >  NC_000913/2967791‑2967983

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