Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,967,887 G→A 100% 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,8870GA88.7% 257.5 / 16.5 95A183V (GCC→GTC) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base G (5/6);  new base A (48/36);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTG  >  NC_000913/2967789‑2967983
                                                                                                  |                                                                                                
gcgcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTTGCATCCTGCCAGCCTTCGGc                                                                                                 <  2:213068/100‑1 (MQ=255)
   ccccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAAt                                                                                              >  2:721382/1‑100 (MQ=255)
     ccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCg                                                                                            >  1:1596417/1‑100 (MQ=255)
       ggTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGcc                                                                                          <  1:2576530/100‑1 (MQ=255)
            gTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCCGCCAGCCTTCGACAATCGCCACAcc                                                                                     >  2:2017382/1‑100 (MQ=255)
              cGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAACGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAg                                                                                   >  2:2226642/1‑100 (MQ=255)
                 tCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtg                                                                                <  1:1356530/100‑1 (MQ=255)
                     gTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtgctgc                                                                            <  1:625378/100‑1 (MQ=255)
                     gTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtgctgc                                                                            <  1:684566/100‑1 (MQ=255)
                      ttCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGcc                                                                           >  1:676991/1‑100 (MQ=255)
                      ttCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGcc                                                                           <  2:150820/100‑1 (MQ=255)
                      ttCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGcc                                                                           <  1:856463/100‑1 (MQ=255)
                       tCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCg                                                                          <  2:349394/100‑1 (MQ=255)
                        ccAGAGCCGGATCCCGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGATGCTGcc                                                                           >  2:178252/1‑98 (MQ=255)
                        ccAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCgg                                                                         <  1:458615/100‑1 (MQ=255)
                        ccAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCgg                                                                         <  2:2092547/100‑1 (MQ=255)
                         cAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGt                                                                        <  1:695201/100‑1 (MQ=255)
                             gCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAAt                                                                     <  2:704839/99‑1 (MQ=255)
                                ggATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATgcgc                                                                 <  1:745122/100‑1 (MQ=255)
                                ggATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATgcgc                                                                 >  1:2364130/1‑100 (MQ=255)
                                ggATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATgcgc                                                                 >  1:272712/1‑100 (MQ=255)
                                ggATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATgcgc                                                                 >  1:498214/1‑100 (MQ=255)
                                  aTCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCgg                                                               >  2:905361/1‑100 (MQ=255)
                                  aTCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCgg                                                               >  2:576161/1‑100 (MQ=255)
                                       gCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcg                                                          >  2:63229/1‑100 (MQ=255)
                                         gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                        <  1:31360/100‑1 (MQ=255)
                                         gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                        <  2:1306369/100‑1 (MQ=255)
                                         gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                        <  1:2143837/100‑1 (MQ=255)
                                         gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                        <  2:1718199/100‑1 (MQ=255)
                                          ttGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGt                                                       >  1:2603461/1‑100 (MQ=255)
                                              tGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTgg                                                   <  1:1389574/100‑1 (MQ=255)
                                              tGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTgg                                                   >  2:32515/1‑100 (MQ=255)
                                                ccTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCg                                                 >  1:459567/1‑100 (MQ=255)
                                                  tGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGat                                               <  2:272709/100‑1 (MQ=255)
                                                   gATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGata                                              >  2:1973265/1‑100 (MQ=255)
                                                    aTACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGata                                              <  2:2596028/99‑1 (MQ=255)
                                                     tACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATAtt                                            >  2:1716701/1‑100 (MQ=255)
                                                     tACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATAtt                                            >  1:189470/1‑100 (MQ=255)
                                                     tACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATAtt                                            >  1:2350406/1‑100 (MQ=255)
                                                      acacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTg                                           >  2:259286/1‑100 (MQ=255)
                                                        acCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGcc                                         <  2:1350822/100‑1 (MQ=255)
                                                         ccTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGccc                                        <  2:551281/100‑1 (MQ=255)
                                                         ccTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGccc                                        >  2:1767202/1‑100 (MQ=255)
                                                            gTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCaaa                                     >  2:646539/1‑100 (MQ=255)
                                                            gTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCaaa                                     >  2:1870048/1‑100 (MQ=255)
                                                            gTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCaaa                                     >  1:2580235/1‑100 (MQ=255)
                                                            gTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCaaa                                     >  2:1911687/1‑100 (MQ=255)
                                                              tCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACa                                   >  2:1924611/1‑100 (MQ=255)
                                                               ccATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACaa                                  >  1:74761/1‑100 (MQ=255)
                                                                cATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAg                                 >  1:1095618/1‑100 (MQ=255)
                                                                 aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCGAAACAAgg                                >  2:1491788/1‑100 (MQ=255)
                                                                 aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                >  1:2425718/1‑100 (MQ=255)
                                                                 aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                >  1:1394538/1‑100 (MQ=255)
                                                                 aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                <  2:89420/100‑1 (MQ=255)
                                                                 aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                <  1:2273376/100‑1 (MQ=255)
                                                                  ttAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGc                               <  1:1945961/100‑1 (MQ=255)
                                                                    aaaGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAg                             <  2:1640775/100‑1 (MQ=255)
                                                                     aaGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGt                            >  1:2480494/1‑100 (MQ=255)
                                                                     aaGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGt                            <  1:1433313/100‑1 (MQ=255)
                                                                      aGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTc                           >  1:585949/1‑100 (MQ=255)
                                                                       ggTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTc                           <  2:1246799/99‑1 (MQ=255)
                                                                       ggTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCa                          >  2:725007/1‑100 (MQ=255)
                                                                       ggTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCa                          <  1:2121731/100‑1 (MQ=255)
                                                                       ggTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCa                          <  2:2506612/100‑1 (MQ=255)
                                                                        gTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCaa                         >  1:1135330/1‑100 (MQ=255)
                                                                         tAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAAc                        <  2:2350393/100‑1 (MQ=255)
                                                                             gTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCg                    >  1:109147/1‑100 (MQ=255)
                                                                             gTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCg                    >  1:378285/1‑100 (MQ=255)
                                                                              tGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGa                   >  2:1973163/1‑100 (MQ=255)
                                                                              tGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGAGATATTGCCCAAACAAGGCAGTCAACTGCg                    >  1:2051711/1‑99 (MQ=255)
                                                                               ggCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGAc                  <  2:2006715/100‑1 (MQ=255)
                                                                               ggCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGAc                  >  1:2167205/1‑100 (MQ=255)
                                                                               ggCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGAc                  <  2:2425705/100‑1 (MQ=255)
                                                                                gCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACt                 >  2:136975/1‑100 (MQ=255)
                                                                                 cATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTg                >  2:1134243/1‑100 (MQ=255)
                                                                                  aTCCTGCCAGCCTTCGACAATCGCCGCACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTg                >  1:1941545/1‑99 (MQ=255)
                                                                                   tCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTgag              <  1:1324950/100‑1 (MQ=255)
                                                                                    ccTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTgaga             >  2:2187560/1‑100 (MQ=255)
                                                                                     cTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaa            <  1:955206/100‑1 (MQ=255)
                                                                                      tGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGCCTGagaaa           >  2:1766950/1‑100 (MQ=255)
                                                                                      tGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaa           <  1:1886576/100‑1 (MQ=255)
                                                                                       tccAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaag          <  2:294642/99‑1 (MQ=255)
                                                                                        ccAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaaga         >  1:294645/1‑100 (MQ=255)
                                                                                          aGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAAt       <  2:2167193/100‑1 (MQ=255)
                                                                                          aGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAAt       <  2:1847694/100‑1 (MQ=255)
                                                                                           gCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATa      >  2:1402453/1‑100 (MQ=255)
                                                                                           gCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATa      >  2:1254484/1‑100 (MQ=255)
                                                                                           gCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATAGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATa      <  2:1600777/100‑1 (MQ=255)
                                                                                            ccTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAg     >  1:1056015/1‑100 (MQ=255)
                                                                                             cTTCGGCAATCGCCACACCTGGTGCTGCCGGTAATGCGCGGAGTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGc    <  2:2076280/100‑1 (MQ=255)
                                                                                             cTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGc    >  1:2581464/1‑100 (MQ=255)
                                                                                             cTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGc    <  2:2086572/100‑1 (MQ=255)
                                                                                             cTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGc    >  1:1529870/1‑100 (MQ=255)
                                                                                             cTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGACAAAGAATAGc    >  2:1717402/1‑100 (MQ=255)
                                                                                               tCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTg  >  1:874279/1‑100 (MQ=255)
                                                                                                  |                                                                                                
GCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTG  >  NC_000913/2967789‑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: