Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,554,539 T→G Y485* (TAT→TAG)  malT → DNA‑binding transcriptional activator MalT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,554,5390TG97.0% 218.5 / ‑2.0 66Y485* (TAT→TAG) malTDNA‑binding transcriptional activator MalT
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (30/34);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

ACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTATAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCATCTGGAACAGCTGCC  >  NC_000913/3554292‑3554784
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
aCAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGc                                                                                                                                                                                                                                                    >  1:101521/1‑251 (MQ=255)
       gCTGGAAAAGCCGCAGTTGGTGGGAGTGCAGGCGGGGCTGAGGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCGGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGAGTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGAAAAACTGGCACTGGAAGAGCTGACGCCGGGCTGGTTCTAGAGCCGCATtg                                                                                                                                                                                                                                             <  2:79325/251‑1 (MQ=255)
         tGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCCAAACTGGCACTGGAAGAGCTGCCCCCGGGCTGGTTCTAGAGCCGCATtgtg                                                                                                                                                                                                                                           >  2:4913/1‑251 (MQ=255)
            aaaaTCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGc                                                                                                                                                                                                                                         >  2:202022/1‑250 (MQ=255)
             aaaTCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAAGTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCGGCCGCCGGGCTGGTTCGAGAGCCGCATTGTGGCaa                                                                                                                                                                                                                                       <  2:26684/251‑1 (MQ=255)
                     aGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGGGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGGAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTg                                                                                                                                                                                                                               >  2:146252/1‑251 (MQ=255)
                         ggTGTTAGTGCAGGCGTGGCTGATGCAAAGCCAACAGCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAATACACCATGCACGCAGAAGTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCGGCCGCCGGTCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTgg                                                                                                                                                                                                                           <  1:170126/251‑1 (MQ=255)
                          gtgtTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTCGCCCGTGCTGCACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTggg                                                                                                                                                                                                                          >  2:259471/1‑251 (MQ=255)
                          gtgtTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTggg                                                                                                                                                                                                                          <  1:50689/251‑1 (MQ=255)
                                     ggCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGtgctgc                                                                                                                                                                                                               <  1:157506/251‑1 (MQ=255)
                                      gCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAATGTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACGCCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAAGCCGGATGAAGCGGAACGGCTGTGAAAACTGGCACTGGGAGAGCTGCCGCCGGGTTGGTTCTAGAGTCGTACTGGGGCAACCTCGGTGCTGGGTGAAGTGctgca                                                                                                                                                                                                              <  1:252041/251‑1 (MQ=255)
                                          ggCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGctgcactgc                                                                                                                                                                                                          >  2:149177/1‑251 (MQ=255)
                                                  cAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGa                                                                                                                                                                                                  <  1:18393/251‑1 (MQ=255)
                                                   aaaGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGaa                                                                                                                                                                                                 <  1:200557/251‑1 (MQ=255)
                                                    aaGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCGCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAAt                                                                                                                                                                                                <  2:203773/251‑1 (MQ=255)
                                                                     gCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTATAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGgcg                                                                                                                                                                               >  1:101404/1‑251 (MQ=255)
                                                                       gAAGTTAACACCCTGCTAGCCCGTGCTGAACATGACATCAAGGAAATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTACTCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCt                                                                                                                                                                             >  1:93748/1‑251 (MQ=255)
                                                                         aGTTAACACCCTGCTAGCCCGCGCTGACCATGACCTCAAGGACATAAGAGAAGACACAATGCCCGCAGACTTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAACCGGAAAGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCACAGGCGAATTGACCCGCTCACTGGCGCTa                                                                                                                                                                            >  2:175104/1‑250 (MQ=255)
                                                                            tAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATgc                                                                                                                                                                        >  1:95096/1‑251 (MQ=255)
                                                                               cacCCTGCTAGCCCGTTCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCGCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGGGCTGCACTGCAAAGGCTAATTGACCCGCTCACTGGCGCTAATgcagc                                                                                                                                                                     <  2:170867/251‑1 (MQ=255)
                                                                                acCCTGCTAGCCCGTGCTGCCCATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTATGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATgcagc                                                                                                                                                                     >  2:243135/1‑250 (MQ=255)
                                                                                     gCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCg                                                                                                                                                               >  1:93422/1‑251 (MQ=255)
                                                                                      cTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCg                                                                                                                                                               >  1:159035/1‑250 (MQ=255)
                                                                                      cTAGCCCGTGCTGAACATGAAATCAAGGACATAAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCg                                                                                                                                                               >  1:128468/1‑250 (MQ=255)
                                                                                        aGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGaa                                                                                                                                                             <  1:213802/250‑1 (MQ=255)
                                                                                           ccGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAAc                                                                                                                                                            >  1:83984/1‑248 (MQ=255)
                                                                                             gTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATg                                                                                                                                                       <  1:198359/251‑1 (MQ=255)
                                                                                                     cATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGGGAAGTGCTGCACTGCAACGGCGAATTGACCCGCTCACTGGCGCTACTGCAGCAAACCGCACAGATGGCACGccc                                                                                                                                               >  2:108193/1‑250 (MQ=255)
                                                                                                             cAAGGACATCAGAGAAGACACCATGCACGCAGAATTGAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCGGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCACCCTCGGTGCTGGGTGAAGTGCTGCACGGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACTCCAGCACGAt                                                                                                                                        <  2:176794/250‑1 (MQ=255)
                                                                                                              aaGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGAt                                                                                                                                                        <  2:257554/233‑1 (MQ=255)
                                                                                                              aaGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGAt                                                                                                                                                        >  1:257554/1‑233 (MQ=255)
                                                                                                              aaGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCCCGCCAGCCCGATGt                                                                                                                                      >  2:137082/1‑251 (MQ=255)
                                                                                                                ggACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCt                                                                                                                                    <  1:263858/251‑1 (MQ=255)
                                                                                                                 gACATCAGAGAAGACACCATGCACGAAGAATGTAACGCGCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGAGGAAGCGGAACGGCTGGCAAAACTGGAACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTg                                                                                                                                   <  2:62229/251‑1 (MQ=255)
                                                                                                                      cAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACt                                                                                                                              >  2:101493/1‑251 (MQ=255)
                                                                                                                         agaagaCACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACg                                                                                                                           <  1:236736/251‑1 (MQ=255)
                                                                                                                            agaCACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTATAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCtt                                                                                                                        <  2:101404/251‑1 (MQ=255)
                                                                                                                             gaCACCATGCACGCAGAAGTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATTGCACGCCAGCACGATGGCTGGCACTACGCttt                                                                                                                       <  1:237647/251‑1 (MQ=255)
                                                                                                                               cacCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTtgg                                                                                                                     >  1:75437/1‑250 (MQ=255)
                                                                                                                                   aTGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTGACTGGCGCTAATGCAGCAAACCGAAAAGATGGCACGCCAGCACGATGTCTGGCACAACGCTTtgagg                                                                                                                   >  1:149130/1‑246 (MQ=255)
                                                                                                                                        cgcaggttttAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGAAATCCGGATGAAGCGGAACGGCTGGCAAAACTTGCACTGGAAGAGCTGCCGCCGGGATGGTTCTAGAGCCGCATTGTGGCGACCTCGGTGCTGGGGGAAGTGCTGCACTGCAAGGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAGCCAAACGGATGGCACGCCAGCACGAGGACTGGCACTACGCTTtgtagcgtgt                                                                                                              >  1:71190/8‑242 (MQ=255)
                                                                                                                                          cAGAATTTAACGCTCGGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATc                                                                                                          <  2:11201/251‑1 (MQ=255)
                                                                                                                                             aaTTTAACGCTCTGCGCGCCCAGGGGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAATTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGTTAATGCAGCAAACCGAACAGATGGCACGCCCGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAg                                                                                                       <  1:25107/251‑1 (MQ=255)
                                                                                                                                                     gCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTg                                                                                                >  2:83641/1‑250 (MQ=255)
                                                                                                                                                          gcgcgcCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAAAGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCGCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGaaa                                                                                             <  2:204082/248‑1 (MQ=255)
                                                                                                                                                           cgcgcCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCt                                                                                         <  1:101493/251‑1 (MQ=255)
                                                                                                                                                                    atgccGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTTGAAGAGCTGCCGCCGGGCTGGGTCTAGAGCCGCATTGTTGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCTCTAATTCAACAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCACAGTGAAATTCTGTTTGccc                                                                                 >  2:79333/5‑250 (MQ=255)
                                                                                                                                                                      ggCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAg                                                                              <  1:202022/251‑1 (MQ=255)
                                                                                                                                                                         gATTAACGATGGTAAGCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAAGGCAGCAAACCGAACAGATGGCACGCCAGCACGATGGCTGGCACGACGCTTTGTGGAGTTGAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGt                                                                           <  1:104532/251‑1 (MQ=255)
                                                                                                                                                                                   ggTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAAc                                                                 <  1:198156/251‑1 (MQ=255)
                                                                                                                                                                                             aTGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGa                                                        >  1:67961/1‑250 (MQ=255)
                                                                                                                                                                                                      aaCGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAACGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTCCGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGACACGCAGGa                                               >  2:90582/1‑250 (MQ=255)
                                                                                                                                                                                                      aaCGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAggg                                               >  1:120114/1‑249 (MQ=255)
                                                                                                                                                                                                             ggCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCAt                                       >  1:265113/1‑251 (MQ=255)
                                                                                                                                                                                                               cAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTc                                     <  2:208743/251‑1 (MQ=255)
                                                                                                                                                                                                                 aaaCTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTcc                                    >  1:122011/1‑250 (MQ=255)
                                                                                                                                                                                                                         acgggAAGAGCTGCCGCCGGGCTGGGTCTAGAGCCGCATTGTGGCAACCGCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCa                           <  2:191467/248‑1 (MQ=255)
                                                                                                                                                                                                                            ggAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACg                        <  1:259471/251‑1 (MQ=255)
                                                                                                                                                                                                                                gagCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGagca                    <  1:149177/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gCTGCCGCCGGGCTGGTTCTAGAGCCGCATTATGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGagcagc                  <  1:243135/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cgccgGGCTGGTTCTAGAGCCGCATTGTGGCAACCTAGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCATCTg             >  1:24068/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cgGGCTGGTTCTAGAGCCGCATTGTGGCAACGTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCATCTGGaa          <  1:169368/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cgGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCGGATCAACGGGCAGCGTCTGGaa          <  2:37138/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ggCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCATCTGGAACa        <  2:122011/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 gTTCTAGAGCCGCATTGGGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCACACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGTTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCATCTGGAACAGCTGc   <  2:55676/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAGCGAGCAGCATCTGGAACAGCTGcc  >  2:198830/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
ACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTATAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCATCTGGAACAGCTGCC  >  NC_000913/3554292‑3554784

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