Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,131,199 G→C 5.0% S61T (AGT→ACT) ‡ 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,1990GC 5.0% 811.0 / 27.2 218S61T (AGT→ACT) ‡flgBflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base G (124/83);  new base C (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 = 8.96e-01

Only 100 of 224 total aligned reads displayed.

ACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCG  >  NC_000913/1131053‑1131348
                                                                                                                                                  |                                                                                                                                                     
acaaGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAgtgt                                                                                                                                                    >  2:2896782/1‑150 (MQ=255)
   agagGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGt                                                                                                                                                 <  2:117433/150‑1 (MQ=255)
      ggCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGc                                                                                                                                              >  2:2522124/1‑150 (MQ=255)
       gcgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCa                                                                                                                                             >  1:2414349/1‑150 (MQ=255)
        cgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTTGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCACACTGGTTGCat                                                                                                                                            <  2:3639402/150‑2 (MQ=255)
        cgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCACACTGGTTGCat                                                                                                                                            <  1:3706324/150‑2 (MQ=255)
         gcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACt                                                                                                                                           >  2:2620919/1‑150 (MQ=255)
            cAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGAc                                                                                                                                        <  2:3922936/150‑1 (MQ=255)
            cAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGAc                                                                                                                                        <  1:168297/150‑1 (MQ=255)
                 tGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGa                                                                                                                                   >  1:1015278/1‑150 (MQ=255)
                     cgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTc                                                                                                                               >  2:611427/1‑150 (MQ=255)
                        ccAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAAc                                                                                                                            >  1:1178503/1‑150 (MQ=255)
                        ccAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAAc                                                                                                                            >  1:2408570/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:2378867/150‑1 (MQ=255)
                           gCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCa                                                                                                                         >  2:2055238/1‑150 (MQ=255)
                              tCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAaca                                                                                                                      >  1:3203278/1‑150 (MQ=255)
                                 ggAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAATAAAGTCATGCAACGTGGACGGGTTGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACAt                                                                                                                   <  1:3655001/150‑1 (MQ=255)
                                 ggAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACAt                                                                                                                   <  1:3674495/150‑1 (MQ=255)
                                        cTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGgcg                                                                                                            <  1:461630/150‑1 (MQ=255)
                                         tGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGgcgc                                                                                                           >  2:1046961/1‑150 (MQ=255)
                                         tGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGgcgc                                                                                                           >  1:1057210/1‑150 (MQ=255)
                                           gcagcaAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAg                                                                                                         >  1:3582222/1‑150 (MQ=255)
                                               caAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgcgc                                                                                                     <  2:2652660/150‑1 (MQ=255)
                                               caAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGgcgc                                                                                                     >  2:2993336/1‑150 (MQ=255)
                                                aaaCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCt                                                                                                    <  1:3384929/150‑1 (MQ=255)
                                                 aaCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTg                                                                                                   >  2:495960/1‑150 (MQ=255)
                                                  aCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGa                                                                                                  <  1:2542090/150‑1 (MQ=255)
                                                   cATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGAc                                                                                                 <  2:1451129/150‑1 (MQ=255)
                                                   cATCGCCAATGCCGATACCCCTCGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGACATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGAc                                                                                                 <  1:3803943/150‑1 (MQ=255)
                                                    aTCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACg                                                                                                >  1:3279202/1‑150 (MQ=255)
                                                     tCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGc                                                                                               <  1:2971571/150‑1 (MQ=255)
                                                     tCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACGTAAAAAAGTCATGCACCGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGc                                                                                               >  2:544950/1‑150 (MQ=255)
                                                        ccAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGcctc                                                                                            <  2:1952249/150‑1 (MQ=255)
                                                            tGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTAc                                                                                        <  1:1248482/150‑1 (MQ=255)
                                                            tGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTGGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTAc                                                                                        <  1:98629/150‑1 (MQ=255)
                                                            tGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTAAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTAc                                                                                        <  1:1716796/150‑1 (MQ=255)
                                                                     cccTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAggtat                                                                               <  2:665791/150‑5 (MQ=255)
                                                                     cccTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACt                                                                               >  2:3738265/1‑150 (MQ=255)
                                                                          gTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAAt                                                                          >  1:3873678/1‑150 (MQ=255)
                                                                          gTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAAAACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAAt                                                                          >  1:1526399/1‑150 (MQ=255)
                                                                                  gcgcgcGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTAtt                                                                  >  1:2674923/1‑150 (MQ=255)
                                                                                    gcgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTcc                                                                >  2:1712892/1‑150 (MQ=255)
                                                                                    gcgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTcc                                                                >  2:3621271/1‑150 (MQ=255)
                                                                                    gcgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTcc                                                                >  2:3499187/1‑150 (MQ=255)
                                                                                     cgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGACCTGCAATACCGTCTTCCg                                                               >  1:684738/1‑150 (MQ=255)
                                                                                     cgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCg                                                               <  2:786524/150‑1 (MQ=255)
                                                                                      gcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCgg                                                              >  1:2550999/1‑150 (MQ=255)
                                                                                       cgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGa                                                             >  1:305150/1‑150 (MQ=255)
                                                                                       cgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGa                                                             >  1:577957/1‑150 (MQ=255)
                                                                                       cgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGa                                                             >  2:2323875/1‑150 (MQ=255)
                                                                                        gagatCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTCCCGCAGAACTGCAATACCGTATTCCGGAc                                                            >  1:2858474/4‑150 (MQ=255)
                                                                                           atCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAg                                                         >  1:1314955/1‑150 (MQ=255)
                                                                                           atCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAg                                                         <  1:2062749/150‑1 (MQ=255)
                                                                                           atCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAg                                                         >  1:1462740/1‑150 (MQ=255)
                                                                                            tCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGc                                                        >  1:1844896/1‑150 (MQ=255)
                                                                                            tCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGc                                                        <  1:2246605/150‑1 (MQ=255)
                                                                                            tCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGc                                                        >  1:935626/1‑150 (MQ=255)
                                                                                                 tttGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCg                                                   <  1:2589749/150‑1 (MQ=255)
                                                                                                    gCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCtt                                                >  2:1753626/1‑150 (MQ=255)
                                                                                                      cAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGa                                              >  2:2118918/1‑150 (MQ=255)
                                                                                                         tGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCACTACCGTCTTCCGGACCAGCCTTCGCTTGACgg                                           >  2:441679/1‑150 (MQ=255)
                                                                                                         tGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACgg                                           >  2:2709692/1‑150 (MQ=255)
                                                                                                           aaCTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACTGTATTCCGGACCAGCCTTCGCTTGACGGTa                                         >  1:664689/1‑150 (MQ=255)
                                                                                                              ttAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATa                                      >  1:2014854/1‑150 (MQ=255)
                                                                                                               tAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAc                                     >  2:3343705/1‑150 (MQ=255)
                                                                                                                aaaaaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCACCGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAcc                                    <  2:2365139/150‑1 (MQ=255)
                                                                                                                   aaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTc                                 <  1:1346339/150‑1 (MQ=255)
                                                                                                                      gTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGat                              >  2:1092802/1‑150 (MQ=255)
                                                                                                                         aTGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATg                           >  2:321297/1‑150 (MQ=255)
                                                                                                                            cAACGTGGACGGGATGCAACCAGTGAGGTTGCACTGACGATGACCTCAACGCAACACATTCAGGCGCAGGCGCTGACGCCTCCAACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACTGTAATACCGTCGATATGGAt                        >  2:2074922/1‑150 (MQ=255)
                                                                                                                             aaCGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATc                       <  2:894745/150‑1 (MQ=255)
                                                                                                                             aaCGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATc                       <  2:3582222/150‑1 (MQ=255)
                                                                                                                               cGTGGACGGGATGCAACCAGTGTGGGTGCACTGACGAGGACCTCCACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTCCCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTACTACCGTCGATATGGATcgc                     <  2:2197565/150‑1 (MQ=255)
                                                                                                                                gTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATcgcg                    >  1:390518/1‑150 (MQ=255)
                                                                                                                                 tGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGa                   >  2:103052/1‑150 (MQ=255)
                                                                                                                                  ggACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGaa                  <  1:1297192/150‑1 (MQ=255)
                                                                                                                                  ggACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGaa                  <  2:1250351/150‑1 (MQ=255)
                                                                                                                                   gACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAAc                 <  1:1744777/150‑1 (MQ=255)
                                                                                                                                   gACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAAc                 >  1:282045/1‑150 (MQ=255)
                                                                                                                                   gACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACAGTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAAc                 <  1:478094/150‑1 (MQ=255)
                                                                                                                                    aCGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACg                >  1:202977/1‑150 (MQ=255)
                                                                                                                                     gtcagtGCAACCACACTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTCATACCGTCGATATGGATCGCGAACGc               >  2:1465820/6‑150 (MQ=255)
                                                                                                                                     cGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGc               >  1:2297069/1‑150 (MQ=255)
                                                                                                                                      gggATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCa              >  2:408769/1‑150 (MQ=255)
                                                                                                                                      gggATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCa              >  2:3420561/1‑150 (MQ=255)
                                                                                                                                        gATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCTGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAcc            >  2:2913801/1‑150 (MQ=255)
                                                                                                                                         gtGCAACCACACTGGTTGCACTGACGATGACCACAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGGCTGCCGGACCAGCCTTCGCTTCACGGTAATACCGTCGCTATGGATCGCGAACGCAccc           >  1:3374567/2‑150 (MQ=255)
                                                                                                                                         aTGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAccc           >  2:2846546/1‑150 (MQ=255)
                                                                                                                                          tGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCa          <  2:3880590/150‑1 (MQ=255)
                                                                                                                                           gCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAg         >  1:842608/1‑150 (MQ=255)
                                                                                                                                           gCAACCAGTGTGGTTGCACTGACGACGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAg         <  2:3077330/150‑1 (MQ=255)
                                                                                                                                             aaCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGtt       >  1:648623/1‑150 (MQ=255)
                                                                                                                                             aaCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGtt       >  1:1604691/1‑150 (MQ=255)
                                                                                                                                              aCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGttt      >  1:2943038/1‑150 (MQ=255)
                                                                                                                                              aCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGttt      >  1:2741910/1‑150 (MQ=255)
                                                                                                                                              aCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGttt      >  2:3163922/1‑150 (MQ=255)
                                                                                                                                                  gtgtGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCg  >  1:3469674/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCG  >  NC_000913/1131053‑1131348

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