Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,233,154 T→A 7.5% intergenic (‑124/+22) dsbB ← / ← nhaB oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I/sodium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,233,1540TA 7.5% 733.4 / 37.8 214intergenic (‑124/+22)dsbB/nhaBoxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (90/108);  new base A (8/8);  total (98/116)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

Only 100 of 222 total aligned reads displayed.

AGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGC  >  NC_000913/1232903‑1233412
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                  
agCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATNATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGAANCAAACATAGGAAGCAATTGGACAATTATCGGGCGAAAAGGCAAACNGNGCATTTAAATGCCCGGGAAggtgttta                                                                                                                                                                                                                                              >  3:39274/1‑272 (MQ=255)
                                    cacaGTGCCGTCATTTCAAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGACCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCAtt                                                                                                                                                                                                         <  1:39274/275‑1 (MQ=255)
                                       aGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGcc                                                                                                                                                                                                                                                            >  4:50903/1‑221 (MQ=255)
                                                                                                       aaCATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCAtt                                                                                                                                                                                                                                                                    >  6:772336/1‑149 (MQ=255)
                                                                                                       aaCATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAATCCGACGAGTTTCAGGACGAGGGGGTAAg                                                                                                                                       >  2:224625/1‑274 (MQ=255)
                                                                                                        aCATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCAttt                                                                                                                                                                                                                                                                   >  8:1022587/1‑149 (MQ=255)
                                                                                                         cATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTta                                                                                                                                                                                                                                                                  <  5:1334304/149‑1 (MQ=255)
                                                                                                          aTTGGTTCAAAAATCTCTACATAATCATTTGCCTGTAATATGCATAGAGTGCGCAGTTTAAACCAATTGATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTCTCAGGCGAAAAGGCAAACCGGGCATTtat                                                                                                                                                                                                                                                                 >  8:903071/1‑149 (MQ=255)
                                                                                                          aTTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtat                                                                                                                                                                                                                                                                 >  8:182467/1‑149 (MQ=255)
                                                                                                          aTTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAgg                                                                                                                                                    >  4:230758/1‑258 (MQ=255)
                                                                                                           ttGGTTCAAAAATTGCAACATAATCATTTCCCTGCAATATGCATAGAGAGGTCAGTTTAAACCAATTCATGCGCTGCGCCACCAGGAGGCAGCAAAAATCGGACGCAATTGGACAATTATCGAACGAAAAGGCAAACCGGGCATTAATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGTAGAGCAAGCCGACGAGTGGCAGGACGAGGGTGTAAGGCAgg                                                                                                                                  <  4:127271/275‑1 (MQ=255)
                                                                                                           ttGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                                                                                                                                                                                >  6:91969/1‑149 (MQ=255)
                                                                                                           ttGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                                                                                                                                                                                <  5:915753/149‑1 (MQ=255)
                                                                                                            tGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTatat                                                                                                                                                                                                                                                               <  8:1001035/149‑1 (MQ=255)
                                                                                                            tGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTatat                                                                                                                                                                                                                                                               <  5:374649/149‑1 (MQ=255)
                                                                                                             ggTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATg                                                                                                                                                                                                                                                              <  5:1373853/149‑1 (MQ=255)
                                                                                                               ttCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATGGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGcc                                                                                                                                                                                                                                                            >  5:1213347/1‑149 (MQ=255)
                                                                                                               ttCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGcc                                                                                                                                                                                                                                                            >  7:308137/1‑149 (MQ=255)
                                                                                                               ttCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGcc                                                                                                                                                                                                                                                            <  8:1053198/149‑1 (MQ=255)
                                                                                                                  aaaaaTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCgg                                                                                                                                                                                                                                                         >  5:369927/1‑149 (MQ=255)
                                                                                                                   aaaaTCGCAACCTAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGt                                                                                                                                                                                                                                                        <  5:1402422/149‑1 (MQ=255)
                                                                                                                   aaaaTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGt                                                                                                                                                                                                                                                        <  6:1240779/149‑1 (MQ=255)
                                                                                                                    aaaTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                                                                                                                                                                                       >  5:654626/1‑149 (MQ=255)
                                                                                                                       tCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgt                                                                                                                                                                                                                                                    <  7:916557/149‑1 (MQ=255)
                                                                                                                        cGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgtt                                                                                                                                                                                                                                                   <  6:969138/149‑1 (MQ=255)
                                                                                                                          cAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttgt                                                                                                                                                                                                                                                 >  6:944649/1‑149 (MQ=255)
                                                                                                                           aaCATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGTCAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGt                                                                                                                                                                                                                                                        >  8:889320/1‑141 (MQ=255)
                                                                                                                           aaCATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGGAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttgtt                                                                                                                                                                                                                                                >  6:8189/1‑149 (MQ=255)
                                                                                                                           aaCATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGt                                                                                                                                                                                                                                                        >  7:902581/1‑141 (MQ=255)
                                                                                                                           aaCATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttgtt                                                                                                                                                                                                                                                >  7:143518/1‑149 (MQ=255)
                                                                                                                           aaCATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttgtt                                                                                                                                                                                                                                                <  7:599414/149‑1 (MQ=255)
                                                                                                                             cATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTAt                                                                                                                                                                                                                                              >  5:381673/1‑149 (MQ=255)
                                                                                                                                 aTCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCaaa                                                                                                                                                                                                                                          <  7:172443/149‑1 (MQ=255)
                                                                                                                                   cATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGc                                                                                                                                                                                                                                        <  5:440508/149‑1 (MQ=255)
                                                                                                                                     tttCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGt                                                                                                                                                                                                                                      <  5:493820/149‑1 (MQ=255)
                                                                                                                                      ttCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGtt                                                                                                                                                                                                                                     <  5:1030718/149‑1 (MQ=255)
                                                                                                                                        cccTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTCATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                                                                                                                                                                                   <  6:511268/149‑1 (MQ=255)
                                                                                                                                         ccTGCAATAGGCATAGAGTGCACAGTTTAAACCAATTCATCGGCTGCGCCACCAGGAGGCAGAAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAgg                                                                                                    <  4:224625/275‑1 (MQ=255)
                                                                                                                                          cTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCACGCGTACACTCGACTCAGAGCAAGCCGCCGAGTTGCAGGTCGAGGGTTTACGGCAGGGCCATCCACACCATGGCGCACAAAGAGAgg                                                                                                    >  4:280949/1‑274 (MQ=255)
                                                                                                                                           tGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTAt                                                                                                                                                                                                                                <  7:495228/149‑1 (MQ=255)
                                                                                                                                                   gCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCa                                                                                                                                                                                                                        <  5:704389/149‑1 (MQ=255)
                                                                                                                                                   gCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCa                                                                                                                                                                                                                        >  6:283142/1‑149 (MQ=255)
                                                                                                                                                    cATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAt                                                                                                                                                                                                                       >  8:435541/1‑149 (MQ=255)
                                                                                                                                                         agTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCa                                                                                                                                                                                                                  <  6:1030535/149‑1 (MQ=255)
                                                                                                                                                           tGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATa                                                                                                                                                                                                                <  7:247019/149‑1 (MQ=255)
                                                                                                                                                            gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTAt                                                                                                                                                                                                                                >  4:370561/1‑132 (MQ=255)
                                                                                                                                                             cgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATaaa                                                                                                                                                                                                              <  8:754876/149‑1 (MQ=255)
                                                                                                                                                             cgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATaaa                                                                                                                                                                                                              <  5:1471768/149‑1 (MQ=255)
                                                                                                                                                             cgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATaaa                                                                                                                                                                                                              <  7:1365217/149‑1 (MQ=255)
                                                                                                                                                                aGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAg                                                                                                                                                                                                                            <  2:370562/132‑1 (MQ=255)
                                                                                                                                                                  tttAAACCAATTCATTCGCTGCGCCTCCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTTTCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCAtt                                                                                                                                                                                                         >  5:1431043/1‑149 (MQ=255)
                                                                                                                                                                   ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                                                                                                                                                                                        <  5:975516/149‑1 (MQ=255)
                                                                                                                                                                         cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTAc                                                                                                                                                                                                  >  8:50615/1‑149 (MQ=255)
                                                                                                                                                                         cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTAc                                                                                                                                                                                                  >  5:1010359/1‑149 (MQ=255)
                                                                                                                                                                         cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTAc                                                                                                                                                                                                  >  5:231288/1‑149 (MQ=255)
                                                                                                                                                                            ttcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACCCCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTcagaag                                                                    >  2:268914/1‑272 (MQ=255)
                                                                                                                                                                            ttcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGGGTGCAGAGGTcagcag                                                                    >  2:287898/1‑272 (MQ=255)
                                                                                                                                                                             tcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAggg                                                                                                                                                                                              <  5:323274/149‑1 (MQ=255)
                                                                                                                                                                              cattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgggg                                                                                                                                                                                             >  5:866257/1‑149 (MQ=255)
                                                                                                                                                                               attcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGc                                                                                                                                                                                            <  5:499446/149‑1 (MQ=255)
                                                                                                                                                                                ttcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATcc                                                                                                                                                                                                                              >  6:1476452/1‑114 (MQ=255)
                                                                                                                                                                                ttcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCa                                                                                                                                                                                           >  6:761455/1‑149 (MQ=255)
                                                                                                                                                                                ttcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCa                                                                                                                                                                                           >  7:746583/1‑149 (MQ=255)
                                                                                                                                                                                 tcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGcaggaaca                                                               <  2:58352/272‑1 (MQ=255)
                                                                                                                                                                                  cGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAg                                                                                                                                                                                         >  8:1242686/1‑149 (MQ=255)
                                                                                                                                                                                   gCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGc                                                                                                                                                                                        <  6:534888/149‑1 (MQ=255)
                                                                                                                                                                                       cgcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAg                                                      >  2:171775/1‑275 (MQ=255)
                                                                                                                                                                                         ccaccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGAAAGTTGTTATCAAAGCGTAGCTATCCAGCCCATTTGCATAAATCATTCGGTTACAGGGGCAGGAGTAAAc                                                                                                                                                                                  >  8:577621/1‑149 (MQ=255)
                                                                                                                                                                                           accaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTc                                                                                                                                                                                >  8:664019/1‑149 (MQ=255)
                                                                                                                                                                                           accaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTc                                                                                                                                                                                >  8:460761/1‑149 (MQ=255)
                                                                                                                                                                                             caGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                                                                                                                                                                                        <  5:448977/123‑1 (MQ=255)
                                                                                                                                                                                             caGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                                                                                                                                                                                        >  6:448977/1‑123 (MQ=255)
                                                                                                                                                                                                     agcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGc                                                                                                                                                                      >  7:1415462/1‑149 (MQ=255)
                                                                                                                                                                                                       caAAAATAGGACGCAATTGGACAATTATCGGGCTAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCaa                                                                                                                                                                    <  8:964626/149‑1 (MQ=255)
                                                                                                                                                                                                       caAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCaa                                                                                                                                                                    <  7:970302/149‑1 (MQ=255)
                                                                                                                                                                                                        aaaaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAt                                                                                                                                                                                                                       >  7:650309/1‑97 (MQ=255)
                                                                                                                                                                                                               ggACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCcgacga                                                                                                                                                            <  7:460761/149‑1 (MQ=255)
                                                                                                                                                                                                               ggACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCcgacga                                                                                                                                                            <  5:193864/149‑1 (MQ=255)
                                                                                                                                                                                                               ggACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCcgacga                                                                                                                                                            >  8:161711/1‑149 (MQ=255)
                                                                                                                                                                                                                     aaTTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCa                                                                                                                                                      >  5:536430/1‑149 (MQ=255)
                                                                                                                                                                                                                      aTTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAg                                                                                                                                                     <  6:440949/149‑1 (MQ=255)
                                                                                                                                                                                                                        tGGAAAAGTATCGGGACAAAAGGAAAACCGGGAATTTATATGCCCGGTAAGTTGTTATCAAAGCGTGGCTATCCAGCCCATGTGCATAAACCATTCGGTTACATGGGCAAGCGTAAACTCGACGCAGGGCAAGCCGACGAGTGTCAGGACGAGGGTGGAAGGCAGGGCCAGCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGAt                     <  4:6204/275‑1 (MQ=255)
                                                                                                                                                                                                                            cAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAg                                                                                                                                               <  7:1274537/149‑1 (MQ=255)
                                                                                                                                                                                                                             aaTTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCCAAGCGTTGCTATCCAGCCCATTTGCATAAACCCTTCGGTTACAGGTGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAgg                                                                                                                                              <  8:1389093/149‑1 (MQ=255)
                                                                                                                                                                                                                             aaTTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAgg                                                                                                                                              <  8:547779/149‑1 (MQ=255)
                                                                                                                                                                                                                             aaTTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAgg                                                                                                                                              >  5:763424/1‑149 (MQ=255)
                                                                                                                                                                                                                               ttATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGgt                                                                                                                                            <  7:372421/149‑1 (MQ=255)
                                                                                                                                                                                                                               ttATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGgt                                                                                                                                            <  6:373637/149‑1 (MQ=255)
                                                                                                                                                                                                                               ttATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCAGTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTAcc              <  4:294810/275‑1 (MQ=255)
                                                                                                                                                                                                                                   cGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag                                                                                                                                                                         <  6:596874/116‑1 (MQ=255)
                                                                                                                                                                                                                                   cGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTaa                                                                                                                                        <  7:869894/149‑1 (MQ=255)
                                                                                                                                                                                                                                   cGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATg                                                                                                                   <  2:71064/170‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAgg                                                                                                                                  >  5:637290/1‑149 (MQ=255)
                                                                                                                                                                                                                                              gCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTaa                                                                                                                                        <  8:414518/138‑1 (MQ=255)
                                                                                                                                                                                                                                              gCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTaa                                                                                                                                        >  7:414518/1‑138 (MQ=255)
                                                                                                                                                                                                                                              gCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCAt                                                                                                                             >  5:1180323/1‑149 (MQ=255)
                                                                                                                                                                                                                                              gCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCAt                                                                                                                             >  6:1363669/1‑149 (MQ=255)
                                                                                                                                                                                                                                              acttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAg                                                                                                                                               <  7:574726/127‑1 (MQ=255)
                                                                                                                                                                                                                                               cttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATc                                                                                                                            <  6:399779/146‑1 (MQ=255)
                                                                                                                                                                                                                                               cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATc                                                                                                                            <  6:1010359/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                  
AGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCATAAGAGAGGCGAATCAATGGCGCGAGTGCAGAGGTCAGCAGGAACAGGAACGCAGCCTGACCGTTCGGCGTAGCGACGGAGGGCAGATTGGTACCGGTATTAATGGC  >  NC_000913/1232903‑1233412

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