Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,770,438 A→T 6.2% L145H (CTT→CAT)  PP_5060 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,770,4380AT 6.2% 188.6 / 2.6 64L145H (CTT→CAT) PP_5060hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/26);  new base T (2/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

ACGCAGCAGCTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTC  >  NC_002947/5770290‑5770586
                                                                                                                                                    |                                                                                                                                                    
aCGCAGCAGCTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCaaa                                                                                                                                                      <  7:109300/149‑1 (MQ=255)
   cagcagCTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAagca                                                                                                                                                   <  6:134907/149‑1 (MQ=255)
   cagcagCTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAagca                                                                                                                                                   <  1:297988/149‑1 (MQ=255)
     gcagcTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCTGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAagcagc                                                                                                                                                 <  1:185031/149‑1 (MQ=255)
        gcTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCAtt                                                                                                                                              <  2:313349/149‑1 (MQ=255)
         cTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGCGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGATCCAATAGTCGATCTGCGCTTGCGGGAAGTGCCAAACCAGCATTg                                                                                                                                             >  6:262667/1‑149 (MQ=255)
            gCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCAGAGGTCCCGGTCGATCGTGCCGTCGACGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGgc                                                                                                                                           <  5:388273/148‑1 (MQ=255)
              gttgttCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGgcag                                                                                                                                         >  4:157217/1‑148 (MQ=255)
                        aTCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTc                                                                                                                               <  5:307381/148‑1 (MQ=255)
                         tCCAGCGCTTCATCCTTCAGGCCATCGAGAAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGGCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAg                                                                                                                             <  4:41191/149‑1 (MQ=255)
                              cgcTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCAt                                                                                                                        <  3:157217/149‑1 (MQ=255)
                                 ttCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAATGCTGCTTTGgcacttcca                                                                                                                                    <  2:197658/134‑7 (MQ=255)
                                 ttCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATATGCGCTTGCTGGAAGTGCCAATGCTGCTTTGgcacttcca                                                                                                                                    >  1:197658/1‑128 (MQ=255)
                                    aTCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCa                                                                                                                  <  3:311853/149‑1 (MQ=255)
                                       cTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCAt                                                                                                                                               <  6:217303/117‑1 (MQ=255)
                                       cTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCAt                                                                                                                                               >  5:217303/1‑117 (MQ=255)
                                        ttCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGgcg                                                                                                              >  8:292694/1‑149 (MQ=255)
                                         tCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGgcgc                                                                                                             >  2:197575/1‑149 (MQ=255)
                                                 tCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTc                                                                                                     >  7:51009/1‑149 (MQ=255)
                                                  cGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTc                                                                                                     >  2:289684/1‑148 (MQ=255)
                                                         gTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCg                                                                                                                      <  2:371750/124‑1 (MQ=255)
                                                         gTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATAGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCg                                                                                                                      >  1:371750/1‑124 (MQ=255)
                                                               gcgcAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCt                                                                                       >  7:94015/1‑149 (MQ=255)
                                                               gcgcAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTCCCAAACCACCATTGGCAGCAACGCGTACAGGCCAACGTACAGGCCCCCGGCCACGAACAAATCGcccg                                                                                       >  2:167857/1‑148 (MQ=255)
                                                                   aacaacTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGc                                                                                   <  5:185214/149‑1 (MQ=255)
                                                                    acaacTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGcc                                                                                  <  7:276834/149‑1 (MQ=255)
                                                                    acaacTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGcc                                                                                  >  3:85135/1‑149 (MQ=255)
                                                                              cAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGc                                                                        >  5:278651/1‑149 (MQ=255)
                                                                              cAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGc                                                                        <  2:369055/149‑1 (MQ=255)
                                                                                aCCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGctc                                                                      >  7:50690/1‑149 (MQ=255)
                                                                                     aGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGctcgactc                                                                 >  7:227920/1‑149 (MQ=255)
                                                                                       gTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGctcgact                                                                  <  2:48810/146‑1 (MQ=255)
                                                                                       gTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGctcgact                                                                  >  1:48810/1‑146 (MQ=255)
                                                                                        tCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGctcgactcga                                                               <  6:278651/148‑1 (MQ=255)
                                                                                        tCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGAt                                                              >  3:329122/1‑149 (MQ=255)
                                                                                         cccGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTAAGGCCATCGTTCAGGCGCGCGGGCTCGATCAAATCGCCGTCAGCCGTTTCGATGCTCGAATAGATg                                                             >  7:325032/1‑149 (MQ=255)
                                                                                          ccGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATgg                                                            >  4:28034/1‑149 (MQ=255)
                                                                                           cGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGc                                                           >  5:379337/1‑149 (MQ=255)
                                                                                            ggTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGc                                                           >  1:227346/1‑148 (MQ=255)
                                                                                             gTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCg                                                          >  3:132049/1‑148 (MQ=255)
                                                                                                  tCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGACAAAGCAGCATTGGCAGCAAAGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAg                                                    >  1:343216/1‑149 (MQ=255)
                                                                                                    gTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAgc                                                   <  8:227920/148‑1 (MQ=255)
                                                                                                               gTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGt                                       <  3:28034/149‑1 (MQ=255)
                                                                                                                ttCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGtt                                      <  4:85135/149‑1 (MQ=255)
                                                                                                                   cAATAGTCGATCTGCGCTTGCTGGAAGTGCGAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAg                                   >  1:30776/1‑149 (MQ=255)
                                                                                                                   cAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAg                                   <  2:30776/149‑1 (MQ=255)
                                                                                                                   cAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAg                                   <  1:197575/149‑1 (MQ=255)
                                                                                                                    aaTAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAgg                                  >  7:114498/1‑149 (MQ=255)
                                                                                                                                gcgcTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAg                      >  4:219140/1‑149 (MQ=255)
                                                                                                                                 cgcTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGttt                                                                              >  8:20560/1‑92 (MQ=255)
                                                                                                                                 cgcTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGttt                                                                              <  7:20560/92‑1 (MQ=255)
                                                                                                                                 cgcTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCg                          <  2:227346/144‑1 (MQ=255)
                                                                                                                                  ggttGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTc                    <  1:167857/147‑1 (MQ=255)
                                                                                                                                  gcTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTa                               <  8:246591/138‑1 (MQ=255)
                                                                                                                                  gcTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTa                               >  7:246591/1‑138 (MQ=255)
                                                                                                                                   cTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAATCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTcg                   <  1:289684/149‑1 (MQ=255)
                                                                                                                                      gCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTcgccc                >  4:135952/1‑149 (MQ=255)
                                                                                                                                         gttgctGCCAATGCTGCTTTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCa                                    >  4:345626/6‑126 (MQ=255)
                                                                                                                                         gttgctGCCAATGCTGCTTTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCa                                    <  3:345626/121‑1 (MQ=255)
                                                                                                                                          gAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCgg                                         <  7:182885/120‑1 (MQ=255)
                                                                                                                                          gAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCgg                                         >  8:182885/1‑120 (MQ=255)
                                                                                                                                            aGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCggg          >  6:12635/1‑149 (MQ=255)
                                                                                                                                            aGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCggg          >  2:304893/1‑149 (MQ=255)
                                                                                                                                              tGCCAAAGCAGTATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGc        >  8:65163/1‑149 (MQ=255)
                                                                                                                                              tGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGCCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGc        >  1:329345/1‑149 (MQ=255)
                                                                                                                                                ccAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTa                       >  7:102447/1‑132 (MQ=255)
                                                                                                                                                ccAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTa                       <  8:102447/132‑1 (MQ=255)
                                                                                                                                                ccAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGa      <  6:379337/149‑1 (MQ=255)
                                                                                                                                                  aaaGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCGGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAgg    <  8:114498/149‑1 (MQ=255)
                                                                                                                                                   aaGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGt   >  2:286040/1‑149 (MQ=255)
                                                                                                                                                    agcagcATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCTCCGGGGCGAGGTc  <  8:51009/149‑1 (MQ=255)
                                                                                                                                                    agcagcATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTc  >  1:296298/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGCAGCAGCTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGCAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGGTCAGGCGCGCGGGCTCGATCAAATCGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGGTC  >  NC_002947/5770290‑5770586

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