Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,131,200 T→A 5.0% S61R (AGT→AGA) ‡ flgB → flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,131,2000TA 5.0% 794.0 / 21.4 218S61R (AGT→AGA) ‡flgBflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base T (124/83);  new base A (6/5);  total (130/88)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 228 total aligned reads displayed.

ACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGA  >  NC_000913/1131053‑1131349
                                                                                                                                                   |                                                                                                                                                     
acaaGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAgtgt                                                                                                                                                     >  2:2896784/1‑150 (MQ=255)
   agagGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGt                                                                                                                                                  <  2:2628760/150‑1 (MQ=255)
      ggCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGc                                                                                                                                               >  2:2189457/1‑150 (MQ=255)
      ggCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGc                                                                                                                                               >  2:2522124/1‑150 (MQ=255)
       gcgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCa                                                                                                                                              >  1:3602799/1‑150 (MQ=255)
        cgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTTGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCACACTGGTTGCat                                                                                                                                             <  2:3639402/150‑2 (MQ=255)
        cgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCACACTGGTTGCat                                                                                                                                             <  1:3706324/150‑2 (MQ=255)
            cAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGAc                                                                                                                                         <  2:3922936/150‑1 (MQ=255)
            cAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGAc                                                                                                                                         <  1:168297/150‑1 (MQ=255)
                 tGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGa                                                                                                                                    >  1:1015278/1‑150 (MQ=255)
                 tGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGa                                                                                                                                    <  1:664995/150‑1 (MQ=255)
                   cgcgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGAcc                                                                                                                                  <  1:694853/150‑1 (MQ=255)
                     cgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTc                                                                                                                                >  2:611427/1‑150 (MQ=255)
                       cccAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTCTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCaa                                                                                                                              <  1:2036565/150‑1 (MQ=255)
                        ccAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCACCGTGGACGGGATGCACCCAGTGTGGTTGCACTGACGATGACCTCacc                                                                                                                             >  1:2408464/1‑148 (MQ=255)
                        ccAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAAc                                                                                                                             >  1:1178503/1‑150 (MQ=255)
                        ccAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAAc                                                                                                                             >  1:2408587/1‑150 (MQ=255)
                        ccAGCGACAGGAAGTGCTGGCAGCAAACAACGCCAATGCCGATACCCCAGGTTAACAGGAGCGCGATATCGATATTGCCAGAGAACTTAAAAAAGTCATGCAACGAGGACGGGATGCAAACAGTGTGGTTGAACTGACCATGACCTCAAc                                                                                                                             >  2:1659833/1‑150 (MQ=255)
                         cAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACg                                                                                                                            >  2:889912/1‑150 (MQ=255)
                         cAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGATATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAAcc                                                                                                                            >  1:3900611/1‑149 (MQ=255)
                          aGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGc                                                                                                                           >  2:1744777/1‑150 (MQ=255)
                           gCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCa                                                                                                                          >  2:2055238/1‑150 (MQ=255)
                             gTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAac                                                                                                                        >  1:968870/1‑150 (MQ=255)
                              tCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAaca                                                                                                                       >  1:3203278/1‑150 (MQ=255)
                                 ggAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACAt                                                                                                                    <  1:3674495/150‑1 (MQ=255)
                                   aaGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCTACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTc                                                                                                                  <  1:2616779/150‑1 (MQ=255)
                                     gTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCg                                                                                                                <  1:929749/150‑1 (MQ=255)
                                         tGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGgcgc                                                                                                            >  1:1057210/1‑150 (MQ=255)
                                          ggCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCa                                                                                                           >  2:1374775/1‑150 (MQ=255)
                                               caAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgcgc                                                                                                      >  2:2342095/1‑150 (MQ=255)
                                               caAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgcgc                                                                                                      <  1:476135/150‑1 (MQ=255)
                                                aaaCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCt                                                                                                     <  1:3384929/150‑1 (MQ=255)
                                                 aaCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTg                                                                                                    >  1:3031098/1‑150 (MQ=255)
                                                 aaCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTg                                                                                                    >  2:495960/1‑150 (MQ=255)
                                                  aCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGa                                                                                                   <  1:2542090/150‑1 (MQ=255)
                                                    aTCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACg                                                                                                 >  1:3279202/1‑150 (MQ=255)
                                                     tCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGc                                                                                                >  2:1159755/1‑150 (MQ=255)
                                                     tCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACGTAAAAAAGTCATGCACCGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGc                                                                                                >  2:544950/1‑150 (MQ=255)
                                                        ccAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctc                                                                                             <  2:1952249/150‑1 (MQ=255)
                                                         cAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctcc                                                                                            <  1:179637/150‑1 (MQ=255)
                                                         cAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctcc                                                                                            <  2:3654588/150‑1 (MQ=255)
                                                         cAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctcc                                                                                            <  1:3849058/150‑1 (MQ=255)
                                                            tGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTAc                                                                                         <  1:1248482/150‑1 (MQ=255)
                                                             gCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTAcc                                                                                        <  1:2118889/150‑1 (MQ=255)
                                                                     cccTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACt                                                                                >  2:645080/1‑150 (MQ=255)
                                                                           ttATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCGTGCAACGTGGACGGGATGCATCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATa                                                                          <  2:3469751/150‑1 (MQ=255)
                                                                             aTCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATAcc                                                                        >  2:1033198/1‑150 (MQ=255)
                                                                                   cgcgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTc                                                                  >  1:3657875/1‑150 (MQ=255)
                                                                                    gcgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTcc                                                                 >  2:3621271/1‑150 (MQ=255)
                                                                                     cgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCg                                                                <  2:2866715/150‑1 (MQ=255)
                                                                                     cgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCg                                                                <  2:786524/150‑1 (MQ=255)
                                                                                       cgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGa                                                              >  1:577957/1‑150 (MQ=255)
                                                                                       cgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGa                                                              >  2:2368359/1‑150 (MQ=255)
                                                                                       cgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGa                                                              >  1:305150/1‑150 (MQ=255)
                                                                                        gagatCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTCCCGCAGAACTGCAATACCGTATTCCGGAc                                                             >  1:2858474/4‑150 (MQ=255)
                                                                                        gATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGAc                                                             >  1:3879267/1‑150 (MQ=255)
                                                                                           atCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAg                                                          >  1:1462740/1‑150 (MQ=255)
                                                                                           atCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAg                                                          >  1:1314955/1‑150 (MQ=255)
                                                                                            tCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGc                                                         <  2:2702875/150‑1 (MQ=255)
                                                                                            tCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGc                                                         <  1:2246605/150‑1 (MQ=255)
                                                                                            tCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGc                                                         >  1:1844896/1‑150 (MQ=255)
                                                                                              gATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCt                                                       >  2:2777825/1‑150 (MQ=255)
                                                                                              gATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCt                                                       >  2:1903650/1‑150 (MQ=255)
                                                                                              gATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCt                                                       >  2:3195798/1‑150 (MQ=255)
                                                                                                      cAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGa                                               <  2:105812/150‑1 (MQ=255)
                                                                                                      cAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGa                                               <  1:663313/150‑1 (MQ=255)
                                                                                                         tGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCACTACCGTAGTCCGGACCAGCCTTCGCTTGACgg                                            >  2:441000/1‑150 (MQ=255)
                                                                                                         tGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACgg                                            >  1:516223/1‑150 (MQ=255)
                                                                                                           aaCTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTa                                          <  1:2473015/150‑1 (MQ=255)
                                                                                                              ttAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATa                                       >  1:2014854/1‑150 (MQ=255)
                                                                                                              ttAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATa                                       >  2:2110383/1‑150 (MQ=255)
                                                                                                               tAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAc                                      >  1:1724438/1‑150 (MQ=255)
                                                                                                                aaaaaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAcc                                     <  2:146466/150‑1 (MQ=255)
                                                                                                                   aaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTc                                  <  1:1346339/150‑1 (MQ=255)
                                                                                                                             aaCGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATc                        <  2:894745/150‑1 (MQ=255)
                                                                                                                               cGTGGACGGGATGCAACCAGTGTGGGTGCACTGACGAGGACCTCCACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTCCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTACTACCGTCGATATGGATcgc                      <  2:2197565/150‑1 (MQ=255)
                                                                                                                                 tGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGa                    >  2:381521/1‑150 (MQ=255)
                                                                                                                                 tGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGa                    >  2:102827/1‑150 (MQ=255)
                                                                                                                                  ggACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGaa                   <  1:1543745/150‑1 (MQ=255)
                                                                                                                                  ggACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGaa                   <  1:1297385/150‑1 (MQ=255)
                                                                                                                                  ggACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGaa                   <  2:1250125/150‑1 (MQ=255)
                                                                                                                                  ggACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGaa                   <  2:325403/150‑1 (MQ=255)
                                                                                                                                   gACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAAc                  >  1:282045/1‑150 (MQ=255)
                                                                                                                                    aCGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACg                 >  1:432688/1‑150 (MQ=255)
                                                                                                                                     gtcagtGCAACCACACTGGTTGCACTGACGATGACCTCAACGCAACTCATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTCTTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATAGGGCTCGCGAACGc                >  2:2424181/6‑150 (MQ=255)
                                                                                                                                     gtcagtGCAACCACACTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTCATACCGTCGATATGGATCGCGAACGc                >  2:1465820/6‑150 (MQ=255)
                                                                                                                                      gggATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCa               >  2:3420561/1‑150 (MQ=255)
                                                                                                                                      gggATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCa               >  2:408769/1‑150 (MQ=255)
                                                                                                                                      gggATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCa               <  1:2582925/150‑1 (MQ=255)
                                                                                                                                         gtGCAACCACACTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGACCGCAccc            >  1:3374604/2‑150 (MQ=255)
                                                                                                                                         aTGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAccc            >  2:313148/1‑150 (MQ=255)
                                                                                                                                          tGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCa           <  2:3512263/150‑1 (MQ=255)
                                                                                                                                           gCAACCAGTGTGGTTGCACTGACGACGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAg          <  2:3077330/150‑1 (MQ=255)
                                                                                                                                             aaCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGtt        >  1:1604691/1‑150 (MQ=255)
                                                                                                                                             aaCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGtt        >  1:2678058/1‑150 (MQ=255)
                                                                                                                                              aCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGttt       >  1:2943038/1‑150 (MQ=255)
                                                                                                                                                cAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGc     <  2:2947920/150‑1 (MQ=255)
                                                                                                                                                  gtgtGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCg   >  1:3469674/1‑150 (MQ=255)
                                                                                                                                                   tgtgGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGa  <  2:2868460/150‑1 (MQ=255)
                                                                                                                                                   tgtgGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGa  >  1:2742846/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGA  >  NC_000913/1131053‑1131349

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