Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,139,790 T→A 5.3% A471A (GCT→GCA flgK → flagellar hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,139,7900TA 5.3% 390.8 / 9.7 113A471A (GCT→GCAflgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base T (62/45);  new base A (4/2);  total (66/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

Only 100 of 113 total aligned reads displayed.

ATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCA  >  NC_000913/1139644‑1139936
                                                                                                                                                  |                                                                                                                                                  
aTGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTa                                                                                                                                                  <  4:851036/149‑1 (MQ=255)
aTGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTa                                                                                                                                                  >  3:427205/1‑149 (MQ=255)
     tGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGcc                                                                                                                                              >  2:473781/1‑148 (MQ=255)
       tATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCg                                                                                                                                           >  4:607889/1‑149 (MQ=255)
        aTTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGt                                                                                                                                          <  4:810892/149‑1 (MQ=255)
          tAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGtt                                                                                                                                         <  2:561948/148‑1 (MQ=255)
            aTCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagt                                                                                                                                      <  4:576278/149‑1 (MQ=255)
             tCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagt                                                                                                                                      >  2:193915/1‑148 (MQ=255)
              cACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtga                                                                                                                                    >  3:872323/1‑149 (MQ=255)
               aCCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtgag                                                                                                                                   >  3:548536/1‑149 (MQ=255)
               aCCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtga                                                                                                                                    >  2:21901/1‑148 (MQ=255)
                ccGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtgagt                                                                                                                                  >  4:605915/1‑149 (MQ=255)
                ccGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtgagt                                                                                                                                  <  3:880314/149‑1 (MQ=255)
                ccGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtgag                                                                                                                                   >  1:949933/1‑148 (MQ=255)
                ccGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtgag                                                                                                                                   <  1:122500/148‑1 (MQ=255)
                 cgacgaAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtg                                                                                                                                     <  1:354307/145‑1 (MQ=255)
                 cgacgaAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtg                                                                                                                                     >  2:354307/1‑145 (MQ=255)
                  gacgaAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAgtgagtga                                                                                                                                <  4:520741/149‑1 (MQ=255)
                  gacgaAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGAGCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAgtgagtga                                                                                                                                >  3:649889/1‑149 (MQ=255)
                    cgaAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGat                                                                                                                               <  1:210063/148‑1 (MQ=255)
                     gaAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGatat                                                                                                                             <  3:843260/149‑1 (MQ=255)
                         cGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCgg                                                                                                                          >  1:745588/1‑148 (MQ=255)
                         cGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGt                                                                                                                         <  3:225746/149‑1 (MQ=255)
                          gAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTCGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTTGTGAGTGATATCGGta                                                                                                                        >  4:775785/1‑149 (MQ=255)
                           aaaaTAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtaa                                                                                                                       >  4:241599/1‑149 (MQ=255)
                             aaTAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtaata                                                                                                                     >  3:718171/1‑149 (MQ=255)
                                aGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATaaaa                                                                                                                  >  4:438971/1‑149 (MQ=255)
                                  cGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAAc                                                                                                                 <  2:122676/148‑1 (MQ=255)
                                    aTGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTag                                                                                                                                       >  2:834324/1‑124 (MQ=255)
                                    aTGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTCag                                                                                                                                       <  1:834324/124‑1 (MQ=255)
                                          agcgaAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGt                                                                                                         >  2:398578/1‑148 (MQ=255)
                                            cgaAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGa                                                                                                      <  3:654030/149‑1 (MQ=255)
                                              aagaagATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaa                                                                                                     <  1:232183/148‑1 (MQ=255)
                                               agaagaTGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaaaa                                                                                                   <  4:83279/149‑1 (MQ=255)
                                               agaagaTGCGCGTGAAAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCGAAACAACAGTAAAACGGTGGGCGGCGCGAAATCCTTTCACGACGCTTATGCCTCGTTAGTGTGTGATATCGGTAATAAAACCGCGACGTTGacaa                                                                                                   >  4:705778/1‑149 (MQ=255)
                                                  agaTGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATAACGGTAATAAAACCGCGACGTTGAAAACCa                                                                                                <  3:775785/149‑1 (MQ=255)
                                                          gTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGcc                                                                                        <  3:607889/149‑1 (MQ=255)
                                                            gATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCa                                                                                       <  1:381568/148‑1 (MQ=255)
                                                               aGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGCGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCa                                                                                   >  4:305632/1‑149 (MQ=255)
                                                                   aTAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAgg                                                                                >  1:434493/1‑148 (MQ=255)
                                                                   aTAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAgg                                                                                >  1:139944/1‑148 (MQ=255)
                                                                    tAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGt                                                                               >  1:227204/1‑148 (MQ=255)
                                                                    tAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTa                                                                              >  3:95758/1‑149 (MQ=255)
                                                                      aCCGCAACGGTCAGGCCATGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGTGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGGGTGAGACCGGTAATAATACCGCGACGTTGAAAACCAGTAGCGCTACGCAAGGTAAt                                                                            <  4:743030/149‑1 (MQ=255)
                                                                       ccGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                           >  3:188797/1‑149 (MQ=255)
                                                                       ccGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAt                                                                            <  1:657265/148‑1 (MQ=255)
                                                                        cGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                           <  2:776977/148‑1 (MQ=255)
                                                                                 cAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTACTAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCa                                                                 >  3:96364/1‑149 (MQ=255)
                                                                                    gCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGc                                                               <  2:949933/148‑1 (MQ=255)
                                                                                    gCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCt                                                              >  4:346449/1‑149 (MQ=255)
                                                                                     cccTGCTGGATCTGCAAATCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTACGTCTCGTTagtg                                                                                                                                     <  1:790838/77‑1 (MQ=255)
                                                                                     cccTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTACGCCTCGTTagtg                                                                                                                                     >  2:790838/1‑77 (MQ=255)
                                                                                       ctgctgGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCttt                                                            >  2:944815/1‑148 (MQ=255)
                                                                                        tgctgGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGta                                                                                                                        <  3:493487/87‑1 (MQ=255)
                                                                                        tgctgGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGta                                                                                                                        >  4:493487/1‑87 (MQ=255)
                                                                                        tgctgGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTc                                                           <  1:21901/148‑1 (MQ=255)
                                                                                           tgGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAAt                                                       >  3:308302/1‑149 (MQ=255)
                                                                                            ggATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATc                                                      >  4:870005/1‑149 (MQ=255)
                                                                                                cTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagca                                                  <  3:872590/149‑1 (MQ=255)
                                                                                                 tGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagcag                                                 >  4:881146/1‑149 (MQ=255)
                                                                                                  gCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagcag                                                 >  2:595698/1‑148 (MQ=255)
                                                                                                   cAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagcagc                                                >  2:19106/1‑148 (MQ=255)
                                                                                                    aaaGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagcagcag                                              >  3:939969/1‑149 (MQ=255)
                                                                                                     aaGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGt                                             <  3:912522/149‑1 (MQ=255)
                                                                                                       gCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCg                                           >  4:534575/1‑149 (MQ=255)
                                                                                                       gCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGAAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTc                                            >  2:586075/1‑148 (MQ=255)
                                                                                                        cAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGTGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGTAGCTTTCCAATCAGCAGCAGTCGa                                          >  3:292587/1‑149 (MQ=255)
                                                                                                        cAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCg                                           <  2:87347/148‑1 (MQ=255)
                                                                                                           cAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGAtt                                        >  2:22657/1‑148 (MQ=255)
                                                                                                             gTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTc                                      >  2:805308/1‑148 (MQ=255)
                                                                                                               aaaaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAAcc                                                                                                 <  2:594090/87‑1 (MQ=255)
                                                                                                                aaaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCgg                                   >  1:489008/1‑148 (MQ=255)
                                                                                                                aaaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCgg                                   <  1:811748/148‑1 (MQ=255)
                                                                                                                aaaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGgt                                  >  3:518168/1‑149 (MQ=255)
                                                                                                                aaaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGgt                                  <  4:718171/149‑1 (MQ=255)
                                                                                                                  aCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGt                                                                                                         >  1:883951/1‑76 (MQ=255)
                                                                                                                   cGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTc                               <  3:534575/149‑1 (MQ=255)
                                                                                                                     gTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTcc                                                          <  2:471423/120‑1 (MQ=255)
                                                                                                                     gTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTcc                                                          >  1:471423/1‑120 (MQ=255)
                                                                                                                      tGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGTTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAAt                            <  4:603492/149‑1 (MQ=255)
                                                                                                                       gggCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAAtc                           >  4:638961/1‑149 (MQ=255)
                                                                                                                         gCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAAtct                          <  1:856269/148‑1 (MQ=255)
                                                                                                                            gTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGa                       <  2:772303/148‑1 (MQ=255)
                                                                                                                              gCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATga                    <  4:427205/149‑1 (MQ=255)
                                                                                                                               cGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATgag                   <  4:308302/149‑1 (MQ=255)
                                                                                                                                gAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATgagg                  >  4:769266/1‑149 (MQ=255)
                                                                                                                                   aTCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATgaggag                >  1:469679/1‑148 (MQ=255)
                                                                                                                                    tCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGt               >  2:852540/1‑148 (MQ=255)
                                                                                                                                      cactAACGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTAc             >  1:470754/4‑148 (MQ=255)
                                                                                                                                       tttAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCATGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACgg           >  3:596675/1‑149 (MQ=255)
                                                                                                                                       tttAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACgg           >  4:401001/1‑149 (MQ=255)
                                                                                                                                        ttAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACgg           >  1:233199/1‑148 (MQ=255)
                                                                                                                                        ttAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACgg           >  1:109221/1‑148 (MQ=255)
                                                                                                                                         tAACGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATga                    <  3:577390/138‑1 (MQ=255)
                                                                                                                                          aaCGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCgg                                   >  1:754516/1‑122 (MQ=255)
                                                                                                                                          aaCGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCgg                                   <  2:754516/122‑1 (MQ=255)
                                                                                                                                          aaCGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGgtg                                 <  2:556032/124‑1 (MQ=255)
                                                                                                                                          aaCGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGgtg                                 >  1:556032/1‑124 (MQ=255)
                                                                                                                                           acgacgCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAAt       >  4:705159/1‑149 (MQ=255)
                                                                                                                                             gacgCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCt     >  3:258576/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCA  >  NC_000913/1139644‑1139936

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