Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,944,738 C→T 7.8% intergenic (+98/+700) 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,7380CT 7.8% 179.4 / 4.3 65intergenic (+98/+700)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base T (3/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

ACAAACTGGTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACCAG  >  NC_002947/4944591‑4944880
                                                                                                                                                      |                                                                                                                                              
aCAAACTGGTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGcgccgc                                                                                                                                                <  4:81403/148‑1 (MQ=255)
   aaCTGGTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACACAATTGAGATTAAGCGCCGCTCa                                                                                                                                             >  6:111544/1‑148 (MQ=255)
      tGGTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACg                                                                                                                                         >  1:242877/1‑149 (MQ=255)
              aTCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGa                                                                                                                                 >  7:58689/1‑149 (MQ=255)
               tattAAACAACTTGGATACTAGAAAATACTGCTGACCGAACAACAAGTATAAGTAGTACCAGAAGGCTATTGGATGTCGTTGGTGTTTTTCGTAGGGTAATTTAT‑‑TGATAACCAAATTGAGATTAATCGACTCTCACACGGCTCTCGAt                                                                                                                                >  8:37250/3‑149 (MQ=255)
               tCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGa                                                                                                                                 >  3:244783/1‑148 (MQ=255)
                  tAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCtt                                                                                                                             >  1:88961/1‑149 (MQ=255)
                  tAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCtt                                                                                                                             >  2:111411/1‑149 (MQ=255)
                   aaaCAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGTTTGAGCGGCGCttaatc                                                                                                                               <  7:162115/146‑6 (MQ=255)
                   aaaCAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGAt                                                                                                                                <  4:244783/145‑1 (MQ=255)
                          ttGGATGCTATCCAATACTGCTGTCCTAACAACAATTATAATTCGTACCATCAGGCTATTTGATGTCGTTGGTGTTTTTCGGATTGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGcgccgc                                                                                                                                                <  8:135771/122‑1 (MQ=255)
                            ggATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATcc                                                                                                                    <  5:111544/148‑1 (MQ=255)
                              aTCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATcc                                                                                                                    >  6:219574/1‑146 (MQ=255)
                                  tAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAAc                                                                                                             >  7:17254/1‑149 (MQ=255)
                                    gCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAAc                                                                                                             <  2:242877/147‑1 (MQ=255)
                                      cAATACTGCTGACCGAACCTCAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTTGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCTAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTc                                                                                                         <  6:359703/149‑1 (MQ=255)
                                               tgACCGAACAACAAGTATAATTCGTACCATCAGGCTCTTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGa                                                                                                                           <  8:76456/122‑1 (MQ=255)
                                                 aCCTAAAAAAAATTATAATTTTTAACATAAGGCTTTTTTATGTCGTTTGTGTTTTTCGGATTTTAATTTATACTCTTCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGaaa                                                                                                   <  1:111411/146‑1 (MQ=255)
                                                  ccGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCa                                                                                             <  8:58689/149‑1 (MQ=255)
                                                   cGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCaa                                                                                            <  2:88961/149‑1 (MQ=255)
                                                   cGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATAACAAAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGAAATAACAGAAACGGTCAAGAAACATTCa                                                                                             >  5:27657/1‑148 (MQ=255)
                                                     aacaacAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACt                                                                                          >  2:207826/1‑149 (MQ=255)
                                                                 aaTTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCa                                                                               >  5:94191/1‑148 (MQ=255)
                                                                      gTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCcg                                                                         >  8:9893/1‑149 (MQ=255)
                                                                         acATCAGGATATTTGATTTCGTTGGTTTTTTTCGTATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGTCTCTCGATATTGATAAATCCCAGAAACGGTCAATAAACTTTCAACTTTCACGCAGCAAATCcggc                                                                       <  7:37250/147‑1 (MQ=255)
                                                                           aTCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCcggcgg                                                                     >  7:325774/1‑148 (MQ=255)
                                                                               ggCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCg                                                                >  2:58546/1‑149 (MQ=255)
                                                                                  tATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACGCCAATTGAGATTAAGCGCCGATCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCgg                                                               >  1:74010/1‑148 (MQ=255)
                                                                                  tATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGa                                                              >  1:323558/1‑148 (MQ=255)
                                                                                       gATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAg                                                        >  3:37679/1‑149 (MQ=255)
                                                                                        aTGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAgg                                                       <  2:323558/149‑1 (MQ=255)
                                                                                        aTGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAg                                                        >  7:207231/1‑148 (MQ=255)
                                                                                             gTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATgg                                                   <  6:94191/148‑1 (MQ=255)
                                                                                             gTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGt                                                  >  4:133808/1‑149 (MQ=255)
                                                                                              ttGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGa                                                              <  4:55783/136‑1 (MQ=255)
                                                                                              ttGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGa                                                              >  3:55783/1‑136 (MQ=255)
                                                                                              ttGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTc                                                 <  6:202642/149‑1 (MQ=255)
                                                                                                 gtgtTTTTCGGATTGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCgg                                               <  6:27657/148‑1 (MQ=255)
                                                                                                 gtgtTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCgg                                               <  2:74010/148‑1 (MQ=255)
                                                                                                    tttttCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAAc                                           >  1:55657/1‑149 (MQ=255)
                                                                                                      tttCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACTGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGACACGCAGCAAATCCGGAGGAATCCGATAGCAGGATGCTCGGGAAc                                           >  5:346823/1‑147 (MQ=255)
                                                                                                        tCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAAcgcg                                        <  7:9893/148‑1 (MQ=255)
                                                                                                             tGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAg                                  <  1:58546/149‑1 (MQ=255)
                                                                                                             tGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAg                                  <  5:128680/149‑1 (MQ=255)
                                                                                                                tAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTc                                <  3:288226/148‑1 (MQ=255)
                                                                                                                  aTTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAg                             <  5:219574/149‑1 (MQ=255)
                                                                                                                      aT‑‑TGATCACCCAATTGAGATTACGCGCCGCTCAACCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGGCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCgg                          <  1:336295/148‑1 (MQ=255)
                                                                                                                      aT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCggg                         >  8:198567/1‑149 (MQ=255)
                                                                                                                            aTCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAg                                                        >  7:75711/1‑114 (MQ=255)
                                                                                                                            aTCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAg                                                        <  8:75711/114‑1 (MQ=255)
                                                                                                                            aTCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGaa                                            <  8:60838/126‑1 (MQ=255)
                                                                                                                            aTCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGaa                                            >  7:60838/1‑126 (MQ=255)
                                                                                                                            aTCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGgcgc                      >  5:170222/1‑148 (MQ=255)
                                                                                                                               aCCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAg                   <  7:198567/148‑1 (MQ=255)
                                                                                                                                cccAATTGAGATTACGCGCCGCTCTATCGGCGGTCGATCTTGAGCAATCGGATAAACGGTCTAGAAACTTTCATCTGTCTCTCCGCAAATCCGGCGGATTCGGAAAGCAGGATTGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAgtc                 <  6:346823/149‑3 (MQ=255)
                                                                                                                                cccAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGc                 <  8:17254/149‑1 (MQ=255)
                                                                                                                                  ccAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCg                >  8:92268/1‑149 (MQ=255)
                                                                                                                                   cAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGt               <  4:37679/149‑1 (MQ=255)
                                                                                                                                      ttGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGtata            <  1:207826/149‑1 (MQ=255)
                                                                                                                                         agaTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATaca          >  8:258690/1‑148 (MQ=255)
                                                                                                                                          gatcgaGCGCCGTTTGAGCGGCGCTCGATCTTGAGTAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATaca          >  8:301336/6‑147 (MQ=38)
                                                                                                                                          gaTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATAcaca        <  7:92268/149‑1 (MQ=255)
                                                                                                                                           atcgaGCGCCGTTTGAGCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCcgg                                                                        <  3:355084/80‑1 (MQ=38)
                                                                                                                                           atcgaGCGCCGTTTGAGCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCcgg                                                                        >  4:355084/5‑84 (MQ=38)
                                                                                                                                           atcgaGCGCCGTTTGAGCGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATaca          >  6:67996/5‑146 (MQ=38)
                                                                                                                                           aTTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGACCGGGCGCAAGGCGTATAcaca        <  8:325774/148‑1 (MQ=255)
                                                                                                                                               aGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACca   <  6:170222/149‑1 (MQ=255)
                                                                                                                                                  gccgcTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACcag  >  1:358617/1‑147 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
ACAAACTGGTACTTATCTTAAACAACTTGGATCCTAGCCAATACTGCTGACCGAACAACAAGTATAATTCGTACCATCAGGCTATTGGATGTCGTTGGTGTTTTTCGGATGGTAATTTAT‑‑TGATCACCCAATTGAGATTAAGCGCCGCTCAAACGGCGCTCGATCTTGAGCAATCCCAGAAACGGTCAAGAAACTTTCAACTGTCACGCCGCAAATCCGGCGGAATCGGAAAGCAGGATGGTCGGGAACGCGGCAAAGTCGAGCGGGCGCAAGGCGTATACACAGACCAG  >  NC_002947/4944591‑4944880

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