Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,972,883 T→A 5.5% Q222L (CAG→CTG)  flgL ← flagellar hook‑associated protein FlgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,972,8830TA 5.5% 229.9 / 2.2 73Q222L (CAG→CTG) flgLflagellar hook‑associated protein FlgL
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base A (2/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTC  >  NC_002947/4972737‑4973026
                                                                                                                                                  |                                                                                                                                               
gCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGGCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGa                                                                                                                                               <  6:39930/149‑1 (MQ=255)
gCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGa                                                                                                                                               >  6:151451/1‑149 (MQ=255)
gCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGa                                                                                                                                               <  2:219601/149‑1 (MQ=255)
 cGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGAt                                                                                                                                              <  4:287447/149‑1 (MQ=255)
 cGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGAt                                                                                                                                              <  4:150444/149‑1 (MQ=255)
   agaTCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTg                                                                                                                                            >  5:316086/1‑149 (MQ=255)
    gaTCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGc                                                                                                                                           >  4:151271/1‑149 (MQ=255)
     atcatcTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGcc                                                                                                                                          >  8:112775/1‑149 (MQ=255)
         tcTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAg                                                                                                                                      >  6:17172/1‑149 (MQ=255)
             gccgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAgg                                                                                                                                     <  5:151451/146‑1 (MQ=255)
              ccgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAATCAGTTTGCCgctg                                                                                                                                      <  5:292914/144‑4 (MQ=255)
              ccgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAATCAGTTTGCCgctg                                                                                                                                      >  6:292914/1‑141 (MQ=255)
               cgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGGGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGgtggtg                                                                                                                                >  3:234078/1‑149 (MQ=255)
                gccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGgtggtg                                                                                                                                <  2:231225/148‑1 (MQ=255)
                gccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGgtggtg                                                                                                                                >  1:231225/1‑148 (MQ=255)
                      ggTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGATGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAgg                                                                                                                         <  1:92947/149‑1 (MQ=255)
                      ggTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAgg                                                                                                                         <  7:241379/149‑1 (MQ=255)
                      ggTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAgg                                                                                                                         <  3:84080/149‑1 (MQ=255)
                      ggTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAg                                                                                                                          <  5:17172/148‑1 (MQ=255)
                             aCGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGGATACGGCTGACCAGCGGCAAACTGATTGCCCTAGGTGGTGGCCGAGGTCAACAg                                                                                                                  >  1:54439/1‑149 (MQ=255)
                               gTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTTGTGAAGGTGACAGAATACGGCTGACCAGCGTCAAACTGATTTCCGTAGTTGGTGGCCGAGGTCAACAGgc                                                                                                                <  6:316852/149‑1 (MQ=255)
                               gTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGgc                                                                                                                <  7:112775/149‑1 (MQ=255)
                                tCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGgcc                                                                                                               >  7:255780/1‑149 (MQ=255)
                                  gAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGgccgc                                                                                                             >  6:313692/1‑149 (MQ=255)
                                    aCTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGAt                                                                                                                                              >  2:40132/1‑114 (MQ=255)
                                    aCTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGAt                                                                                                                                              <  1:40132/114‑1 (MQ=255)
                                    aCTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACa                                                                                                                   >  3:71281/1‑141 (MQ=255)
                                    aCTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACa                                                                                                                   <  4:71281/141‑1 (MQ=255)
                                    aCTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCa                                                                                                           <  5:143575/149‑1 (MQ=255)
                                    aCTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCa                                                                                                           <  7:26074/149‑1 (MQ=255)
                                      ttgcccttgccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCg                                                                                                                            >  6:220132/1‑130 (MQ=255)
                                      ttgcccttgccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCg                                                                                                                            <  5:220132/130‑1 (MQ=255)
                                         cccttgccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAg                                                                                                      <  4:192593/149‑1 (MQ=255)
                                            ttgccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAATCAGTTTGCCgct                                                                                                                                       <  4:324205/113‑3 (MQ=255)
                                            ttgccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAATCAGTTTGCCgct                                                                                                                                       >  3:324205/1‑111 (MQ=255)
                                            ttgccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCaa                                                                                                   >  2:223235/1‑149 (MQ=255)
                                            ttgccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCaa                                                                                                   >  2:191317/1‑149 (MQ=255)
                                              gccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAAt                                                                                                  <  4:64572/148‑1 (MQ=255)
                                               ccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATAc                                                                                                >  8:261839/1‑149 (MQ=255)
                                                 ggcggTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCt                                                                                              >  7:12788/1‑149 (MQ=255)
                                                      tttCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCg                                                                                         <  1:191317/149‑1 (MQ=255)
                                                           ttGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGa                                                                                                         >  1:182396/1‑128 (MQ=255)
                                                           ttGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGa                                                                                                         <  2:182396/128‑1 (MQ=255)
                                                               tGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGc                                                                                <  3:184065/149‑1 (MQ=255)
                                                                gATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCg                                                                               >  5:155809/1‑149 (MQ=255)
                                                                   gTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACg                                                                            >  3:136752/1‑149 (MQ=255)
                                                                      aTTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGt                                                                          >  5:1093/1‑148 (MQ=255)
                                                                          ccGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTAc                                                                     <  6:104996/149‑1 (MQ=255)
                                                                           cggcggCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTAcc                                                                    <  6:155809/149‑1 (MQ=255)
                                                                              cggcATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAgg                                                                 <  4:234078/149‑1 (MQ=255)
                                                                                 cATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTc                                                              <  1:223235/149‑1 (MQ=255)
                                                                                   tCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGa                                                            >  6:1193/1‑149 (MQ=255)
                                                                                     ccTCACGCTGTACTCGGTGGCACTGGTGAAGGTGATAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACt                                                          >  6:193120/1‑149 (MQ=255)
                                                                                     ccTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGAc                                                           <  6:316086/148‑1 (MQ=255)
                                                                                        cACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTggg                                                       >  8:31297/1‑149 (MQ=255)
                                                                                         aCGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGt                                                      <  4:176792/149‑1 (MQ=255)
                                                                                                    ggTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCtgtt                                             <  8:12788/147‑1 (MQ=255)
                                                                                                         cACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTtggtg                                      <  8:255780/149‑1 (MQ=255)
                                                                                                          aCTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTtggtgg                                     >  5:258459/1‑149 (MQ=255)
                                                                                                              gTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTg                                 >  6:135385/1‑149 (MQ=255)
                                                                                                                 aaGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTc                              <  5:313692/149‑1 (MQ=255)
                                                                                                                  aGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTc                              >  5:257927/1‑148 (MQ=255)
                                                                                                                    gTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGa                            >  6:226439/1‑148 (MQ=255)
                                                                                                                      gACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGaaaa                         <  1:74026/149‑1 (MQ=255)
                                                                                                                       aCAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCtgttg                                            >  6:49717/1‑129 (MQ=255)
                                                                                                                        cAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATg                       <  7:31297/149‑1 (MQ=255)
                                                                                                                         aGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGc                      >  8:181562/1‑149 (MQ=255)
                                                                                                                                 gCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCg              <  5:1193/149‑1 (MQ=255)
                                                                                                                                   tGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGgt            >  1:209361/1‑149 (MQ=255)
                                                                                                                                        aGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGGTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGt       <  2:54439/149‑1 (MQ=255)
                                                                                                                                          cGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTg     <  5:135385/149‑1 (MQ=255)
                                                                                                                                            gCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGt   >  1:132462/1‑149 (MQ=255)
                                                                                                                                             cAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTc  >  6:181052/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTC  >  NC_002947/4972737‑4973026

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