Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,650,987 G→C 6.2% A62P (GCC→CCC) ‡ PP_2323 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,650,9870GC 6.2% 211.0 / 4.4 65A62P (GCC→CCC) ‡PP_2323hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/27);  new base C (2/2);  total (36/29)
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.49e-01

ATTACAATCGTCCAATGGTCGAAGCGGCAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGTGGGATCGCTGCCGCCCTGACCGAGCGGG  >  NC_002947/2650841‑2651134
                                                                                                                                                  |                                                                                                                                                   
aTTACAATCGTCCAATGGTCGAAGCGGCAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGcc                                                                                                                                                   <  5:195335/149‑1 (MQ=255)
 ttACAATCGTCCAATGGTCGAAGCGGCAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCg                                                                                                                                                  <  4:85607/149‑1 (MQ=255)
   aCAATCGTCCAATGGTCGAAGCGGCAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGt                                                                                                                                                >  3:19331/1‑149 (MQ=255)
            cAATGGTCGAAGCGGCAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTc                                                                                                                                       >  4:37889/1‑149 (MQ=255)
                 gTCGAAGCGGCAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGAc                                                                                                                                  <  5:357318/149‑1 (MQ=255)
                     aaGCGGCAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAAc                                                                                                                              <  8:388559/149‑1 (MQ=255)
                          gCAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCgaga                                                                                                                                   >  5:169436/1‑139 (MQ=255)
                          gCAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCgaga                                                                                                                                   <  6:169436/139‑1 (MQ=255)
                          gCAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGa                                                                                                                         >  5:330156/1‑149 (MQ=255)
                           cAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGAc                                                                                                                        <  6:330156/149‑1 (MQ=255)
                                    gCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGt                                                                                                               <  3:37889/149‑1 (MQ=255)
                                       tAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGc                                                                                                             >  4:3675/1‑148 (MQ=255)
                                       tAGAGGATTCTGTGCAGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGCCCCTCCACCATGACCAGGCGGGTCATCAGTTCGAGCGCAGCGTGGAGGGGGATCGTGcg                                                                                                            >  8:68942/1‑148 (MQ=255)
                                            gATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACtt                                                                                                       >  4:116274/1‑149 (MQ=255)
                                                 tgtgCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACg                                                                                                   >  4:135886/1‑148 (MQ=255)
                                                   tgCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACg                                                                                                   <  3:135886/146‑1 (MQ=255)
                                                       ggATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGtcatca                                                                                                                                           <  2:235880/102‑1 (MQ=255)
                                                       ggATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGtcatca                                                                                                                                           >  1:235880/1‑102 (MQ=255)
                                                        gATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGc                                                                                                             >  6:45335/1‑131 (MQ=255)
                                                        gATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGc                                                                                                             <  5:45335/131‑1 (MQ=255)
                                                             accaccATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTgg                                                                                      >  8:11155/1‑149 (MQ=255)
                                                              ccaccaTCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTggg                                                                                     <  4:53810/149‑1 (MQ=255)
                                                              ccaccaTCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTggg                                                                                     <  8:42939/149‑1 (MQ=255)
                                                               caccaTCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGt                                                                                                  <  8:99149/135‑1 (MQ=255)
                                                               caccaTCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGt                                                                                                  >  7:99149/1‑135 (MQ=255)
                                                               caccaTCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTggg                                                                                     >  2:214309/1‑148 (MQ=255)
                                                                 ccaTCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTggac                                                                                    >  3:76566/1‑145 (MQ=255)
                                                                    tCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCa                                                                               <  2:226563/149‑1 (MQ=255)
                                                                      ggTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGa                                                                             <  7:11155/149‑1 (MQ=255)
                                                                       gTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGAc                                                                            >  2:14286/1‑149 (MQ=255)
                                                                       gTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGACAAGCAGAc                                                                            <  4:76566/149‑1 (MQ=255)
                                                                        ttAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACg                                                                           <  4:19331/149‑1 (MQ=255)
                                                                          aaggtaaggGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGTCAGCAGACGCt                                                                         >  1:141322/1‑149 (MQ=255)
                                                                           aggtaaggGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTa                                                                                                          <  5:18918/115‑1 (MQ=255)
                                                                           aggtaaggGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTa                                                                                                          >  6:18918/1‑115 (MQ=255)
                                                                            ggtaaggGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGc                                                                                <  4:108491/140‑1 (MQ=255)
                                                                            ggtaaggGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGc                                                                                >  3:108491/1‑140 (MQ=255)
                                                                             gtaaggGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTgg                                                                       >  8:257599/1‑148 (MQ=255)
                                                                             gtaaggGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGa                                                                      <  3:3675/149‑1 (MQ=255)
                                                                                   gTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTa                                                                >  3:135076/1‑149 (MQ=255)
                                                                                    tCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTAt                                                               >  8:249727/1‑149 (MQ=255)
                                                                                        ccGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGt                                                                                                  >  8:222391/1‑110 (MQ=255)
                                                                                        ccGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGt                                                                                                  <  7:222391/110‑1 (MQ=255)
                                                                                         cGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACAATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATATGGGAAAGCAGACGCTGGACATCTATCGCAc                                                          >  3:401303/1‑149 (MQ=255)
                                                                                                     cGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCaa                                              <  4:135076/149‑1 (MQ=255)
                                                                                                           gtagtGATGAGCGAAGCGCTGACCATCCACCATGACGAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGt                                                    <  6:331449/137‑1 (MQ=255)
                                                                                                           gtagtGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCAc                                                          <  4:146275/131‑1 (MQ=255)
                                                                                                           gtagtGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCAc                                                          >  3:146275/1‑131 (MQ=255)
                                                                                                           gtagtGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGt                                                    >  5:331449/1‑137 (MQ=255)
                                                                                                            tagtGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTg                                       >  2:139324/1‑149 (MQ=255)
                                                                                                                 atgaGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACATTCGTGCCCAACGCCCTGCGCgg                                  >  3:304141/1‑149 (MQ=255)
                                                                                                                  tgaGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCgg                                  <  7:249727/148‑1 (MQ=255)
                                                                                                                    aGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCg                               >  8:317299/1‑149 (MQ=255)
                                                                                                                              tGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGTGGGATCGCt                     >  8:294773/1‑149 (MQ=255)
                                                                                                                                  cATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGTGGGATCGCTgccg                 >  3:311129/1‑149 (MQ=255)
                                                                                                                                   aTCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGTGGGATCGCTgccgc                >  7:245983/1‑149 (MQ=255)
                                                                                                                                   aTCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGTGGGATCGCTgccgc                >  1:410965/1‑149 (MQ=255)
                                                                                                                                        tgaTGACCGGCCTGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGt                              >  3:100972/3‑130 (MQ=255)
                                                                                                                                        tgaTGACCGGCCTGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGt                              <  4:100972/128‑1 (MQ=255)
                                                                                                                                        ccaTGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGTGGGATCGCTGCCGCCCTGa           >  5:154426/1‑149 (MQ=255)
                                                                                                                                          aTGACCGGCCTGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTg                                                   >  6:29920/1‑107 (MQ=255)
                                                                                                                                          aTGACCGGCCTGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTg                                                   <  5:29920/107‑1 (MQ=255)
                                                                                                                                             aCCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGTGGGATCGc                      <  1:19375/133‑1 (MQ=255)
                                                                                                                                             aCCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGTGGGATCGc                      >  2:19375/1‑133 (MQ=255)
                                                                                                                                               cAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGTGGGATCGCTGCCGCCCTGACCGAGCg    <  2:410965/149‑1 (MQ=255)
                                                                                                                                                 ggCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGTGGGATCGCTGCCGCCCTGACCGAGCggg  <  7:257599/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATTACAATCGTCCAATGGTCGAAGCGGCAACAGGTGGCATAGAGGATTCTGTGCCGGATAAACCACCATCGGTTAAGGTAAGGGTCTGCCGAGAAGACCTGCGAGGAGTAGTGATGAGCGAAGCGCTGACCATCCACCATGACCAGGCCGGTCATCAGTTCGAGACCAACGTGGACGGTCATCGTGCCTACTTGACGTACATGGATCTGGGCAAGCAGACGCTGGACATCTATCGCACCTTCGTGCCCAACGCCCTGCGCGGCCGTGGGATCGCTGCCGCCCTGACCGAGCGGG  >  NC_002947/2650841‑2651134

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