Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,577 G→C 6.7% intergenic (+67/‑1233) 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,5770GC 6.7% 187.5 / 3.9 60intergenic (+67/‑1233)PP_3774/PP_5627hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base C (2/2);  total (31/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.27e-01

CCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTCTGCC  >  NC_002947/4301429‑4301722
                                                                                                                                                     |                                                                                                                                                 
ccAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGcgcg                                                                                                                                                   >  3:17962/1‑149 (MQ=255)
 cAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCgg                                                                                                                                                  >  6:208838/1‑149 (MQ=255)
       gACTACACAGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTATCCCTCATCCATCGCTGAACGTCGATCCCCAGCCCTGCTATCAAAACCCCTAAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctgg                                                                                                                                            >  5:16325/1‑149 (MQ=255)
             acCGACGCAACCGCCGCCGATGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCATAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCt                                                                                                                                      <  2:302869/149‑1 (MQ=255)
              ccGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCCCCGCGCTGGGCttttt                                                                                                                                     <  7:93954/149‑4 (MQ=255)
                      aCCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTgg                                                                                                                             <  7:108299/149‑1 (MQ=255)
                         cccGCCGATGACGGCTCGATGTGCATCATCAGCGAGCTCTGGCTCATCCTTCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTTGGCTTCTATCAGCTGGGCt                                                                                                                          <  6:16325/149‑1 (MQ=255)
                          ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCgg                                                                                                                                                  >  6:137937/1‑124 (MQ=255)
                          ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCgg                                                                                                                                                  <  5:137937/124‑1 (MQ=255)
                          ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTg                                                                                                                         <  5:140118/149‑1 (MQ=255)
                          ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTAACAGCTGGGCt                                                                                                                          >  8:69935/1‑148 (MQ=255)
                              cgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACa                                                                                                                     <  8:11023/149‑1 (MQ=255)
                                aGGACGGCTCGATCTGCATCATCAGCGGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAg                                                                                                                    >  1:197411/1‑149 (MQ=255)
                                      gCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCTTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTcc                                                                                                             >  6:290965/1‑149 (MQ=255)
                                       cTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGaa                                                                                                                       >  2:380664/1‑138 (MQ=255)
                                       cTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGaa                                                                                                                       <  1:380664/138‑1 (MQ=255)
                                                catcatCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGc                                                                                                       <  7:7446/145‑1 (MQ=255)
                                                catcatCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCt                                                                                                     <  4:153457/147‑1 (MQ=255)
                                                     tcaGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTacaca                                                                                              >  1:105842/1‑149 (MQ=255)
                                                     tcaGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTacaca                                                                                              <  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)
                                                                                                                                                     |                                                                                                                                                 
CCAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTCTGCC  >  NC_002947/4301429‑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: