Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 731,964 G→A W794* (TGG→TGA rhsC → rhs element protein RhsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913731,9640GA98.3% 194.5 / NA 59W794* (TGG→TGArhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/31);  minor base C (1/0);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGG  >  NC_000913/731677‑732234
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                              
tacgACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATTAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGc                                                                                                                                                                                                                                                                    <  2:398576/297‑1 (MQ=255)
tATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCTAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGc                                                                                                                                                                                                                                                                    <  2:138868/300‑1 (MQ=255)
tATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACATTAAAGTCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGc                                                                                                                                                                                                                                                                    <  2:332797/300‑1 (MQ=255)
            gACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGTTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGgctggctg                                                                                                                                                                                                                                                        <  1:348642/300‑1 (MQ=255)
              ccGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTc                                                                                                                                                                                                                                                      >  1:21283/1‑300 (MQ=255)
                gCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGAAAGGCATCAGCCATATCAGCGAAGAGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTAGCCAATGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGACGAACCGCTGTATACCGAACAGCCTGACCGCAGTGGAATGACTGAACTACGAAAGCGGCTGGATGTCAg                                                                                                                                                                                                                                                    >  2:770575/1‑300 (MQ=21)
                     aCGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGACGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATg                                                                                                                                                                                                                                               >  1:120910/1‑300 (MQ=255)
                            gcacagTGAATGGCGAAACCGCCGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGATGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACAACTGACGGTGCATCATCCGCCGACGAATAAACTGCCCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATAACGAACAGCCTGCCCGACGTGGAATGCCTGACCTACGACAGCGGCTGGTCGTCAGGCATGAAACTCg                                                                                                                                                                                                                                        >  2:346117/1‑300 (MQ=21)
                                              ccGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGNGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCtggtgg                                                                                                                                                                                                                      >  2:742740/1‑300 (MQ=255)
                                              ccGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCtggtgg                                                                                                                                                                                                                      >  1:474423/1‑300 (MQ=255)
                                                gCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGTTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAg                                                                                                                                                                                                                    <  2:768431/300‑1 (MQ=255)
                                                   gagCGGTTGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTac                                                                                                                                                                                                                 <  1:819511/300‑1 (MQ=255)
                                                   gagCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTac                                                                                                                                                                                                                 <  2:315268/300‑1 (MQ=255)
                                                   gagCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCAGATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTac                                                                                                                                                                                                                 <  2:802777/300‑1 (MQ=255)
                                                      cGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGTGCACCGGTTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACAcc                                                                                                                                                                                                              <  2:649191/300‑1 (MQ=255)
                                                        gTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCcg                                                                                                                                                                                                            <  1:470460/300‑1 (MQ=255)
                                                            cAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCACGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAGACTCGGCGACACACCGCTGGTGGAGTACACCCGCGAc                                                                                                                                                                                                        >  2:43654/1‑300 (MQ=255)
                                                                     gAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCNTCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCAc                                                                                                                                                                                               <  1:608179/300‑1 (MQ=255)
                                                                               ggctgACAGACATCCGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGc                                                                                                                                                                                     >  1:226275/1‑300 (MQ=255)
                                                                               ggctgACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGc                                                                                                                                                                                     <  1:640665/300‑1 (MQ=255)
                                                                               ggctgACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAAACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCAGTGGAGTGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGc                                                                                                                                                                                     >  1:226381/1‑300 (MQ=255)
                                                                                         caTCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGNAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCtt                                                                                                                                                                           <  2:640115/300‑1 (MQ=255)
                                                                                                  tatCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTa                                                                                                                                                                  >  1:656528/1‑300 (MQ=255)
                                                                                                             ggCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTcacca                                                                                                                                                       >  1:265226/1‑300 (MQ=255)
                                                                                                                   gggTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGNAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCCCCGCtt                                                                                                                                                 >  2:203020/1‑300 (MQ=255)
                                                                                                                   gggTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCtt                                                                                                                                                 >  1:797207/1‑300 (MQ=255)
                                                                                                                        aCGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCCGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATAcc                                                                                                                                            >  2:550577/1‑300 (MQ=255)
                                                                                                                        aCGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTtatccc                                                                                                                                            >  1:423513/1‑297 (MQ=255)
                                                                                                                        aCGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATAcc                                                                                                                                            <  1:701125/300‑1 (MQ=255)
                                                                                                                        aCGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCAGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTCTGAACTCACCACCGCTTATAcc                                                                                                                                            >  1:472060/1‑300 (MQ=255)
                                                                                                                                 tATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCggg                                                                                                                                   <  1:459395/300‑1 (MQ=255)
                                                                                                                                    ggTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAg                                                                                                                                >  1:823901/1‑300 (MQ=255)
                                                                                                                                    ggTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCAGGGCAg                                                                                                                                >  2:537827/1‑300 (MQ=255)
                                                                                                                                            cAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACagag                                                                                                                        <  1:707581/300‑1 (MQ=255)
                                                                                                                                            cAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACagag                                                                                                                        <  2:332961/300‑1 (MQ=255)
                                                                                                                                            cAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTACTCTGGCAGCATGAGACCATACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACagag                                                                                                                        <  2:12585/300‑1 (MQ=255)
                                                                                                                                                         cTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACNGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAAc                                                                                                           <  1:609827/300‑1 (MQ=255)
                                                                                                                                                                gTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCctgc                                                                                                    <  1:677469/300‑1 (MQ=255)
                                                                                                                                                                  gaacaccttaatgTGCATCATCCGCATACGAATGAACTGCTCTGTCAGCATTAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGTTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACGCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCctgctg                                                                                                  <  2:388561/288‑1 (MQ=14)
                                                                                                                                                                       ccTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTACGAAGCTTCGGCCGTTATGAACTCACCAACGCCTCTACCCATGCCGAGCAGTTAACGAGCCCGCACCTGAACAGCCTGCTGTCTGa                                                                                             >  2:826414/1‑300 (MQ=14)
                                                                                                                                                                          gACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACAGCTGGTGAAGTACACCCGCGAACGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTAACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCttgccg                                                                                          >  2:102012/1‑295 (MQ=255)
                                                                                                                                                                              gTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGAt                                                                                      >  1:118678/1‑300 (MQ=255)
                                                                                                                                                                                  atcatcCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACNCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTNCCGGGCAGTTACAGAGCCAGCACCTGCACAGCCTGCTGTCTGACCGCGATTaca                                                                                  >  2:793209/1‑300 (MQ=255)
                                                                                                                                                                                   tcatcCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTacac                                                                                 <  2:37316/300‑1 (MQ=255)
                                                                                                                                                                                   gcatcCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTacac                                                                                 <  1:187539/299‑1 (MQ=255)
                                                                                                                                                                                             aCGAATGGACTGCTCTGGCAGCATTAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTTGCTGTCAGGCATGAAACTCGGCGACACACNGCTGGTGGAGTACACCCGCGACCGCCNGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGAc                                                                       <  2:568918/300‑1 (MQ=255)
                                                                                                                                                                                                   gAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCACGATTACACCTGGGACGACAATGGc                                                                 >  2:325542/1‑300 (MQ=255)
                                                                                                                                                                                                                  caTGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCANGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGTAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACTGCGAACTCATCCGCATc                                                  <  1:617138/300‑1 (MQ=255)
                                                                                                                                                                                                                       gaCCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATcagcag                                             >  1:232549/1‑300 (MQ=39)
                                                                                                                                                                                                                         ccAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGcc                                           <  1:712411/300‑1 (MQ=39)
                                                                                                                                                                                                                              cATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCgcgc                                      >  1:184287/1‑300 (MQ=39)
                                                                                                                                                                                                                                tGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGAACGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCa                                    >  1:738399/1‑300 (MQ=255)
                                                                                                                                                                                                                                    tACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGAcc                                <  1:223433/300‑1 (MQ=39)
                                                                                                                                                                                                                                                 taCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTtacagct                   <  2:340219/299‑1 (MQ=25)
                                                                                                                                                                                                                                                 gACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGGGGAGTACACCCGCGACCGCCTGCACCGGGNAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTtacagct                   <  2:196949/300‑1 (MQ=37)
                                                                                                                                                                                                                                                        gAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCTGAGTTACAGCTACAGcac            <  1:569173/300‑1 (MQ=25)
                                                                                                                                                                                                                                                             gCTGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCg       >  1:185060/1‑300 (MQ=25)
                                                                                                                                                                                                                                                               tGTATACCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGc     <  2:540389/300‑1 (MQ=25)
                                                                                                                                                                                                                                                                  ataCCGGACAGCCTGCCCGCCGTGGAATGACTGACCTACGGAAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTCCACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGAACGCGATTACAACTGGAACGACAACGGCGAACCCATCCGCATCAGCAGCCAACGCCAGACCCGGAGTTACAGCTACCGCACCACCGGCAgg  >  2:509621/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                              
TATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATACCCCTGCCGGGCAGTTACAGAGCCAGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGG  >  NC_000913/731677‑732234

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