Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 436,702 G→A 6.0% 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 6.0% 540.5 / 15.4 164G38G (GGG→GGA) pgpAphosphatidylglycerophosphatase A
Reads supporting (aligned to +/- strand):  ref base G (75/79);  new base A (5/5);  total (80/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

GCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATG  >  NC_000913/436559‑436850
                                                                                                                                               |                                                                                                                                                     
gCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGGAGTGGATTAAGCCCGATCGTTCCTGGGACGAt                                                                                                                                                  <  1:684091/149‑1 (MQ=255)
gCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAt                                                                                                                                                  <  4:1099418/149‑1 (MQ=255)
             aTTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCa                                                                                                                                    >  2:1008744/1‑149 (MQ=255)
             aTTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCa                                                                                                                                    >  4:294934/1‑149 (MQ=255)
             aTTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGAACGAt                                                                                                                                                  >  1:78192/1‑136 (MQ=255)
             aTTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGAACGAt                                                                                                                                                  <  2:78192/136‑1 (MQ=255)
              ttAGACTGGAAAGGATATGACCATTTTGCCCCGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGc                                                                                                                                             <  2:979649/139‑1 (MQ=255)
              ttAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAG‑AATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGc                                                                                                                                  >  4:78799/1‑149 (MQ=255)
                      gAAAGGTTATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑gg                                                                                                                                               >  2:1039629/1‑129 (MQ=255)
                      gAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAt                                                                                                                                                  >  2:249101/1‑127 (MQ=255)
                      gAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAt                                                                                                                                                  <  1:249101/127‑1 (MQ=255)
                      gAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCC‑ATCGTCCCAGGAACGAT‑cgggc                                                                                                                                            <  3:17838/131‑3 (MQ=255)
                      gAAAGGATATGACCATTTCGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCg                                                                                                                           >  3:511394/1‑149 (MQ=255)
                        aaGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑ggg                                                                                                                                              <  1:529578/128‑1 (MQ=255)
                            atatGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTg                                                                                                                      >  3:765734/1‑148 (MQ=255)
                            atatGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTg                                                                                                                      <  4:765734/148‑1 (MQ=255)
                              atGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTgg                                                                                                                                      >  2:829840/1‑130 (MQ=255)
                                gACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTg                                                                                                                      >  2:514407/1‑144 (MQ=255)
                                   cATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtga                                                                                                              <  2:1349624/149‑1 (MQ=255)
                                   cATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtga                                                                                                              >  4:1453672/1‑149 (MQ=255)
                                   cATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtga                                                                                                              <  1:826405/149‑1 (MQ=255)
                                    aTTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgat                                                                                                             <  1:1221798/149‑1 (MQ=255)
                                    aTTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgat                                                                                                             <  4:1322922/149‑1 (MQ=255)
                                    aTTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgat                                                                                                             <  3:1453672/149‑1 (MQ=255)
                                    aTTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgat                                                                                                             >  1:540338/1‑149 (MQ=255)
                                           cACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCttt                                                                                                      >  3:571073/1‑149 (MQ=255)
                                            aCGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCtttt                                                                                                     <  3:798551/149‑1 (MQ=255)
                                            aCGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCC‑ATCGTCCCAGGAACGAT‑cgggc                                                                                                                                            >  1:499286/1‑107 (MQ=38)
                                            aCGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCC‑ATCGTCCCAGGAACGAT‑cgggc                                                                                                                                            <  2:499286/109‑3 (MQ=38)
                                             cGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCttttt                                                                                                    <  2:540338/149‑1 (MQ=255)
                                                cATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGtt                                                                                                                         <  3:1404989/125‑1 (MQ=255)
                                                 aTAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGccc                                                                                                <  4:221768/149‑1 (MQ=255)
                                                  tAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCt                                                                                               <  1:1008744/149‑1 (MQ=255)
                                                    aaGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgat                                                                                                             <  3:130315/133‑1 (MQ=255)
                                                              gAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTa                                                                                                                   >  2:663596/1‑117 (MQ=255)
                                                                 aaGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGc                                                                                <  3:983429/149‑1 (MQ=255)
                                                                 aaGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGc                                                                                <  3:294934/149‑1 (MQ=255)
                                                                   gTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtg                                                                              <  1:1281807/149‑1 (MQ=255)
                                                                    tCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTg                                                                                                                                       <  3:782441/91‑1 (MQ=255)
                                                                    tCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTg                                                                                                                                       >  4:782441/1‑91 (MQ=255)
                                                                    tCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtgg                                                                             >  2:1280772/1‑149 (MQ=255)
                                                                    tCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtgg                                                                             >  4:231405/1‑149 (MQ=255)
                                                                        cTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGctc                                                                                       <  2:1392419/135‑1 (MQ=255)
                                                                        cTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGctc                                                                                       >  1:1392419/1‑135 (MQ=255)
                                                                          gAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAg                                                                                          <  2:368120/130‑1 (MQ=255)
                                                                          gAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAg                                                                                          >  1:368120/1‑130 (MQ=255)
                                                                              aTGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCt                                                                   >  3:278083/1‑149 (MQ=255)
                                                                                gAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGGCGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtgg                                                                             >  2:404528/1‑137 (MQ=255)
                                                                                gAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTgg                                                                 <  2:1077131/149‑1 (MQ=255)
                                                                                 aGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATc                                                                                                                 >  4:1432428/1‑100 (MQ=255)
                                                                                  gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgat                                                                                                             <  4:963476/103‑1 (MQ=255)
                                                                                  gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgat                                                                                                             >  3:1136860/1‑103 (MQ=255)
                                                                                  gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgat                                                                                                             <  4:1136860/103‑1 (MQ=255)
                                                                                  gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgat                                                                                                             >  3:963476/1‑103 (MQ=255)
                                                                                  gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTgggg                                                               >  1:545941/1‑149 (MQ=255)
                                                                                     aTCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATc                                                            >  2:1114745/1‑149 (MQ=255)
                                                                                      tCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCt                                                           <  1:1114745/149‑1 (MQ=255)
                                                                                         gTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTa                                                        >  3:1028986/1‑149 (MQ=255)
                                                                                          tGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGc                                                                                         >  2:643145/1‑115 (MQ=255)
                                                                                          tGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGc                                                                                         <  1:643145/115‑1 (MQ=255)
                                                                                              aTCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGgtggtg                                                                        >  1:399125/1‑128 (MQ=255)
                                                                                                cTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGt                                                 <  2:545941/149‑1 (MQ=255)
                                                                                                  aCTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCt                                               <  4:511394/149‑1 (MQ=255)
                                                                                                  aCTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCt                                               <  4:1028986/149‑1 (MQ=255)
                                                                                                      gCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCt                                                           >  1:1026984/1‑133 (MQ=255)
                                                                                                        tGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCttt                                         <  2:38727/149‑1 (MQ=255)
                                                                                                         gTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTg                                        >  4:366974/1‑149 (MQ=255)
                                                                                                            ggATTCGGAAGTGGATTAAGCCCGATCGATCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGtca                                     >  3:934706/1‑149 (MQ=255)
                                                                                                               ttCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGtcat                                    >  1:440943/1‑147 (MQ=255)
                                                                                                               ttCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGATggg                                                                <  1:310951/119‑1 (MQ=255)
                                                                                                               ttCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGATggg                                                                >  2:310951/1‑119 (MQ=255)
                                                                                                                  ggAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCt                                                                                               <  3:1124802/85‑1 (MQ=255)
                                                                                                                  ggAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAAc                               <  4:1286722/149‑1 (MQ=255)
                                                                                                                        ggATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCtttt                                                                                                     <  2:742421/73‑1 (MQ=255)
                                                                                                                        ggATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCtttt                                                                                                     >  1:742421/1‑73 (MQ=255)
                                                                                                                         gATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtggtg                                                                           >  3:457679/1‑98 (MQ=255)
                                                                                                                         gATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtggtg                                                                           <  4:457679/98‑1 (MQ=255)
                                                                                                                          aTTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTa                                                                                     >  3:229125/1‑87 (MQ=255)
                                                                                                                          aTTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGCGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCaaa                                >  1:105606/1‑140 (MQ=255)
                                                                                                                              aGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCt                                                                                                                       >  1:854442/1‑49 (MQ=255)
                                                                                                                               gCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCt                                                           <  4:645671/108‑1 (MQ=255)
                                                                                                                               gCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCt                                                           >  3:645671/1‑108 (MQ=255)
                                                                                                                               gCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATg                  >  4:534050/1‑149 (MQ=255)
                                                                                                                                cccGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATgg                 <  2:764868/149‑1 (MQ=255)
                                                                                                                                cccGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATgg                 >  1:666398/1‑149 (MQ=255)
                                                                                                                                 ccGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCt                                                           >  1:465026/1‑106 (MQ=255)
                                                                                                                                 ccGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCt                                                           <  2:465026/106‑1 (MQ=255)
                                                                                                                                   gATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTg                                                                                                                      >  1:1314059/1‑45 (MQ=255)
                                                                                                                                   gATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTg                                                                                                                      <  2:1314059/45‑1 (MQ=255)
                                                                                                                                    aTCGTCCCAGGAACGATCGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTgggg                                                               >  1:345948/1‑100 (MQ=255)
                                                                                                                                    aTCGTCCCAGGAACGATCGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTgggg                                                               <  2:345948/100‑1 (MQ=255)
                                                                                                                                     tCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTAt                                             <  2:311576/116‑1 (MQ=255)
                                                                                                                                        ttCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTGCTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATg                  >  4:882282/1‑140 (MQ=255)
                                                                                                                                            tGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATc     <  1:316523/149‑1 (MQ=255)
                                                                                                                                             gggACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGaaa                        <  2:735486/129‑1 (MQ=255)
                                                                                                                                             gggACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGaaa                        >  1:735486/1‑129 (MQ=255)
                                                                                                                                             gggACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGaaa                        <  2:888634/129‑1 (MQ=255)
                                                                                                                                             gggACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGa                      <  1:770406/131‑1 (MQ=255)
                                                                                                                                             gggACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCa    >  3:844204/1‑149 (MQ=255)
                                                                                                                                             gggACGAT‑GGGCTCGCTGGCAGAGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGaaa                        >  1:888634/1‑129 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATG  >  NC_000913/436559‑436850

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑