Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,228,308 G→C 14.4% pseudogene (300/1287 nt) mdtQ ← pseudogene, putative channel/filament proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,3080GC14.4% 141.4 / 2.2 70pseudogene (300/1287 nt)mdtQpseudogene, putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base C (3/7);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTG  >  NC_000913/2228083‑2228542
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                
gTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCGCCCGGGCTTTAACCGTACCCCGGCGGGGCGTAACCTCCCCCCGCTTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAg                                                                                                                                                                                                                                     >  1:1120914/1‑238 (MQ=255)
         taGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGTCGTTAAGCCAAAAGTAc                                                                                                                                                                                                                           <  2:801045/239‑1 (MQ=255)
           gCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGCTTTCTCGCGTTCCGCCGCCCGTGCTTTAACCCTCCCCAGGCGGGCCGTAAACTCCGCCCCATTCTTTCCCCCCATTTCCCGATGCCCACCCGCCGTTAAAcccacaaaacc                                                                                                                                                                                                                          >  2:551156/1‑229 (MQ=255)
            cTGGCGATCGGTCGCGATAAGGTTGTTCTGCTCTTTTTCGATTTGCTGAAAGCCCGTTTTTACCGCCGCCTGGGTTTGCCACCCCCAGTACAGGCGGGCTCCGCTGCCACCCAACAATGGGGGGGTTTCCTCGCGTCCCGCCCCCCGTGTTTTACCCGTACCCAGGGGGGCAGTACCCTCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCg                                                                                                                                                                                                                         <  2:890938/238‑1 (MQ=255)
               gCGATCGGTCGGAATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGCCCGTGTTTACCGCCGCCTGGGTTTGCCACCCCGAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGTTTTGCCCGCGTCCCGCCGCCCGTGCTTTACCCGTACCCGGGCGGGAAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTg                                                                                                                                                                                                                      <  2:1120914/238‑1 (MQ=255)
                  aTCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTTACCGTACCCCGGCGGGCAGTAACCTCCCCCCGATTCTTTCCCCCGATATCGAGATGGCACCCCGCCGTTAAGACAAAAGTACCGTTGgggt                                                                                                                                                                                                                  >  1:706395/1‑239 (MQ=255)
                          cgttattgggttttgtcctgtttcCTTTTTGCTGAAAGACCGTGTTTGACGCCGCCGGGGTTGGCAACTCCCAGGACGGGCGGGCTGCGTGCCCCGCCGCGAATTGGGGGGTTTGTCCGGGTTCCCCCGCCCGTGTTTTACCCGGCCCCAGGGGGGCTGTACCCTCCGCCGGTTTTTTCCCCCAGATTTCGGGAGGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACggg                                                                                                                                                                                                        <  2:767712/219‑1 (MQ=255)
                            aaagggggttttgtgtttttttCTATTGGTGGCAGACCCGGGTTAAGCCCCGCCGGGTTTGCCCCCCCCAGTGACGGGGGGGCTGCCCGTCCGGCCAGCTATTGGCGGGTTTGCCCGCGTCCCCCCGCCCGTGCTTTACCCTGCCCCAGGCGGGAAGAACCCTCCGCCCGATTCTTCCCCCATATACCGAAATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                                                                                                                                                      <  2:1169282/222‑1 (MQ=255)
                                     ttCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCCGCCAGCAATTGGCCGGCTTGCTCCCGTGGCCCCCCCCCTGCTTTTACCCGTCCCACGGGGGCAGCGAACCCCGCCCCCCTCTTTCCCCCACCATCCAGATGCCCGCCCCCCCCTAAGCCAAAAGCAACGTTTGTGTGACACGGGGCCGTccgac                                                                                                                                                                                             >  2:1107462/1‑234 (MQ=255)
                                     ttCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTTCGCCGCCCGTGCTTTAACCCTACCCCGGCGGGCAGTAACCTCCCCCCCGTTCTTTCCCCCGATATCGAGATGCCCGCCCCCCCGTAAGCCACAAGCACCCCTTGTGTACCACGGGccccgcgtac                                                                                                                                                                                              >  2:1085285/1‑230 (MQ=255)
                                       cTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCCCCCCCCCGGCTTTAACCCTACCCCGGCGGGCCGTAACCCCCCCCCGATTCTTTCCCCCGATATCGAGATACCAGGCCGCCCTTAAGCCCAAAACACCGTTGGTGTACCACGGGCCCGGCGTAccc                                                                                                                                                                                          >  2:1065712/1‑237 (MQ=255)
                                          ctctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCGCCGCCCGGGCTTTTACCCGACCCAGGCGGGGAGGAAGCTCCCCCCGCTTTTTTCCCCCGATTTCCGGAGGGCAAGCCCGCCGTAAACCACAAATACCACTGGTTGACCCCGGGCCCGGCGTAACTGct                                                                                                                                                                                       >  2:888835/1‑237 (MQ=255)
                                               tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTTCGCCGCCCCTGCTTTTAACGTACCCCAGCGGGCAGTAACCTCCGCCCGGTTCTTTTCCCCCATATTGAGATGCCAGCCCGCCGCTAAGCCAAAAATACCGTTGGTGTACCCCGGGCCGGGCGGACCCACGGCCCGa                                                                                                                                                                                 >  2:90209/1‑239 (MQ=255)
                                                       tgctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGGTCCGCCCCCCCTGCCTTTAACCTACCCAAGCGGGCAGGAACCTCCCCCCCGTTCTTTTCCCCGATATCGAGATGCCCGCCCCCCGTTTAACCAAAAGGACCGCTGGGGTACCCCGGGGCGGTCGTGCCTTCGGCCCGGTCGTTTag                                                                                                                                                                         >  2:745530/1‑239 (MQ=255)
                                                                gccGTGTTGACCGCCGCCGGGGTTTCCCCCTCCACCTGCCGGCGGGCTACGCCGCCACCCAGCACGGGGGGGGTTTGCTCGCGTCCCGCCCCCCGTGCTTAAACGGTCCCCGGGCGGGGAGTAACCCCCGCCCGATTTTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACgg                                                                                                                                                               <  2:701825/238‑1 (MQ=255)
                                                                   gtgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTTGCGGTCCGCCGCCCGTGCTTTAAACGAACCCAGGGGGGCAGTAACCCCCCCCCCATTTTTTTCCCCGATATCGAGATGCCCGCCCCCCCGTAAGCCAAAAATACCACTGGTGTAACACCGGGCCGGCCGACCCTCGGCCGGATCGTTccaaaccaacgcgc                                                                                                                                                              >  2:649427/1‑225 (MQ=255)
                                                                   gtgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTTCCCCGCCCCGGCTTTTACCCAAACCCCGCCGGCCGTAACCTCCCCCCGCTTCTTTTCCCCGCTATCCAGATGCCCACCCCCCCGTAAACCCAAAAGACCGCTTGTGGAACACCGGCCCGTCGTTCCCgccgcccgggcgctccaccgaaaacgcg                                                                                                                                                              >  2:387458/1‑216 (MQ=255)
                                                                    tgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTTCGCCGCCCCTGCCTTTACCGTACCCCGGCCGGCCGTAAACTCCCCCCCCATTTTTTCCCCGGTATTGAGAAGCCCACCCCCCGGTAAACCAAAAAGACCCCTTGGGTACCACCGGCCGGGCGTGACCTCGGCCCGGCCGTTCCGAGCCAACCGccc                                                                                                                                                            >  2:370020/1‑238 (MQ=255)
                                                                                ccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTTCTCCCGGTCCGCCCCCCGCGCCTTTACCCGACCCCGGCCGGGAGTAAACTCCCCCCCGCTTTTTTCCCCAACATCCAGAGGCCCACCCCCCCTTTAACCCAAAAGACCCTTTGCGTAACCCAGGGCCGGCCGTCCCTCCGCCGGGTTGGTTCGCGCAAAACGGCCCCATTAGCCTTc                                                                                                                                                >  2:275915/1‑238 (MQ=255)
                                                                                ccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGTTCGCGTTCCGCCCCCCGCGCTTTTACCCGACCCAGGGCGGGAGTAAACCCCCCCCCGTTCTTTTCCCCCCTATTGCGGTGCCCACCCCCCCGTAAACCACAAACACCCTTGGGGGACCACCGGGCCGGCCGACCCGCGGCCCGGACGTTCCGAGCAAACGGGCCCCtttacacctct                                                                                                                                               >  2:932476/1‑230 (MQ=255)
                                                                                     ttttGCCCCTCCCAGTACAGGGGGGCTACGCTGCCAGCCAGCAATTGGGGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGAAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCg                                                                                                                                           <  1:436116/237‑1 (MQ=255)
                                                                                      ttgtccccTCCCAGGCGAGGCGGGCTACGCTACCAGCCGGCAATTGGGGGTTTTGCTCACTTGCCGCGGCCCGTGCTTTAACCGTCCCCAGGGGGTCAGTAACCTCCCCCCGATTCTTCCCCCAGATATCGAGTTGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGa                                                                                                                                          <  2:863926/231‑1 (MQ=255)
                                                                                               cccAGTACGGGCGGGCTCCGCGCCCAGCCAGCAATTGGCGGGTTTGCTCGCTTCCCCCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATACCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGc                                                                                                                                <  1:223047/239‑1 (MQ=255)
                                                                                               cacAGTACAGGCGGGCTCCGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCTTTCCGCCGCCCGTGTTTTACCCGTCCCCAGGCGGGCAGTAACCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGc                                                                                                                                <  1:711313/237‑1 (MQ=255)
                                                                                                         ggCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTTTCCCGGTTCCCCGCCCCTGCTTTTACCCGACCCAAGCCGGGAGTAAACCCCCCCCCGTTCTTTTCCCCGATATCCAGAGGCCCGCCCCCCGGTTAAGCAAAAATACCCTTGGGGTACCACGGGCCGGCCGGTCCTGCGGCCGGTGCGGTCTAAGCAAACAGCCCCCATAACGCTTCTGCCCACCTTTTTTGCCGCGccttt                                                                                                                       >  2:201642/1‑236 (MQ=255)
                                                                                                             ggCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGTTCCGCCGCCCCTGGCTTTAACCTAACCCCGCGGGGAGTAAACCCCGCCCGATTTTTTTCCCCAAAATCCAGAATCCAGCCCCCCCTTAAACCCAAAAAAACGCTGGGGGACCACCGGGCCGGCGGTACCGACGCCGGATCGTTCTTAGCAAACGGCCCCCTTTAGCCTTCCTCCGCCCttttttttcccccccccatc                                                                                                                    >  2:14310/1‑224 (MQ=255)
                                                                                                                    cTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCGCCCGGGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGCTTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTACGCCAAAAGTAACCGTGGGGTACCCCGGGCCGGTCGTACCCGCGGCCCGAACGTTCAGAGCAAAAGGGCCCAATAAGCCCTCTGCCGGCCTTTTTTTCCGCCCCAAAACCGCCGcaaa                                                                                                            >  2:311032/1‑239 (MQ=255)
                                                                                                                       cccacgaaaatgggggggttttGCTCGCTTCCCGCCCCCTGTTTTTTACCCCTCCACAGGCGGACGTACACCCCCCCCCTATTCTTCCCCCAGATACCGGGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTc                                                                                                         <  2:1242777/220‑1 (MQ=255)
                                                                                                                                 aTTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGTTTTACCCGTACCCAGGCGGGCAGTAACCTCCGCCTGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGa                                                                                               <  1:732795/239‑1 (MQ=255)
                                                                                                                                    ggCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGGGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTTCCGCCCCATATCCGCCGAAAAGTCGATCTGCCGCCCa                                                                                            >  2:198529/1‑239 (MQ=255)
                                                                                                                                                                  gCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCa                                                                                                                         >  1:553862/1‑181 (MQ=255)
                                                                                                                                                                   tcttaacccGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTggg                                                                <  1:201642/231‑1 (MQ=255)
                                                                                                                                                                        aaCCGTCCCAGGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCag                                                          <  1:387458/239‑1 (MQ=255)
                                                                                                                                                                          cctcccccGGGCGGGCAGTAACCCCCGCCCGATTCTTTCCCAAGATAGCGAGATCCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagct                                                        <  1:145633/234‑1 (MQ=255)
                                                                                                                                                                           cGTACCCGGCGGGCAACTACCCTCCCCCCTATTTTTCCCCAAGATATCGAAGCGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                                                       <  1:1059770/239‑1 (MQ=255)
                                                                                                                                                                            gTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                                                       >  2:781683/1‑238 (MQ=255)
                                                                                                                                                                              aCCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                                                       >  1:173510/1‑236 (MQ=255)
                                                                                                                                                                              aCCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                                                       >  1:758394/1‑236 (MQ=255)
                                                                                                                                                                               cccAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATc                                                                                          >  1:306728/1‑200 (MQ=255)
                                                                                                                                                                               cccAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagc                                                         <  1:784653/233‑1 (MQ=255)
                                                                                                                                                                               cccAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagc                                                         >  2:784653/1‑233 (MQ=255)
                                                                                                                                                                                    gCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATgc                                              >  2:804852/1‑239 (MQ=255)
                                                                                                                                                                                       ggCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGCCATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATgc                                              >  1:18506/1‑236 (MQ=255)
                                                                                                                                                                                       ggCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATgc                                              >  2:1014849/1‑236 (MQ=255)
                                                                                                                                                                                       ggCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCt                                           >  1:225979/1‑239 (MQ=255)
                                                                                                                                                                                        ggacattacccTCCGCCCGATTCTTTCCCCAGATAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATgcg                                             <  1:275915/229‑1 (MQ=255)
                                                                                                                                                                                               acccTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGtt                                        <  1:781683/234‑1 (MQ=255)
                                                                                                                                                                                                cccTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGc                                    <  1:235010/238‑1 (MQ=255)
                                                                                                                                                                                                accccccgccgCATTCTTCCCCAAGATATCGAGAGGCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTCCCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGtt                                        <  2:96190/231‑1 (MQ=255)
                                                                                                                                                                                                accccccccccGTTTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCa                                      <  1:745530/230‑1 (MQ=255)
                                                                                                                                                                                                     ccgccgGTTTCTTGCCCAAGATATCGAGAGGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCAcc                                <  1:198529/239‑1 (MQ=255)
                                                                                                                                                                                                     ccgccgGATTCTTCCCCCAGATACCGGGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCAcc                                <  1:1019025/239‑1 (MQ=255)
                                                                                                                                                                                                     ccGCCCGATTCTTTCCCCAGTTACGGGGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCAcc                                <  1:1246289/239‑1 (MQ=255)
                                                                                                                                                                                                     ccGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCa                                  >  1:96190/1‑237 (MQ=255)
                                                                                                                                                                                                      cccccGAATTTTTCCCCAGGATACCGAAGTGCACGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCAcc                                <  1:649427/237‑1 (MQ=255)
                                                                                                                                                                                                      cGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCa                                  >  2:484460/1‑236 (MQ=255)
                                                                                                                                                                                                               cTTCCCCCAGATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCg                      <  1:932476/239‑1 (MQ=255)
                                                                                                                                                                                                                tttCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGCTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCgg                     >  1:746114/1‑239 (MQ=255)
                                                                                                                                                                                                                 tcccccAGATAGCGGAATGCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGt                    <  2:928752/237‑1 (MQ=255)
                                                                                                                                                                                                                 tcccccAAATATGGAGGTTCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGGCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGt                    <  1:1196060/237‑1 (MQ=255)
                                                                                                                                                                                                                  tCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGt                    >  1:976107/1‑238 (MQ=255)
                                                                                                                                                                                                                    cccAGATTGGGAGATGCCTGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTCCCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCg                 <  1:14310/239‑1 (MQ=255)
                                                                                                                                                                                                                    cccAGATATCGAGATGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCg                 <  1:92495/239‑1 (MQ=255)
                                                                                                                                                                                                                    cccAGATACCGAGATGCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCg                 <  1:1065712/239‑1 (MQ=255)
                                                                                                                                                                                                                          aaatCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTg            <  1:370020/237‑1 (MQ=255)
                                                                                                                                                                                                                              agcgagcgcccAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTa          <  2:173510/230‑1 (MQ=255)
                                                                                                                                                                                                                               gagaTGCCACCCCGCCGTTAAGCCAAAAGCACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGt        <  1:1085285/238‑1 (MQ=255)
                                                                                                                                                                                                                                     gCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCAtt   >  2:342504/1‑239 (MQ=255)
                                                                                                                                                                                                                                     gCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTg  >  1:589036/1‑239 (MQ=255)
                                                                                                                                                                                                                                     gCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGAAGTGCATTg  >  2:1053926/1‑239 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                
GTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGTGCATTG  >  NC_000913/2228083‑2228542

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: