Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,237,105 G→T P217T (CCG→ACG)  mglC ← D‑galactose/methyl‑galactoside ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,237,1050GT98.0% 137.6 / NA 49P217T (CCG→ACG) mglCD‑galactose/methyl‑galactoside ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (30/18);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGA  >  NC_000913/2236864‑2237350
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
gTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCCGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTCCCGcc                                                                                                                                                                                                                                                >  2:235086/1‑251 (MQ=255)
 tAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGcc                                                                                                                                                                                                                                                >  1:58734/1‑250 (MQ=255)
            aCCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGCTACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGATCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAAt                                                                                                                                                                                                                                    >  2:77462/1‑251 (MQ=255)
                 ggTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGCACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCt                                                                                                                                                                                                                               >  1:14849/1‑251 (MQ=255)
                  gacaccacGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCtt                                                                                                                                                                                                                              >  2:404721/3‑251 (MQ=255)
                     accacGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTAcc                                                                                                                                                                                                                           >  1:38651/1‑251 (MQ=255)
                     accacGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTAcc                                                                                                                                                                                                                           >  1:54638/1‑251 (MQ=255)
                     accacGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATACCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTAcc                                                                                                                                                                                                                           >  2:11550/1‑251 (MQ=255)
                         cgccaacccccccccccACCCCGCCGCGGAACGATCCCCCGCCTCCGACGCACGCGGCAATA‑GCATCCAGCGCATACATAAAGCCGAGGTTTTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACTCGCCAGAAAACGCGTAAATCAGCAGCAGGTTCAGGCCTAC‑GTTGACACCAGATTCTTTTGCCGCTGCCGTGGTACCGCCAAGGGCAAAAATGTTCTTACCGAAg                                                                                                                                                                                                                       <  1:376895/238‑1 (MQ=255)
                          gCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGc                                                                                                                                                                                                                      >  1:322080/1‑251 (MQ=255)
                            caaccccccccccccccccgccgCTGACCGATCCCCCGCCTACCCCGCACCCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATCCGTCCGGCTTCAAACACCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGACCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTTTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCgg                                                                                                                                                                                                                    <  1:170661/234‑1 (MQ=255)
                                 cccgaccccACCCCGCCGCTGAACGATCCCCCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGtttt                                                                                                                                                                                                               <  1:330336/246‑1 (MQ=255)
                                                     gAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACg                                                                                                                                                                                           >  1:405284/1‑250 (MQ=255)
                                                       ccGATACCCCGCCTCCCACGCCCGCCGCGATA‑GCACCCAGCTCATACATAAACCCGAGGTTGTGGGTGGCAGAGCCGATACGCCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGTTATTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAAc                                                                                                                                                                                        <  2:379582/250‑1 (MQ=255)
                                                                 gCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCa                                                                                                                                                                              >  1:218548/1‑250 (MQ=255)
                                                                     aCCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCTCCAAAGTCATAGAACTCTCCAGACAACTCGTAGATCATCAGCAGGTTCATTCCGAC‑GTTGACACCAGATATTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACTT‑CACCGCAATc                                                                                                                                                                           <  1:383972/249‑1 (MQ=255)
                                                                         cgcacgcCGCG‑TACGCATCCAGCTCATAAATAAAGCCGGGGTTGTGGGTGGCAGAGCCGATACGTCCGGCTTCTAACACCCCGCAAAAGGCAGAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACCCCAGATCCTTTTGCCGCTTCCGTGTTCCCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAAcgcg                                                                                                                                                                     <  2:170442/251‑1 (MQ=255)
                                                                             cgccgcGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAAC‑CCCACGACGAACGC‑CACCGCAATCCACGCGt                                                                                                                                                                    >  1:211712/1‑248 (MQ=255)
                                                                             cgccgcGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGACAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGt                                                                                                                                                                    >  2:88124/1‑248 (MQ=255)
                                                                                   gATA‑GCATCCAGCTCATACATAAACCCGAGGTTGTTGGGGGCAGAGCCGATACTTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACCCGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCATATACTTTGGCCCCTTCCGTTTTTCCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCTT‑GACGA‑CGTTAAGCGCAATCAACGCGTAGAAGGTGa                                                                                                                                                           <  1:24065/251‑1 (MQ=255)
                                                                                         cacccAGCTCATACAAAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTGGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTaa                                                                                                                                                      <  1:404721/248‑1 (MQ=255)
                                                                                         cATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGGTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCACcgcgaagcaggcgaggta                                                                                                                                                       >  1:115569/1‑236 (MQ=255)
                                                                                             cAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGGGGAGAAGGGGATGTAag                                                                                                                                                     >  2:222033/1‑248 (MQ=255)
                                                                                               gCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCACGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGGGATGTAagcg                                                                                                                                                   >  2:367484/1‑246 (MQ=255)
                                                                                                 aaaaacaaaaaGCCGAGTTTGTTGGGGACAGCGCCGATACGCCCGGCTACTAACACCCCGCCAAAGGCATCAACCACGCCAGAAAACGCGAAGATCATCAGCAGGTTCAGGCCGACTGTT‑ACACCAGATACTTTTGCCGCTTCCGTGGTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTGTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACg                                                                                                                                              <  2:62046/243‑1 (MQ=255)
                                                                                                            cccgatttggtgggtggCAGAGCCGATACGCCCGGCTTCTAACACCCCGCCAAAGGCATAGACCACGCCAGAAAACGCGTAGACCATCAGCAGGTTCAGGCCGAC‑GTTGACCCCAGATCCTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTccc                                                                                                                                     <  2:82313/237‑1 (MQ=255)
                                                                                                              cGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAgc                                                                                                                                 >  1:93490/1‑251 (MQ=255)
                                                                                                                  gttgttGGTGGCAGAGCCGATCCGTCCGGCTTCTAACATCCCTCCAAAGGCATAGAACCCGCCAGACAACGCGTAGATCATCAGCGGGTTCAGGCCGCC‑GTTGACTCCAGATCCTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGa                                                                                                                             <  1:85276/251‑1 (MQ=255)
                                                                                                                    tgttgGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTccc                                                                                                                                     >  1:82313/1‑241 (MQ=255)
                                                                                                                    tgtgggtggCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCTTCAGCAGTATCAGGCCGCC‑GTTGACACCAGATACTTTTGCCGCTTCCTTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGAc                                                                                                                            <  1:426923/246‑1 (MQ=255)
                                                                                                                                    gATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAgg                                                                                                           >  1:236169/1‑251 (MQ=255)
                                                                                                                                                 tCTAACATCCCGCCAAAGGCATAGCACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGGTCTTACCGAAGCGGGTTTTGTTCCACAACACCAA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGCGATGTAAGAGAGACGGAAACTCCCCACCGCGACAAAGCCCTGCGCAAAGGGAAAGAAGCCACt                                                                                              >  2:307268/1‑251 (MQ=255)
                                                                                                                                                                               ccAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGACCGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGAGCGACAAAGCCCTGAGCAACGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGAc                                                                >  2:193465/1‑251 (MQ=255)
                                                                                                                                                                                 aGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACa                                                               >  1:261348/1‑250 (MQ=255)
                                                                                                                                                                                    cAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCAGCCAAAgg                                                           >  2:404996/1‑250 (MQ=255)
                                                                                                                                                                                      aCGCGTAGATCATCGGCAGGTTCAGGCCGAC‑GTTGACACCAGATTCTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCg                                                         <  1:193465/251‑1 (MQ=255)
                                                                                                                                                                                         ggtagacaatcaccggtttcggtcagAC‑GTGGACACAAGTTCCTTTTGCCGCTTCCTTGTTCCCGCCAATGGCAAAAATGTTCTTGCCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAAAGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGTCCCTACAAAGTCGTaa                                                      <  2:249922/226‑1 (MQ=255)
                                                                                                                                                                                         cgTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGCCAAAGTCGTa                                                       >  2:320587/1‑250 (MQ=255)
                                                                                                                                                                                                  ccagcagGTTCAGGCCGCT‑GTTGACACCAGATCCTTTTCCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAg                                             <  2:132638/250‑1 (MQ=255)
                                                                                                                                                                                                        ggTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATg                                       >  1:328883/1‑251 (MQ=255)
                                                                                                                                                                                                             aGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCAt                                   >  1:89984/1‑250 (MQ=255)
                                                                                                                                                                                                               gCCGAC‑TTTGACACCAGATACTTTTGCCGCTCCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAAGTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGa                                <  1:367484/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCCCAAATGTTCTTACCGAAGCGGGTTTTGTTCCACCACACCCA‑GACGCACGC‑CACCGCAATCAACGCGTCGAAGGTGATGTAAGAGCGACGGACACTCCCCAGCGCGACCAAGCCCTGAGCAAAGGTAGAAAAGCCACTGTCCAAACCAGAAATTGGCGACGCCCCGACAAAGTAGTAATAGAGCGCGTTGATGACATAGACg                              >  1:284688/1‑251 (MQ=255)
                                                                                                                                                                                                                       ttGACGCCAGATCCTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACgatg                           <  1:77462/249‑1 (MQ=255)
                                                                                                                                                                                                                                 aTACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCg               >  1:14374/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aCTTATGCCGCTTCCGGGTAACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAg             <  2:38651/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGt       >  2:31904/1‑250 (MQ=255)
                                                                                                                                                                                                                                              ttCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTACACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCAACTGTCAAACCCAGAAATTGGCGACGCCCCGACAACGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATg   >  1:68490/1‑250 (MQ=255)
                                                                                                                                                                                                                                               tCCGTGTTACCGCCAATGGCAAAAATGTGCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCTAGCGTGGTAATGa  <  2:61880/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
GTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATA‑GCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGAC‑GTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCA‑GACGAACGC‑CACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGA  >  NC_000913/2236864‑2237350

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