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,8870GA99.1% 367.2 / ‑4.6 108A183V (GCC→GTC) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (61/46);  total (61/47)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 108 total aligned reads displayed.

GCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCT  >  NC_000913/2967791‑2967982
                                                                                                |                                                                                               
gcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACaa                                                                                              >  2:1590302/1‑100 (MQ=255)
gcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACaa                                                                                              >  2:2068950/1‑100 (MQ=255)
         aGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCacac                                                                                     >  2:2635121/1‑100 (MQ=255)
            cGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAg                                                                                  >  1:2102244/1‑100 (MQ=255)
            cGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAg                                                                                  >  1:2652227/1‑100 (MQ=255)
            cGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAg                                                                                  >  2:2620492/1‑100 (MQ=255)
             gTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATGCACCTGTTCCATTAAAGGTAACGTGGCAGCCTGCCAGCCTTCGACAATCGCCACACCAgg                                                                                 <  1:906920/100‑1 (MQ=255)
             gTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAgg                                                                                 <  1:2330336/100‑1 (MQ=255)
               tCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtg                                                                               <  1:380852/100‑1 (MQ=255)
               tCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtg                                                                               >  2:1747852/1‑100 (MQ=255)
                cgcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtgc                                                                              >  2:2027492/1‑100 (MQ=255)
                cgcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtgc                                                                              >  2:2021393/1‑100 (MQ=255)
                 gcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtgct                                                                             >  2:2116264/1‑100 (MQ=255)
                 gcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGtgct                                                                             >  1:2717342/1‑100 (MQ=255)
                    ttCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGcc                                                                          <  2:565079/100‑1 (MQ=255)
                     tCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCg                                                                         >  1:3018449/1‑100 (MQ=255)
                      ccAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCgg                                                                        <  1:2157189/100‑1 (MQ=255)
                       cAGAGCCGGATCGAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGt                                                                       >  1:2003022/1‑100 (MQ=255)
                       cAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGt                                                                       <  2:2508478/100‑1 (MQ=255)
                       cAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGt                                                                       <  1:326699/100‑1 (MQ=255)
                       cAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGt                                                                       >  1:2508512/1‑100 (MQ=255)
                          agCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAAt                                                                    <  1:1246164/100‑1 (MQ=255)
                          agCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAAt                                                                    <  1:237654/100‑1 (MQ=255)
                           gCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATg                                                                   >  2:1622725/1‑100 (MQ=255)
                           gCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATg                                                                   >  1:2640043/1‑100 (MQ=255)
                           gCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATg                                                                   >  1:2898843/1‑100 (MQ=255)
                           gCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATg                                                                   >  2:137456/1‑100 (MQ=255)
                           gCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATg                                                                   >  2:1527586/1‑100 (MQ=255)
                           gCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATg                                                                   <  2:3018408/100‑1 (MQ=255)
                              ggATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTATCGTGGCATCCTGCCAGCCTTCGACAATAGCCACACCAGGTGCTGCTGGTAATgtgc                                                                >  2:2523053/1‑100 (MQ=255)
                               gATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGcgcg                                                               <  1:3523/100‑1 (MQ=255)
                               gATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGcgcg                                                               <  2:2507412/100‑1 (MQ=255)
                               gATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGcgcg                                                               >  2:2669505/1‑100 (MQ=255)
                               gATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGcgcg                                                               <  1:1004836/100‑1 (MQ=255)
                                aTCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCgg                                                              <  2:976311/100‑1 (MQ=255)
                                aTCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCgg                                                              <  2:1652827/100‑1 (MQ=255)
                                  ccAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGAt                                                            <  1:2673314/100‑1 (MQ=255)
                                  ccAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGAt                                                            >  2:1173845/1‑100 (MQ=255)
                                    tgCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTc                                                          <  1:1892358/99‑1 (MQ=255)
                                     gCGTTGATGCCTGCTACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTca                                                         >  2:1892336/1‑99 (MQ=255)
                                       gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                       <  1:3001548/100‑1 (MQ=255)
                                       gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                       >  1:1366471/1‑100 (MQ=255)
                                       gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                       <  1:1381070/100‑1 (MQ=255)
                                       gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                       <  2:1410320/100‑1 (MQ=255)
                                       gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                       <  2:1071098/100‑1 (MQ=255)
                                            tGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCATACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTgg                                                  >  2:2054276/1‑100 (MQ=255)
                                             gCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGc                                                 >  1:1056569/1‑100 (MQ=255)
                                             gCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGc                                                 <  1:347328/100‑1 (MQ=255)
                                              ccTGATACACCTGTTCCATTAAATGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCg                                                <  2:664719/100‑1 (MQ=255)
                                               cTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCg                                                >  1:1002306/1‑99 (MQ=255)
                                               cTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGa                                               <  2:2620953/100‑1 (MQ=255)
                                                tGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGat                                              >  1:1413404/1‑100 (MQ=255)
                                                tGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGat                                              >  1:1541013/1‑100 (MQ=255)
                                                tGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGat                                              >  1:1412910/1‑100 (MQ=255)
                                                tGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGat                                              <  2:2652729/100‑1 (MQ=255)
                                                 gATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGata                                             >  2:1749909/1‑100 (MQ=255)
                                                     cacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                         >  2:1110795/1‑100 (MQ=255)
                                                     cacCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                         <  2:1481409/100‑1 (MQ=255)
                                                       ccTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGccc                                       >  1:159460/1‑100 (MQ=255)
                                                       ccTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGccc                                       <  2:877940/100‑1 (MQ=255)
                                                        cTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCa                                      <  2:1205138/100‑1 (MQ=255)
                                                         tGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCaa                                     >  1:3027015/1‑100 (MQ=255)
                                                          gTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCaaa                                    <  2:1905879/100‑1 (MQ=255)
                                                           ttCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAAc                                   >  2:739531/1‑100 (MQ=255)
                                                            tCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACa                                  >  1:1371394/1‑100 (MQ=255)
                                                            tCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACa                                  >  1:1864191/1‑100 (MQ=255)
                                                             ccATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACaa                                 >  2:400337/1‑100 (MQ=255)
                                                             ccATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACaa                                 >  1:807060/1‑100 (MQ=255)
                                                              cATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAg                                >  1:2018049/1‑100 (MQ=255)
                                                              cATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAg                                <  1:2123866/100‑1 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                               <  1:2500306/100‑1 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                               >  2:197435/1‑100 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                               >  1:2077576/1‑100 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                               >  2:2912462/1‑100 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                               >  2:1129770/1‑100 (MQ=255)
                                                               aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                               >  2:2626731/1‑100 (MQ=255)
                                                                ttAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGc                              <  2:1422041/100‑1 (MQ=255)
                                                                 tAAAGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCa                             >  2:2018689/1‑100 (MQ=255)
                                                                  naaGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAg                            >  2:2796499/2‑100 (MQ=255)
                                                                   aaGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGt                           >  2:790824/1‑100 (MQ=255)
                                                                    aGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTc                          >  1:130858/1‑100 (MQ=255)
                                                                    aGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTc                          <  1:2642574/100‑1 (MQ=255)
                                                                    aGGTAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTc                          >  1:8277/1‑100 (MQ=255)
                                                                       tAACGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAAc                       <  1:197442/100‑1 (MQ=255)
                                                                        aaCGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACt                      >  2:224114/1‑100 (MQ=255)
                                                                          cgcggCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACttc                    <  1:2071498/97‑3 (MQ=255)
                                                                          cGTGGCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGc                    <  1:2613630/100‑1 (MQ=255)
                                                                          cGTGGCATCCTGCAAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGc                    >  1:2073457/1‑100 (MQ=255)
                                                                             ggCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGAc                 <  2:1563879/100‑1 (MQ=255)
                                                                             ggCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGAc                 <  2:221150/100‑1 (MQ=255)
                                                                             ggCATCCTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGAc                 <  2:2102216/100‑1 (MQ=255)
                                                                                   cTGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaa           >  2:49413/1‑100 (MQ=255)
                                                                                    tGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaa          >  1:2753300/1‑100 (MQ=255)
                                                                                    tGCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaa          >  2:1491596/1‑100 (MQ=255)
                                                                                     gCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaag         <  1:2916074/100‑1 (MQ=255)
                                                                                     gCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaag         <  2:24112/100‑1 (MQ=255)
                                                                                     gCCAGCCTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaag         <  1:1561416/100‑1 (MQ=255)
                                                                                          ccTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAg    >  1:2629702/1‑100 (MQ=255)
                                                                                           cTTCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGc   >  1:2467781/1‑100 (MQ=255)
                                                                                            ttCGACAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCt  <  1:2213078/100‑1 (MQ=255)
                                                                                                |                                                                                               
GCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCT  >  NC_000913/2967791‑2967982

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: