Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,967,786 C→T 5.2% G107S (GGT→AGT)  PP_4377 ← flagellin FlaG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,967,7860CT 5.2% 237.9 / 2.1 79G107S (GGT→AGT) PP_4377flagellin FlaG
Reads supporting (aligned to +/- strand):  ref base C (42/30);  new base T (2/2);  total (46/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

GGAGCCGATACCACTGACTGTTGTACTTGCCATTGCTTTATATCCCTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAACTGCCGCTGCGAAGCT  >  NC_002947/4967644‑4967934
                                                                                                                                              |                                                                                                                                                    
ggAGCCGATACCACTGACTGTTGTACTTGCCATTGCTTTATATCCCTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCAcc                                                                                                                                                <  2:263600/149‑1 (MQ=255)
          ccACTGACTGTTGTACTTGCCATTGCTTTATATCCCTGATCAGCGAATGCTAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGa                                                                                                                                      <  2:325980/149‑1 (MQ=255)
                actgTTGTACTTGCCATTGCTTTATATCCCTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGc                                                                                                                                  <  2:204200/147‑1 (MQ=255)
                  tgttgtACTTGCCATTGCTTTATATCCCTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGc                                                                                                                              <  1:340238/149‑1 (MQ=255)
                  tgttgtACTTGCCATTGCTTTATATCCCTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGc                                                                                                                              <  1:143596/149‑1 (MQ=255)
                    ttgtACTTGCCATTGCTTTATATCCCTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCa                                                                                                                                                     <  7:346818/124‑1 (MQ=255)
                    ttgtACTTGCCATTGCTTTATATCCCTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCa                                                                                                                                                     >  8:346818/1‑124 (MQ=255)
                            gCCATTGCTTTATATCCCTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAgg                                                                                                                    >  6:242270/1‑149 (MQ=255)
                               aTTGCTTTATATCCCTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAg                                                                                                                     >  5:397334/1‑145 (MQ=255)
                               aTTGCTTTATATCCCTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAg                                                                                                                     <  6:397334/145‑1 (MQ=255)
                                           cccTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACAATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTgg                                                                                                      >  1:305650/1‑148 (MQ=255)
                                            ccTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTgg                                                                                                      >  6:146716/1‑147 (MQ=255)
                                                 tCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGa                                                                                                                                <  7:249242/116‑1 (MQ=255)
                                                 tCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGa                                                                                                                                >  8:249242/1‑116 (MQ=255)
                                                 tCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGt                                                                                               >  5:227253/1‑149 (MQ=255)
                                                   aGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGATGTc                                                                                              >  8:368668/1‑148 (MQ=255)
                                                     cGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTc                                                                                              <  3:140774/146‑1 (MQ=255)
                                                     cGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTc                                                                                              >  4:140774/1‑146 (MQ=255)
                                                     cGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGaa                                                                                           >  8:46124/1‑149 (MQ=255)
                                                     cGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGaa                                                                                           <  8:116347/149‑1 (MQ=255)
                                                     cGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGaa                                                                                           >  8:39948/1‑149 (MQ=255)
                                                      gAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAAt                                                                                          <  7:39948/149‑1 (MQ=255)
                                                      gAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAAt                                                                                          <  7:46124/149‑1 (MQ=255)
                                                      gAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAAt                                                                                          <  5:242270/149‑1 (MQ=255)
                                                        aTGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGa                                                                                        <  4:149928/149‑1 (MQ=255)
                                                         tGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTGCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGAc                                                                                       <  5:1113/149‑1 (MQ=255)
                                                            cAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCATCACTTGGAAGCTGTCGAATGACTTc                                                                                    >  8:29610/1‑149 (MQ=255)
                                                                   ttGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACt                                                                              >  2:403413/1‑148 (MQ=255)
                                                                   ttGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACAATCACCCTCCGCCAGACTGCGACCCAGGCTCAGCGCTGCATCACATGGAAGCTGTCGAATGACTTCGCCACTa                                                                             >  8:17532/1‑149 (MQ=255)
                                                                    tGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTAc                                                                            >  5:329745/1‑149 (MQ=255)
                                                                        gATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCCTAAATTAAAAAACAATCACCCTCCGACATACTGCGAGCCAGGCTCAGGGCTGACTCACTTGTAAGCTGTCTAATGACTTCGCCACTACCgtt                                                                        <  5:211962/149‑3 (MQ=255)
                                                                        gATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGt                                                                        >  6:366722/1‑149 (MQ=255)
                                                                        gATGAATGAGCGAGACACAATCCGAGGCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGt                                                                        <  5:211961/149‑1 (MQ=255)
                                                                          tgaatgaGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTgg                                                                      <  2:305650/149‑1 (MQ=255)
                                                                            aatgaGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAAt                                                                                          >  6:176848/1‑127 (MQ=255)
                                                                              tgaGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAAt                                                                  <  6:227253/149‑1 (MQ=255)
                                                                               gaGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAAt                                                                  >  4:98870/1‑148 (MQ=255)
                                                                                aGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATa                                                                 <  7:368668/148‑1 (MQ=255)
                                                                                    agaCACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCtt                                                            <  5:366722/149‑1 (MQ=255)
                                                                                       cacaATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTTac                                                         >  2:343562/1‑149 (MQ=255)
                                                                                        acaATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGac                                                         >  5:389287/1‑148 (MQ=255)
                                                                                          aaTCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGaccac                                                      >  8:38013/1‑149 (MQ=255)
                                                                                            tCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGaccac                                                      >  4:64076/1‑147 (MQ=255)
                                                                                              cGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAgg                                                  <  6:389287/149‑1 (MQ=255)
                                                                                                  tCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTcc                                              >  6:73156/1‑149 (MQ=255)
                                                                                                    aaGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTcccc                                            >  4:25605/1‑149 (MQ=255)
                                                                                                    aaGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTcccc                                            >  4:220288/1‑149 (MQ=255)
                                                                                                    aaGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTcccc                                            >  1:395162/1‑149 (MQ=255)
                                                                                                       aCTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTg                                         >  3:159337/1‑149 (MQ=255)
                                                                                                         tCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTgg                                        >  5:141717/1‑148 (MQ=255)
                                                                                                           aTGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAAt                                     >  5:278923/1‑149 (MQ=255)
                                                                                                           aTGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAAt                                     >  2:411072/1‑149 (MQ=255)
                                                                                                            tGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATc                                    >  8:156425/1‑149 (MQ=255)
                                                                                                              ccAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGt                                  <  7:38013/149‑1 (MQ=255)
                                                                                                              ccAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCATCAATTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGTAATCGt                                  >  6:331269/1‑149 (MQ=255)
                                                                                                               cAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTc                                 <  5:73156/149‑1 (MQ=255)
                                                                                                                aaCATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACAACCAGGGTCCCCGTGGAATCGTc                                 >  3:358103/1‑148 (MQ=255)
                                                                                                                   aTTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACg                             >  8:165541/1‑149 (MQ=255)
                                                                                                                    ttACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACgg                            <  3:98870/149‑1 (MQ=255)
                                                                                                                          aTTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAAt                      >  5:129828/1‑149 (MQ=255)
                                                                                                                           ttGTCATCAAGTAAAAATGCATCACACTCCGCCAGAGTGCGAGCCAGTCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCAGTACCGGTGGGAATAACCTTGTCCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATc                     <  7:17532/149‑1 (MQ=255)
                                                                                                                               catcaAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAAc                 >  6:310815/1‑149 (MQ=255)
                                                                                                                                 tcaAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTc                                                                                              >  1:169565/1‑70 (MQ=255)
                                                                                                                                 tcaAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTc                                                                                              <  2:169565/70‑1 (MQ=255)
                                                                                                                                 tcaAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGAccacca                                                    >  3:191433/1‑112 (MQ=255)
                                                                                                                                 tcaAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGAccacca                                                    <  4:191433/112‑1 (MQ=255)
                                                                                                                                    aGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATc                                    <  4:189923/125‑1 (MQ=255)
                                                                                                                                    aGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATc                                    >  3:189923/1‑125 (MQ=255)
                                                                                                                                      tAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAACTGCCGCt          >  8:161309/1‑149 (MQ=255)
                                                                                                                                       aaaaaaCCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTcc                                <  2:391178/126‑1 (MQ=255)
                                                                                                                                       aaaaaaCCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTcc                                >  1:391178/1‑126 (MQ=255)
                                                                                                                                         aaaaCCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAACTGCCGCTGc        >  2:127393/1‑148 (MQ=255)
                                                                                                                                          aaaCCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAACTGCCGCTGCg       >  8:198295/1‑148 (MQ=255)
                                                                                                                                           aaCTATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAACTGc              <  4:103358/140‑1 (MQ=255)
                                                                                                                                           aaCTATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAACTGc              >  3:103358/1‑140 (MQ=255)
                                                                                                                                           aaCTATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAACTGc              >  3:97615/1‑140 (MQ=255)
                                                                                                                                           aaCTATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAACTGc              <  4:97615/140‑1 (MQ=255)
                                                                                                                                           aaCCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCACCGTGGAATCGTCCACGGAAAAATCCAACTGCCGCTgag       >  5:326034/1‑145 (MQ=255)
                                                                                                                                           aaCCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTCACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAACTGCCGCTGCGaa     <  1:343562/149‑1 (MQ=255)
                                                                                                                                              cATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAACTGCCGCTGCGAAGCt  >  4:97892/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGAGCCGATACCACTGACTGTTGTACTTGCCATTGCTTTATATCCCTGATCAGCGAATGCCAATTGCTTGACGATGAATGAGCGAGACACAATCCGAGTCAAGACTCATGCCAACATTACACATTGTCATCAAGTAAAAAACCATCACCCTCCGCCAGACTGCGAGCCAGGCTCAGGGCTGCCTCACTTGGAAGCTGTCGAATGACTTCGCCACTACCGGTGGCAATAACCTTGACCACCAGGGTCCCCGTGGAATCGTCCACGGAAAAATCCAACTGCCGCTGCGAAGCT  >  NC_002947/4967644‑4967934

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: ‑