Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 546,698 A→T 5.9% K17* (AAG→TAG) ‡ rpsG → 30S ribosomal protein S7

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947546,6980AT 5.9% 208.7 / 2.8 68K17* (AAG→TAG) ‡rpsG30S ribosomal protein S7
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base T (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

CAGACATCCATTTTTTTGAGTCGATAAGAGTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTC  >  NC_002947/546555‑546845
                                                                                                                                               |                                                                                                                                                   
cAGACATCCATTTTTTTGAGTCGATAAGAGTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTAc                                                                                                                                                >  4:58852/1‑149 (MQ=255)
 aGAAATCCATTTTTTTGAGTCTATAAGTGTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACg                                                                                                                                               <  8:187806/149‑1 (MQ=255)
     aTCCATTTTTTTGAGTCGATAAGAGTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATc                                                                                                                                           >  5:254883/1‑149 (MQ=255)
       ccATTTTTTTGAGTCGATAAGAGTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATccc                                                                                                                                         <  8:93714/149‑1 (MQ=255)
        cATTTTTTTGAGTCTATAAGAGTATGGGGGGGCTGGGGTCTAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGTACTTCGGATcgtc                                                                                                                                        <  4:187777/149‑4 (MQ=255)
        cATTTTTTTGAGTCGATAAGAGTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATcc                                                                                                                                          >  4:97873/1‑147 (MQ=255)
          tttttttGAGTCGATAAGAGTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGa                                                                                                                                      <  3:97873/149‑1 (MQ=255)
            tttttGATTACATAACAGAAAGGCCGGTCTGGGCCCGACCACCCACGTCACGTGATAAACTGAAGACCTTTTGACGGCTTATCATGCCAATACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGAt                                                                                                                                     <  6:343261/148‑1 (MQ=255)
                      gATAAGAGTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGTACTTCGGATcgt                                                                                                                                         <  4:286694/134‑3 (MQ=255)
                      gATAAGAGTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGTACTTCGGATcgt                                                                                                                                         >  3:286694/1‑132 (MQ=255)
                           gagTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCAt                                                                                                                      >  8:89596/1‑148 (MQ=255)
                            agTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATg                                                                                                                     >  8:10577/1‑148 (MQ=255)
                            agTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATg                                                                                                                     >  4:102252/1‑148 (MQ=255)
                             gTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGa                                                                                                                    <  3:58852/148‑1 (MQ=255)
                              tAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGCCGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGAATCCCAGATTCTCGCGAAGTTCATGAAc                                                                                                                  >  3:249145/1‑149 (MQ=255)
                                 ggTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCAc                                                                                                               >  6:142910/1‑149 (MQ=255)
                                      ggCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGATGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGAt                                                                                                          <  4:82904/149‑1 (MQ=255)
                                       gCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATcccg                                                                                                                                        >  1:139877/1‑117 (MQ=255)
                                       gCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATcccg                                                                                                                                        <  2:139877/118‑2 (MQ=255)
                                               cGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCgg                                                                                                 <  8:47472/149‑1 (MQ=255)
                                                      ccACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGAt                                                                                                                                            <  5:44622/99‑1 (MQ=255)
                                                      ccACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGAt                                                                                                                                            >  6:44622/1‑99 (MQ=255)
                                                      ccACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCaagaag                                                                                          <  3:91046/149‑1 (MQ=255)
                                                      ccACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCaagaag                                                                                          >  2:199833/1‑149 (MQ=255)
                                                      ccACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCaagaag                                                                                          >  1:351630/1‑149 (MQ=255)
                                                      ccACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGCCaagaag                                                                                          >  3:357105/1‑149 (MQ=255)
                                                      ccACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGCATCCCAGATTCTCGAGAAGTTCATGAACCACGTGATGGAAAGCGGCaagaag                                                                                          >  6:140103/1‑149 (MQ=255)
                                                                         gAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATtctc                                                                                                                                 >  2:358056/1‑91 (MQ=255)
                                                                         gAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATtctc                                                                                                                                 <  1:358056/91‑1 (MQ=255)
                                                                         gAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCg                                                                       >  8:99016/1‑149 (MQ=255)
                                                                          aaGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGt                                                                      >  7:105791/1‑149 (MQ=255)
                                                                             aCCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTa                                                                   <  4:357105/149‑1 (MQ=255)
                                                                                gTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCAAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACgg                                                                >  4:140212/1‑149 (MQ=255)
                                                                                 tttGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGt                                                               <  5:142910/149‑1 (MQ=255)
                                                                                      gggCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGGACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCt                                                          <  3:102252/149‑1 (MQ=255)
                                                                                       ggCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCaagaa                                                                                           <  1:151799/115‑1 (MQ=255)
                                                                                       ggCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCaagaa                                                                                           >  2:151799/1‑115 (MQ=255)
                                                                                       ggCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTg                                                         >  3:141058/1‑149 (MQ=255)
                                                                                          ttATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGa                                                       >  6:142846/1‑148 (MQ=255)
                                                                                            aTCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATAc                                                    >  8:86297/1‑149 (MQ=255)
                                                                                             tCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATAc                                                    <  6:254883/148‑1 (MQ=255)
                                                                                               aTGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGt                                                 >  2:113296/1‑149 (MQ=255)
                                                                                               aTGCCAAGACGTCGTGTAGCAGCAAAACGTGACAACCTTGACGATCCGAACTACGGATCCCAGAAAATCGCCAAGTTCATCAACCACGTGATGGAAAGCGCCAACAAGGCCGTAGCCCAGCGCATCGTTTACGGAGCCCTGGAAAccct                                                 >  4:315548/1‑147 (MQ=255)
                                                                                                   cAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCaa                                                                                              <  8:296159/100‑1 (MQ=255)
                                                                                                   cAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCaa                                                                                              >  7:296159/1‑100 (MQ=255)
                                                                                                   cAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCaaa                                             <  2:351630/149‑1 (MQ=255)
                                                                                                   cAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCa                                               >  1:328788/1‑147 (MQ=255)
                                                                                                         gtcgtGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGc                                       <  7:89596/149‑1 (MQ=255)
                                                                                                                 gcagcaAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCaaa                                             >  8:88965/1‑135 (MQ=255)
                                                                                                                 gcagcaAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCaaa                                             <  7:88965/135‑1 (MQ=255)
                                                                                                                 gcagcaAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAg                               <  7:99016/149‑1 (MQ=255)
                                                                                                                     caAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGAc                           >  4:352378/1‑149 (MQ=255)
                                                                                                                            tGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGc                                       >  5:182302/1‑130 (MQ=255)
                                                                                                                            tGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGc                                       <  6:182302/130‑1 (MQ=255)
                                                                                                                                aTCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGaaa                 >  4:131439/1‑148 (MQ=255)
                                                                                                                                aTCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAAt                >  7:86589/1‑149 (MQ=255)
                                                                                                                                aTCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAAt                >  3:173280/1‑149 (MQ=255)
                                                                                                                                  ccTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCaaa                                             <  6:166592/118‑1 (MQ=255)
                                                                                                                                  ccTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCaaa                                             >  5:166592/1‑118 (MQ=255)
                                                                                                                                   cTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCtt             <  2:328788/149‑1 (MQ=255)
                                                                                                                                     tGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAg                               >  6:45911/1‑129 (MQ=255)
                                                                                                                                     tGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAg                               <  5:45911/129‑1 (MQ=255)
                                                                                                                                      tgggATCCGTACTTCGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTc            >  3:144521/4‑147 (MQ=255)
                                                                                                                                      gACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTc            <  7:10577/147‑1 (MQ=255)
                                                                                                                                      gACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTc            <  3:140212/147‑1 (MQ=255)
                                                                                                                                          aTCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGaaa      <  4:173280/149‑1 (MQ=255)
                                                                                                                                          aTCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGaaa      <  7:86297/149‑1 (MQ=255)
                                                                                                                                              gAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGctc  <  4:141058/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CAGACATCCATTTTTTTGAGTCGATAAGAGTAAGGTCGGGCAGGGGTCGAAGACCCACGTCCCGGGCTAACCTGAAGACCGTTTGAGGGCTTATCATGCCAAGACGTCGTGTAGCAGCAAAACGTGAGATCCTTGACGATCCGAAGTACGGATCCCAGATTCTCGCGAAGTTCATGAACCACGTGATGGAAAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTC  >  NC_002947/546555‑546845

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: