Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,575 G→C 6.4% intergenic (+65/‑1235) PP_3774 → / → PP_5627 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,5750GC 6.4% 188.9 / 3.9 62intergenic (+65/‑1235)PP_3774/PP_5627hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base C (2/2);  total (33/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 = 6.06e-01

AACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTCTGCC  >  NC_002947/4301427‑4301722
                                                                                                                                                     |                                                                                                                                                   
aaCCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCCATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAgcg                                                                                                                                                     >  2:15092/1‑149 (MQ=255)
  ccAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGcgcg                                                                                                                                                   >  3:17962/1‑149 (MQ=255)
   cAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCgg                                                                                                                                                  >  6:208838/1‑149 (MQ=255)
         gACTACACAGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGC‑GAGCTATCCCTCATCCATCGCTGAACGTCGATCCCCAGCCCTGCTATCAAAACCCCTAAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctgg                                                                                                                                            >  5:16325/1‑149 (MQ=255)
               acCGACGCAACCGCCGCCGATGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCATAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCt                                                                                                                                      <  2:302869/149‑1 (MQ=255)
                ccGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCCCCGCGCTGGGCttttt                                                                                                                                     <  7:93954/149‑4 (MQ=255)
                        aCCCCCGCCGAGGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTgg                                                                                                                             <  7:108299/149‑1 (MQ=255)
                           cccGCCGATGACGGCTCGATGTGCATCATCAGC‑GAGCTCTGGCTCATCCTTCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTTGGCTTCTATCAGCTGGGCt                                                                                                                          <  6:16325/149‑1 (MQ=255)
                            ccgccgAGGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCgg                                                                                                                                                  >  6:137937/1‑124 (MQ=255)
                            ccgccgAGGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCgg                                                                                                                                                  <  5:137937/124‑1 (MQ=255)
                            ccgccgAGGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTg                                                                                                                         <  5:140118/149‑1 (MQ=255)
                            ccgccgAGGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTAACAGCTGGGCt                                                                                                                          >  8:69935/1‑148 (MQ=255)
                                cgAGGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACa                                                                                                                     <  8:11023/149‑1 (MQ=255)
                                  aGGACGGCTCGATCTGCATCATCAGCGGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAg                                                                                                                    >  1:197411/1‑149 (MQ=255)
                                        gCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCTTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTcc                                                                                                             >  6:290965/1‑149 (MQ=255)
                                         cTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGaa                                                                                                                       >  2:380664/1‑138 (MQ=255)
                                         cTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGaa                                                                                                                       <  1:380664/138‑1 (MQ=255)
                                                  catcatCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGc                                                                                                       <  7:7446/145‑1 (MQ=255)
                                                  catcatCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCt                                                                                                     <  4:153457/147‑1 (MQ=255)
                                                       tcaGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTacaca                                                                                              >  1:105842/1‑149 (MQ=255)
                                                       tcaGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTacaca                                                                                              <  5:208838/149‑1 (MQ=255)
                                                               gCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACt                                                                                       <  4:178025/149‑1 (MQ=255)
                                                                   tcGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACt                                                                                       >  6:207495/1‑145 (MQ=255)
                                                                   tcGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACt                                                                                       <  5:207495/145‑1 (MQ=255)
                                                                     gCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGa                                                                                    <  4:17962/146‑1 (MQ=255)
                                                                        catccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGGTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTAc                                                                              >  2:127158/1‑149 (MQ=255)
                                                                         atccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGaa                                                                                   <  4:295496/143‑1 (MQ=255)
                                                                         atccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGaa                                                                                   >  3:295496/1‑143 (MQ=255)
                                                                          tccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCACTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGaa                                                                                                                       >  7:162892/1‑106 (MQ=255)
                                                                           ccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTc                                                                            >  5:102488/1‑148 (MQ=255)
                                                                            catcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTc                                                                            >  4:377503/1‑147 (MQ=255)
                                                                                gCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTgg                                                                      <  4:3673/149‑1 (MQ=255)
                                                                                 cTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctgggc                                                                                                                                          <  8:154414/80‑1 (MQ=255)
                                                                                 cTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctgggc                                                                                                                                          >  7:154414/1‑80 (MQ=255)
                                                                                        tCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTtgt                                                              <  1:15092/149‑1 (MQ=255)
                                                                                          gATGCCCAGCCCTGCTATCAAAACCCCTCACGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGa                                                            >  1:31852/1‑149 (MQ=255)
                                                                                              cccAGCCTTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAAc                                                        <  5:290965/149‑1 (MQ=255)
                                                                                               ccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTTTCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAAc                                                        >  2:98030/1‑148 (MQ=255)
                                                                                               ccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAAc                                                        >  5:125191/1‑148 (MQ=255)
                                                                                                   cccTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATaa                                                   <  2:105842/149‑1 (MQ=255)
                                                                                                      tGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACAATAACGc                                                >  3:151030/1‑149 (MQ=255)
                                                                                                        cTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTc                                              >  3:68534/1‑149 (MQ=255)
                                                                                                        cTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGAGGCTGGGCTTCTATCAGCTGGGCTGAACAGTAATTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTc                                              >  7:238684/1‑149 (MQ=255)
                                                                                                          aTCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTtg                                                               >  6:261252/1‑130 (MQ=255)
                                                                                                          aTCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTtg                                                               <  5:261252/130‑1 (MQ=255)
                                                                                                           tCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCgg                                            >  8:311768/1‑148 (MQ=255)
                                                                                                              tatccccTAAGGGCGTCTCATGTATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCAtt                                        <  7:69935/146‑1 (MQ=255)
                                                                                                               atccccTCAGGGCTTCGCATGATTGCAAAAAGCGCAGCGGGGGGCTGGGCTTCTATCAGCTGGGCTGAACa                                                                                                                     <  6:320088/69‑1 (MQ=255)
                                                                                                               aaCCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTTGCCTACTATCAGCAGGGCTGAACa                                                                                                                     >  5:320088/1‑71 (MQ=37)
                                                                                                               aaCCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTg                                       >  6:33998/1‑149 (MQ=255)
                                                                                                                    cTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGc                                                                                                                           >  3:182733/1‑60 (MQ=255)
                                                                                                                    cTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGc                                                                                                                           <  4:182733/60‑1 (MQ=255)
                                                                                                                         ggCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCcgc                              >  1:186511/1‑148 (MQ=255)
                                                                                                                              cGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATggg                        >  6:108108/1‑149 (MQ=255)
                                                                                                                                 aTGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGAt                                    >  6:277516/1‑134 (MQ=255)
                                                                                                                                 aTGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGAt                                    <  5:277516/134‑1 (MQ=255)
                                                                                                                                   gAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGCGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCcgg                   >  6:265597/1‑149 (MQ=255)
                                                                                                                                      tgatagaaGCCCAGCCCCGCGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCcggcgg                >  5:141725/7‑149 (MQ=255)
                                                                                                                                       gatagaaGCCCAGCCCCGCGCTGGGCTTTTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCggg               >  1:223833/6‑149 (MQ=255)
                                                                                                                                        cAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGc              >  1:295626/1‑149 (MQ=255)
                                                                                                                                         aaaaaGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACAATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGcc             <  4:151030/149‑1 (MQ=255)
                                                                                                                                           aaaGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTg           <  5:33998/149‑1 (MQ=255)
                                                                                                                                            aaGCCCAGCCCCGCGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGa          <  6:141725/149‑1 (MQ=255)
                                                                                                                                                  aGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTctg    <  6:125191/149‑1 (MQ=255)
                                                                                                                                                   gcgcGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTctgc   <  6:102488/149‑1 (MQ=255)
                                                                                                                                                    cgcgGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTctgcc  <  2:186511/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AACCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGC‑GAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTCTGCC  >  NC_002947/4301427‑4301722

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