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

GTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGA  >  NC_000913/2236864‑2237350
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
gTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCCGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTCCCGcc                                                                                                                                                                                                                                                >  2:235086/1‑251 (MQ=255)
 tAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGcc                                                                                                                                                                                                                                                >  1:58734/1‑250 (MQ=255)
            aCCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGCTACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGATCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAAt                                                                                                                                                                                                                                    >  2:77462/1‑251 (MQ=255)
                 ggTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGCACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCt                                                                                                                                                                                                                               >  1:14849/1‑251 (MQ=255)
                  gacaccacGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCtt                                                                                                                                                                                                                              >  2:404721/3‑251 (MQ=255)
                     accacGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTAcc                                                                                                                                                                                                                           >  1:38651/1‑251 (MQ=255)
                     accacGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTAcc                                                                                                                                                                                                                           >  1:54638/1‑251 (MQ=255)
                     accacGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATACCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTAcc                                                                                                                                                                                                                           >  2:11550/1‑251 (MQ=255)
                         cgccaacccccccccccACCCCGCCGCGGAACGATCCCCCGCCTCCGACGCACGCGGCAATAGCATCCAGCGCATACATAAAGCCGAGGTTTTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACTCGCCAGAAAACGCGTAAATCAGCAGCAGGTTCAGGCCTACGTTGACACCAGATTCTTTTGCCGCTGCCGTGGTACCGCCAAGGGCAAAAATGTTCTTACCGAAg                                                                                                                                                                                                                       <  1:376895/238‑1 (MQ=255)
                          gCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGc                                                                                                                                                                                                                      >  1:322080/1‑251 (MQ=255)
                            caaccccccccccccccccgccgCTGACCGATCCCCCGCCTACCCCGCACCCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATCCGTCCGGCTTCAAACACCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGACCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTTTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCgg                                                                                                                                                                                                                    <  1:170661/234‑1 (MQ=255)
                                 cccgaccccACCCCGCCGCTGAACGATCCCCCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGtttt                                                                                                                                                                                                               <  1:330336/246‑1 (MQ=255)
                                                     gAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACg                                                                                                                                                                                           >  1:405284/1‑250 (MQ=255)
                                                       ccGATACCCCGCCTCCCACGCCCGCCGCGATAGCACCCAGCTCATACATAAACCCGAGGTTGTGGGTGGCAGAGCCGATACGCCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGTTATTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAAc                                                                                                                                                                                        <  2:379582/250‑1 (MQ=255)
                                                                 gCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCa                                                                                                                                                                              >  1:218548/1‑250 (MQ=255)
                                                                     aCCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCTCCAAAGTCATAGAACTCTCCAGACAACTCGTAGATCATCAGCAGGTTCATTCCGACGTTGACACCAGATATTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACTTCACCGCAATc                                                                                                                                                                           <  1:383972/249‑1 (MQ=255)
                                                                         cgcacgcCGCGTACGCATCCAGCTCATAAATAAAGCCGGGGTTGTGGGTGGCAGAGCCGATACGTCCGGCTTCTAACACCCCGCAAAAGGCAGAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACCCCAGATCCTTTTGCCGCTTCCGTGTTCCCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAAcgcg                                                                                                                                                                     <  2:170442/251‑1 (MQ=255)
                                                                             cgccgcGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACCCCACGACGAACGCCACCGCAATCCACGCGt                                                                                                                                                                    >  1:211712/1‑248 (MQ=255)
                                                                             cgccgcGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGACAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGt                                                                                                                                                                    >  2:88124/1‑248 (MQ=255)
                                                                                   gATAGCATCCAGCTCATACATAAACCCGAGGTTGTTGGGGGCAGAGCCGATACTTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACCCGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCATATACTTTGGCCCCTTCCGTTTTTCCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCTTGACGACGTTAAGCGCAATCAACGCGTAGAAGGTGa                                                                                                                                                           <  1:24065/251‑1 (MQ=255)
                                                                                         cacccAGCTCATACAAAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTGGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTaa                                                                                                                                                      <  1:404721/248‑1 (MQ=255)
                                                                                         cATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGGTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCACcgcgaagcaggcgaggta                                                                                                                                                       >  1:115569/1‑236 (MQ=255)
                                                                                             cAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGGGGAGAAGGGGATGTAag                                                                                                                                                     >  2:222033/1‑248 (MQ=255)
                                                                                               gCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCACGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGGGATGTAagcg                                                                                                                                                   >  2:367484/1‑246 (MQ=255)
                                                                                                 aaaaacaaaaaGCCGAGTTTGTTGGGGACAGCGCCGATACGCCCGGCTACTAACACCCCGCCAAAGGCATCAACCACGCCAGAAAACGCGAAGATCATCAGCAGGTTCAGGCCGACTGTTACACCAGATACTTTTGCCGCTTCCGTGGTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTGTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACg                                                                                                                                              <  2:62046/243‑1 (MQ=255)
                                                                                                            cccgatttggtgggtggCAGAGCCGATACGCCCGGCTTCTAACACCCCGCCAAAGGCATAGACCACGCCAGAAAACGCGTAGACCATCAGCAGGTTCAGGCCGACGTTGACCCCAGATCCTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTccc                                                                                                                                     <  2:82313/237‑1 (MQ=255)
                                                                                                              cGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAgc                                                                                                                                 >  1:93490/1‑251 (MQ=255)
                                                                                                                  gttgttGGTGGCAGAGCCGATCCGTCCGGCTTCTAACATCCCTCCAAAGGCATAGAACCCGCCAGACAACGCGTAGATCATCAGCGGGTTCAGGCCGCCGTTGACTCCAGATCCTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGa                                                                                                                             <  1:85276/251‑1 (MQ=255)
                                                                                                                    tgttgGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTccc                                                                                                                                     >  1:82313/1‑241 (MQ=255)
                                                                                                                    tgtgggtggCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCTTCAGCAGTATCAGGCCGCCGTTGACACCAGATACTTTTGCCGCTTCCTTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGAc                                                                                                                            <  1:426923/246‑1 (MQ=255)
                                                                                                                                    gATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAgg                                                                                                           >  1:236169/1‑251 (MQ=255)
                                                                                                                                                 tCTAACATCCCGCCAAAGGCATAGCACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGGTCTTACCGAAGCGGGTTTTGTTCCACAACACCAAGACGAACGCCACCGCAATCAACGCGTAGAAGGCGATGTAAGAGAGACGGAAACTCCCCACCGCGACAAAGCCCTGCGCAAAGGGAAAGAAGCCACt                                                                                              >  2:307268/1‑251 (MQ=255)
                                                                                                                                                                               ccAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGACCGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGAGCGACAAAGCCCTGAGCAACGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGAc                                                                >  2:193465/1‑251 (MQ=255)
                                                                                                                                                                                 aGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACa                                                               >  1:261348/1‑250 (MQ=255)
                                                                                                                                                                                    cAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCAGCCAAAgg                                                           >  2:404996/1‑250 (MQ=255)
                                                                                                                                                                                      aCGCGTAGATCATCGGCAGGTTCAGGCCGACGTTGACACCAGATTCTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCg                                                         <  1:193465/251‑1 (MQ=255)
                                                                                                                                                                                         ggtagacaatcaccggtttcggtcagACGTGGACACAAGTTCCTTTTGCCGCTTCCTTGTTCCCGCCAATGGCAAAAATGTTCTTGCCGAAGCGGGTTTTGTTCCACAACACCCAGACGAAAGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGTCCCTACAAAGTCGTaa                                                      <  2:249922/226‑1 (MQ=255)
                                                                                                                                                                                         cgTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGCCAAAGTCGTa                                                       >  2:320587/1‑250 (MQ=255)
                                                                                                                                                                                                  ccagcagGTTCAGGCCGCTGTTGACACCAGATCCTTTTCCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAg                                             <  2:132638/250‑1 (MQ=255)
                                                                                                                                                                                                        ggTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATg                                       >  1:328883/1‑251 (MQ=255)
                                                                                                                                                                                                             aGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCAt                                   >  1:89984/1‑250 (MQ=255)
                                                                                                                                                                                                               gCCGACTTTGACACCAGATACTTTTGCCGCTCCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAAGTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGa                                <  1:367484/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGACGTTGACACCAGATACTTTTGCCGCTTCCGTGTTACCGCCAATGGCCCAAATGTTCTTACCGAAGCGGGTTTTGTTCCACCACACCCAGACGCACGCCACCGCAATCAACGCGTCGAAGGTGATGTAAGAGCGACGGACACTCCCCAGCGCGACCAAGCCCTGAGCAAAGGTAGAAAAGCCACTGTCCAAACCAGAAATTGGCGACGCCCCGACAAAGTAGTAATAGAGCGCGTTGATGACATAGACg                              >  1:284688/1‑251 (MQ=255)
                                                                                                                                                                                                                       ttGACGCCAGATCCTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACgatg                           <  1:77462/249‑1 (MQ=255)
                                                                                                                                                                                                                                 aTACTTTTGCCGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCg               >  1:14374/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aCTTATGCCGCTTCCGGGTAACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAg             <  2:38651/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGCTTCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGt       >  2:31904/1‑250 (MQ=255)
                                                                                                                                                                                                                                              ttCCGTGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTACACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCAACTGTCAAACCCAGAAATTGGCGACGCCCCGACAACGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATg   >  1:68490/1‑250 (MQ=255)
                                                                                                                                                                                                                                               tCCGTGTTACCGCCAATGGCAAAAATGTGCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCTAGCGTGGTAATGa  <  2:61880/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
GTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGCCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGA  >  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: