Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 436,702 G→A 5.8% G38G (GGG→GGA pgpA → phosphatidylglycerophosphatase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,7020GA 5.8% 326.1 / 6.4 104G38G (GGG→GGApgpAphosphatidylglycerophosphatase A
Reads supporting (aligned to +/- strand):  ref base G (54/42);  new base A (2/4);  total (56/48)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 111 total aligned reads displayed.

GTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATGG  >  NC_000913/436553‑436851
                                                                                                                                                     |                                                                                                                                                      
gTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTggg                                                                                                                                                        >  1:758710/1‑150 (MQ=255)
gTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTggg                                                                                                                                                        <  1:1402454/150‑1 (MQ=255)
 tGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCATCGTCCCAGGAAc                                                                                                                                                      <  2:1699235/150‑1 (MQ=255)
      gCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATg                                                                                                                                                 >  1:1937740/1‑150 (MQ=255)
       cGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATgg                                                                                                                                                >  2:1259093/1‑150 (MQ=255)
       cGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATgg                                                                                                                                                <  1:728126/150‑1 (MQ=255)
       cGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCATCGTCCCAGGAACGATcgg                                                                                                                                               >  2:219327/1‑150 (MQ=255)
          aCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCATCGTCCCAGGAACGATCGGGCt                                                                                                                                             <  1:219327/150‑1 (MQ=255)
          aCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCATCGTCCCAGGAACGATCGGGCt                                                                                                                                             <  1:1211186/150‑1 (MQ=255)
            cTGTTACATTAGACTGGAAAGTATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCTTCGCATATACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCTATCGTTCCTGGTACGATGTTCTCg                                                                                                                                           <  2:1178902/150‑1 (MQ=255)
            cTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCATCGTCCCAGGAACGATCGGGCtta                                                                                                                                           <  2:1489878/150‑3 (MQ=255)
             tGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGc                                                                                                                                          >  1:1725241/1‑150 (MQ=255)
              gTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCt                                                                                                                                         >  2:1597434/1‑150 (MQ=255)
                    ttAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGc                                                                                                                                   >  1:1406311/1‑150 (MQ=255)
                    ttAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGc                                                                                                                                   <  1:689244/150‑1 (MQ=255)
                    ttAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGc                                                                                                                                   >  1:332343/1‑150 (MQ=255)
                     tAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCg                                                                                                                                  <  1:1878642/150‑1 (MQ=255)
                      aGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGa                                                                                                                                 <  1:1597434/150‑1 (MQ=255)
                       gACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGAt                                                                                                                                >  2:877863/1‑150 (MQ=255)
                          tGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCACGCTGGCAGCGATTcc                                                                                                                             <  1:1473596/150‑1 (MQ=255)
                           ggAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCg                                                                                                                            >  2:1995569/1‑150 (MQ=255)
                               aGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCt                                                                                                                        <  2:1033136/150‑1 (MQ=255)
                                ggATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTg                                                                                                                       <  1:1069355/150‑1 (MQ=255)
                                ggATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTg                                                                                                                       >  2:1198119/1‑150 (MQ=255)
                                ggATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTg                                                                                                                       >  2:1922765/1‑150 (MQ=255)
                                     tGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGCTCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATc                                                                                                                  >  1:328734/1‑150 (MQ=255)
                                       aCCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtg                                                                                                                <  2:860795/150‑1 (MQ=255)
                                       aCCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtg                                                                                                                >  1:1796146/1‑150 (MQ=255)
                                        ccATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtga                                                                                                               >  1:148143/1‑150 (MQ=255)
                                         cATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgat                                                                                                              <  1:1259093/150‑1 (MQ=255)
                                           ttttGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgatga                                                                                                            <  2:1088771/150‑1 (MQ=255)
                                            tttGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGAc                                                                                                           <  1:1237640/150‑1 (MQ=255)
                                              tGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCt                                                                                                         >  2:848024/1‑150 (MQ=255)
                                                ccACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCttt                                                                                                       >  1:6127/1‑150 (MQ=255)
                                                  aCGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCttttt                                                                                                     >  1:1193965/1‑150 (MQ=255)
                                                    gCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGc                                                                                                   <  1:1702550/150‑1 (MQ=255)
                                                     ccATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCCTCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTGGTGGTAGCTGATGCCCTTTTTGcc                                                                                                  >  2:489374/1‑150 (MQ=255)
                                                         aaaGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTgg                                                                                              >  2:11501/1‑150 (MQ=255)
                                                          aaGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGc                                                                                             <  1:795197/150‑1 (MQ=255)
                                                             aTGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGc                                                                                          <  2:681520/150‑1 (MQ=255)
                                                             aTGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGc                                                                                          <  2:1796146/150‑1 (MQ=255)
                                                                tCGCGAAAAGTCGCCTGAATATTAGTAATCCGTTGCATCTACTTGCTGTCGGTTTCGGAATTGGATTAAGCCCGATCTTTCCTTGTACTATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGctct                                                                                       <  2:1849749/150‑1 (MQ=255)
                                                                tCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGctct                                                                                       <  2:257180/150‑1 (MQ=255)
                                                                   cgAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACt                                                                                    <  1:895288/150‑1 (MQ=255)
                                                                         gTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtgg                                                                              >  1:758977/1‑150 (MQ=255)
                                                                           cGCCTGAAGATGAGTAATCGGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtggtg                                                                            <  2:1685667/150‑1 (MQ=255)
                                                                           cGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtggtg                                                                            >  1:1661896/1‑150 (MQ=255)
                                                                           cGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTATGGTATCTGATGACCTTTTTGCCATGGCAGCTATACTCGCtggtg                                                                            >  1:1493860/1‑150 (MQ=255)
                                                                            gCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtggtgg                                                                           <  1:504824/150‑1 (MQ=255)
                                                                             ccTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTggtggt                                                                          >  1:144139/1‑150 (MQ=255)
                                                                              cTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGgtggtg                                                                         <  2:1718158/150‑1 (MQ=255)
                                                                                gAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGAt                                                                       >  2:365680/1‑150 (MQ=255)
                                                                                   gaagAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCGGGTGGTGATGcg                                                                    >  2:792235/4‑149 (MQ=255)
                                                                                   gATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGCCCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCt                                                                    >  2:23679/1‑150 (MQ=255)
                                                                                   gATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCt                                                                    >  1:336061/1‑150 (MQ=255)
                                                                                       aGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTgggg                                                                >  1:1730476/1‑150 (MQ=255)
                                                                                       aGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTgggg                                                                <  1:111266/150‑1 (MQ=255)
                                                                                        gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGa                                                               >  2:1375517/1‑150 (MQ=255)
                                                                                              cGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTa                                                         <  2:144139/150‑1 (MQ=255)
                                                                                                tGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTg                                                                                                                                        <  2:928033/69‑1 (MQ=255)
                                                                                                tGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTg                                                                                                                                        >  1:928033/1‑69 (MQ=255)
                                                                                                 ggCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCg                                                      >  2:505810/1‑150 (MQ=255)
                                                                                                             cTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCttt                                          >  2:1777387/1‑150 (MQ=255)
                                                                                                             cTGTCGGATTCGCAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCCATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGCGGCTCTGTATCGGCGTCTATCttt                                          >  2:830049/1‑150 (MQ=255)
                                                                                                              tGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTg                                         <  1:505810/150‑1 (MQ=255)
                                                                                                              tGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTg                                         <  2:1725241/150‑1 (MQ=255)
                                                                                                              tGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGAAGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTg                                         <  2:1100848/150‑1 (MQ=255)
                                                                                                               gTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGt                                        <  2:1730476/150‑1 (MQ=255)
                                                                                                               gTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGt                                        >  2:43800/1‑150 (MQ=255)
                                                                                                               gTCGGATACGGAAGTGGATTAAGCCCGATCGTTCCTGCGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGt                                        >  2:751186/1‑150 (MQ=255)
                                                                                                                 cGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGtca                                      <  1:364293/150‑1 (MQ=255)
                                                                                                                 cGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGtca                                      >  2:225571/1‑150 (MQ=255)
                                                                                                                  ggATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTTTGGTATCTGATGACCTTTTTTCCGTGTCATCTCTAGTCGCTGTTGGTTATGCTGGTGATCTGTATCGGGGTCTATATTTGtcat                                     <  2:332343/150‑1 (MQ=255)
                                                                                                                    aTTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGtcatca                                   >  1:1694427/1‑150 (MQ=255)
                                                                                                                    aTTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGtcatca                                   >  1:1687505/1‑150 (MQ=255)
                                                                                                                       cGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAAc                                >  1:1175808/1‑150 (MQ=255)
                                                                                                                       cGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATAGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAAc                                >  1:1937888/1‑150 (MQ=255)
                                                                                                                        ggAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACg                               >  2:820248/1‑150 (MQ=255)
                                                                                                                            gTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGa                           >  2:1675705/1‑150 (MQ=255)
                                                                                                                             tGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGaa                          <  2:328734/150‑1 (MQ=255)
                                                                                                                              ggATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGaaa                         >  2:1564241/1‑150 (MQ=255)
                                                                                                                               gATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAg                        <  2:672029/150‑1 (MQ=255)
                                                                                                                                   aaGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAg                        <  1:446092/146‑1 (MQ=255)
                                                                                                                                   aaGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAg                        >  2:446092/1‑146 (MQ=255)
                                                                                                                                    aGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATg                   <  1:1198119/150‑1 (MQ=255)
                                                                                                                                     gCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATgg                  >  2:913666/1‑150 (MQ=255)
                                                                                                                                     agcccaTCGTCCCAGGAACGATCGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATg                   >  2:1202206/6‑150 (MQ=255)
                                                                                                                                      cccGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATggg                 <  1:751186/150‑1 (MQ=255)
                                                                                                                                      cccGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATggg                 <  2:1193965/150‑1 (MQ=255)
                                                                                                                                        cGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGgtg               >  2:344986/1‑150 (MQ=255)
                                                                                                                                          aTCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGtgtg             <  1:1777387/150‑1 (MQ=255)
                                                                                                                                             gTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCAc          >  2:33642/1‑150 (MQ=255)
                                                                                                                                             gTTCCTGGGACGATAGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCAc          <  2:1937888/150‑1 (MQ=255)
                                                                                                                                                ccTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGAt       >  1:1901844/1‑150 (MQ=255)
                                                                                                                                                  tGGGTCTATGGGCTCGCTGGCAGCGTTTCATTTCTGGTATCTGATTACCTTTTTGCCCTGGCAGGTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCa     <  2:1937740/150‑1 (MQ=255)
                                                                                                                                                   gggACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTTTCTGATGACCTTTTTGCCCTGGCAGTTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCg                                                   >  2:595291/1‑103 (MQ=255)
                                                                                                                                                   gggACGATGGGCTCGCTGGCAGAGATTCAGTTCTGGTATCTGATGACCTTTTTGCACTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCAt    >  2:914669/1‑150 (MQ=255)
                                                                                                                                                   gggACGAGGGGCTCGCTGGCAGCGATTCCGTTCTGGTTTCTGATGACCTTTTTGCCCTGGCAGTTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCg                                                   <  1:595291/103‑1 (MQ=255)
                                                                                                                                                    ggACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCGGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATg   <  2:1205425/150‑1 (MQ=255)
                                                                                                                                                     gACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATgg  <  2:1406311/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGATGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATGG  >  NC_000913/436553‑436851

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