Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,944,743 A→G 5.5% intergenic (+103/+695) PP_5673 → / ← flhB D‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,944,7430AG 5.5% 216.9 / 2.1 74intergenic (+103/+695)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (31/38);  new base G (2/2);  total (33/41)
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.75e-01

ACTGGTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACCAGCAGCACCAC  >  NC_002947/4944595‑4944889
                                                                                                                                                     |                                                                                                                                                  
aCTGGTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCaaa                                                                                                                                                    <  5:272660/149‑1 (MQ=255)
    gTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGgc                                                                                                                                                <  2:3195/149‑1 (MQ=255)
       cTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGgcgc                                                                                                                                              <  5:82447/148‑1 (MQ=255)
           tCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGTTTGAGCGGCGCttaa                                                                                                                                          <  1:102916/148‑4 (MQ=255)
               aaaCAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCtt                                                                                                                                      >  6:11001/1‑148 (MQ=255)
                aacaacTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTGGTACAATCAGGCTATTGGATGTGGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCGAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATGtt                                                                                                                                      >  2:176012/1‑147 (MQ=255)
                aacaacTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGa                                                                                                                                    <  7:32447/149‑1 (MQ=255)
                 acaacTTGGATCCTACCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATAATGCTTTTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCCGCCGCTCAAACGGCGCTCGATCTTGa                                                                                                                                    <  1:15371/149‑1 (MQ=255)
                     cTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGTCGCTCGATCTTGAGCAAt                                                                                                                               >  3:182396/1‑149 (MQ=255)
                        ggATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATc                                                                                                                              <  1:184567/147‑1 (MQ=255)
                          aTCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCa                                                                                                                           <  5:126157/148‑1 (MQ=255)
                                gCCAATACTGCTGACCGAACAACAAGTATAATTCGAACCATCAGGCTATAGCATGTCGATGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACgg                                                                                                                    >  3:84214/1‑149 (MQ=255)
                                  cAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACgg                                                                                                                                                 <  6:333495/118‑1 (MQ=255)
                                   aaTACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGa                                                                                                                                    <  1:201714/130‑1 (MQ=255)
                                   aaTACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGa                                                                                                                                    >  2:201714/1‑130 (MQ=255)
                                         gctgACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAAc                                                                                                           >  4:45598/1‑149 (MQ=255)
                                               cGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCaa                                                                                                     <  3:340052/149‑1 (MQ=255)
                                               cGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCaa                                                                                                     >  5:273980/1‑149 (MQ=255)
                                                gAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAAc                                                                                                    >  4:183022/1‑149 (MQ=255)
                                                   caacaaTTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTg                                                                                                  <  4:84214/148‑1 (MQ=255)
                                                        aGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGACGCTCAAACGGCGCTCGATCTTGAGCAATCCAAGAAACGGTCAAGAAACTTTCAACTGTCAcga                                                                                            >  5:316097/1‑148 (MQ=255)
                                                             aaTTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTc                                                                                                <  2:298740/140‑1 (MQ=255)
                                                             aaTTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTc                                                                                                >  1:298740/1‑140 (MQ=255)
                                                                      cATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACCGTCACGCCGCAAATCcggcg                                                                               >  3:196425/1‑148 (MQ=255)
                                                                       aTCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCcggcgg                                                                              >  8:54236/1‑148 (MQ=255)
                                                                        tCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGaa                                                                            <  6:273980/149‑1 (MQ=255)
                                                                          aGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCaa                                                                                                                >  4:97102/1‑111 (MQ=255)
                                                                          aGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCaa                                                                                                                <  3:97102/111‑1 (MQ=255)
                                                                          aGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACt                                                                                                   >  2:118380/1‑124 (MQ=255)
                                                                          aGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACt                                                                                                   <  1:118380/124‑1 (MQ=255)
                                                                          aGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATc                                                                          <  5:11001/149‑1 (MQ=255)
                                                                          aGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACCGTCACGCCGCAAATCCGGCGGAATc                                                                          <  4:196425/149‑1 (MQ=255)
                                                                           ggCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCg                                                                         >  6:334203/1‑149 (MQ=255)
                                                                           ggCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCg                                                                         >  2:170080/1‑149 (MQ=255)
                                                                                ttGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGaaa                                                                     <  7:96121/148‑1 (MQ=255)
                                                                                   gATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCa                                                                  <  2:46595/148‑1 (MQ=255)
                                                                                   gATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAg                                                                 <  3:183022/149‑1 (MQ=255)
                                                                                     tGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGa                                                               <  1:170080/149‑1 (MQ=255)
                                                                                                tttttCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATg                                                             <  3:221277/140‑1 (MQ=255)
                                                                                                tttttCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATg                                                             >  4:221277/1‑140 (MQ=255)
                                                                                                          ggTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGaaa                                                                     >  7:143868/1‑122 (MQ=255)
                                                                                                          ggTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGaaa                                                                     <  8:143868/122‑1 (MQ=255)
                                                                                                            tAATTTATTGATCACCCAATTGAGATTAAGAGCCGCTCAAACGGCGCTCGATCTTGAGCAATACCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCGGGATGGTCGGGAACGCGGCAAAGTc                                         >  8:200221/1‑148 (MQ=255)
                                                                                                              aTTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGa                                                                       <  5:116209/116‑1 (MQ=255)
                                                                                                              aTTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGa                                                                       >  6:116209/1‑116 (MQ=255)
                                                                                                              aTTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCgg                                                <  2:262454/139‑1 (MQ=255)
                                                                                                              aTTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCgg                                                >  1:262454/1‑139 (MQ=255)
                                                                                                               tttATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCTGGATGGTCGGGAACGCGGCAAAGTCGAGc                                     <  7:200221/149‑1 (MQ=255)
                                                                                                                  aTTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAgg                                                                >  7:302661/1‑119 (MQ=255)
                                                                                                                  aTTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAgg                                                                <  8:302661/119‑1 (MQ=255)
                                                                                                                   ttGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCggg                                  <  5:334203/148‑1 (MQ=255)
                                                                                                                   ttGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGgc                                 >  5:185041/1‑149 (MQ=255)
                                                                                                                   ttGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGgc                                 <  4:182396/149‑1 (MQ=255)
                                                                                                                     gATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAAcgc                                                  >  6:364485/1‑130 (MQ=255)
                                                                                                                     gATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAAcgc                                                  <  5:364485/130‑1 (MQ=255)
                                                                                                                      aTCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGt                                          >  7:208230/1‑137 (MQ=255)
                                                                                                                      aTCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGt                                          <  8:208230/137‑1 (MQ=255)
                                                                                                                         aCCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTc                                                          >  6:86544/1‑118 (MQ=255)
                                                                                                                         aCCCAATTGAGATTAAGCGCCGCCCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGATCGGGCGCAAg                            <  4:162102/148‑1 (MQ=255)
                                                                                                                             aaTTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAANTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGta                       <  1:111860/149‑1 (MQ=255)
                                                                                                                             aaTTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGta                       >  3:232777/1‑149 (MQ=255)
                                                                                                                             aaTTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGta                       <  1:115533/149‑1 (MQ=255)
                                                                                                                                 gagaTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATaca                   >  3:182194/1‑149 (MQ=255)
                                                                                                                                  agaTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATaca                   <  7:54236/148‑1 (MQ=255)
                                                                                                                                    atcgaGCGCCGTTTGAGCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATac                    >  6:354237/5‑145 (MQ=38)
                                                                                                                                    atcgaGCGCCGTTTGAGCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATac                    <  5:354237/141‑1 (MQ=38)
                                                                                                                                    aTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATaca                   >  4:277076/1‑146 (MQ=255)
                                                                                                                                       gaGCGCCGTTTGAGCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAAc                                                    >  4:153973/2‑110 (MQ=255)
                                                                                                                                       aaGCGCCTCTCAATCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGTCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGAc              <  6:316097/148‑1 (MQ=255)
                                                                                                                                       aaGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGAc              <  3:45598/148‑1 (MQ=255)
                                                                                                                                       aaGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGAc              <  3:277076/148‑1 (MQ=255)
                                                                                                                                          ccgccgcTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGGCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACcag           <  1:176012/148‑1 (MQ=255)
                                                                                                                                            gccgcTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACcag           >  1:150693/1‑147 (MQ=255)
                                                                                                                                                cTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACca            >  2:44874/1‑142 (MQ=255)
                                                                                                                                                cTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACca            <  1:44874/142‑1 (MQ=255)
                                                                                                                                                   aaaCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGc                          >  2:80840/1‑125 (MQ=255)
                                                                                                                                                   aaaCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACCAGCAGcaccac  >  5:229221/1‑149 (MQ=255)
                                                                                                                                                   aaaCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACCAGCAGcaccac  >  6:344718/1‑149 (MQ=255)
                                                                                                                                                   aaaCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTACAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGc                          <  1:80840/125‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACTGGTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTATTGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACCAGCAGCACCAC  >  NC_002947/4944595‑4944889

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: