Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,080,989 T→A 14.9% intergenic (+36/+26) PP_5332 → / ← PP_5333 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,080,9890TA14.9% 151.7 / 17.5 66intergenic (+36/+26)PP_5332/PP_5333hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base A (6/4);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

TGGCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGG  >  NC_002947/6080843‑6081129
                                                                                                                                                  |                                                                                                                                            
tggCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGc                                                                                                                                            <  7:285039/149‑1 (MQ=255)
   cAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCcg                                                                                                                                         >  3:179264/1‑149 (MQ=255)
    aGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                        >  7:229365/1‑149 (MQ=255)
     gCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                        >  3:176601/1‑148 (MQ=255)
         aTGCTGGTCATCTAGAGAGCGTATTTGTTGGCGGAGCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATATGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGTCCTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                   <  2:91778/149‑1 (MQ=255)
          tGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAagag                                                                                                                                  <  2:363657/149‑1 (MQ=255)
                  aTCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                              <  6:206187/145‑1 (MQ=255)
                  aTCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                              >  5:206187/1‑145 (MQ=255)
                        aCAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAAtt                                                                                                                      >  3:298595/1‑147 (MQ=255)
                        aCAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAAtt                                                                                                                      <  4:298595/147‑1 (MQ=255)
                          aGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATc                                                                                                                   >  5:366503/1‑148 (MQ=255)
                          aGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATc                                                                                                                   >  5:356749/1‑148 (MQ=255)
                           gCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCa                                                                                                                  >  1:333251/1‑148 (MQ=255)
                             gAATGTGTGGGATGACCTGGTGGAGTTTAAGCCTGGGGCGAGTCAATCGCTGGGAGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTATACCGGCCTCTTCGCGGGCTTGACCGCGAAGAGGCAGCGGCAATTATCAAAg                                                                                                               >  7:175975/1‑149 (MQ=255)
                             gAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAg                                                                                                               >  4:259386/1‑149 (MQ=255)
                             gAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAg                                                                                                               <  1:13120/149‑1 (MQ=255)
                               ttgtgggGGTTGACCTGGTGGTGTTCATCCCTGGGGCTTGTCATTCGCTGTGCGTCCAGATCGGTTGTGAAAATATCTCAGCCCTCATCTTTTTTTACCGGCCTCTCCGCGGTCTGGCCCTCGAATAGGCCGGGGCAATTATCAAAg                                                                                                               <  4:316248/146‑1 (MQ=255)
                               aTGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGTCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGt                                                                                                             <  1:78489/149‑1 (MQ=255)
                                    tgGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCg                                                                                                        <  6:262297/149‑1 (MQ=255)
                                     gggATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCACGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGt                                                                                                       >  3:225542/1‑149 (MQ=255)
                                           aCCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGGCTGCGTCTTcc                                                                                                  <  6:366503/148‑1 (MQ=255)
                                            ccTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTccac                                                                                                <  3:259386/149‑1 (MQ=255)
                                                 tggAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCAccc                                                                                           >  6:352059/1‑149 (MQ=255)
                                                     gTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGGGGGGTTTACCGGGGATGGGGCGGGGCAATTATCAAATGGGTGCGTGTTGCCCCACCCTCAc                                                                                       >  3:17959/1‑149 (MQ=255)
                                                     gTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGcc                                                                                                                                           >  4:114469/1‑97 (MQ=255)
                                                     gTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGcc                                                                                                                                           <  3:114469/97‑1 (MQ=255)
                                                                 ggCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAg                                                                                                               <  6:153428/113‑1 (MQ=255)
                                                                 ggCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAg                                                                                                               >  5:153428/1‑113 (MQ=255)
                                                                         aTTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCa                                                                    >  7:238480/1‑148 (MQ=255)
                                                                         aTTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCa                                                                    >  1:330396/1‑148 (MQ=255)
                                                                          ttCGCTGGGCGTCCATATCTGTTTTTAAATTTTTTTAGCCCTGCTCTTTTTTTACGGGCCTCTTCGCGGGCTTGCCCGCTAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATa                                                                  <  2:322569/149‑1 (MQ=255)
                                                                             gCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                     <  1:24345/79‑1 (MQ=255)
                                                                             gCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                     >  2:24345/1‑79 (MQ=255)
                                                                                   cGTCCAGATCGGTTGTGAAAATCTCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATc                                                         <  4:225542/149‑1 (MQ=255)
                                                                                        aGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGc                                                      >  7:9633/1‑147 (MQ=255)
                                                                                              gTTGTGAAAATATCTGAGCCCTGCTCATTTTTTACCGGCCTCTTCGCGGGCTAGCCCGCGAAGAGGCCGGTGCAATTATCAAAGGTCTGCGTCTGCAACAACCCTCCCCTTGTCTGCAAACAGCGCATAGGCCGCATAGGCCAGGTCAt                                              >  4:294856/1‑149 (MQ=255)
                                                                                               ttGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCAtt                                             <  5:352059/149‑1 (MQ=255)
                                                                                                     aaaTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGAc                                       <  6:131780/149‑1 (MQ=255)
                                                                                                           cTGAGCCCTGCTCTTTTTTTATCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAgc                                                                       <  4:75918/111‑1 (MQ=255)
                                                                                                           cTGAGCCCTGCTCTTTTTTTATCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCAACCTCACCTTGTCTGCATCCAgc                                                                       >  3:75918/1‑111 (MQ=255)
                                                                                                             gAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCAt                                                                   <  6:109901/113‑1 (MQ=255)
                                                                                                              aGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCAc                              <  2:330396/149‑1 (MQ=255)
                                                                                                                        tttttttACCGGCCTCTTCGCTGGCTTGCCCGCGAAGAGGCCGGTGCAATTATCAAATGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCTGATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCt                      <  4:137248/147‑1 (MQ=255)
                                                                                                                        tttttttACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTg                    >  4:356964/1‑149 (MQ=255)
                                                                                                                         aattgccccGGCCTCTTCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGc                   >  6:82471/8‑149 (MQ=255)
                                                                                                                          tttttACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGc                   <  5:71138/148‑1 (MQ=255)
                                                                                                                           ttgccccGGCCTCTTCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCg                 <  5:82471/144‑1 (MQ=255)
                                                                                                                                 cGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCAc                                                                                       <  6:344584/73‑1 (MQ=255)
                                                                                                                                 cGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCAc                                                                                       >  5:344584/1‑73 (MQ=255)
                                                                                                                                 cGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGAcc           >  7:276789/1‑149 (MQ=255)
                                                                                                                                    ccTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCcac        <  6:356749/149‑1 (MQ=255)
                                                                                                                                    ccTCTTCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGa             <  8:155188/144‑1 (MQ=255)
                                                                                                                                    ccTCTTCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGAcagat        >  7:155188/1‑145 (MQ=255)
                                                                                                                                     ctctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCcaca       <  4:179264/149‑1 (MQ=255)
                                                                                                                                      tctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCa                         >  6:97893/1‑130 (MQ=255)
                                                                                                                                      tctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCa                         <  5:97893/130‑1 (MQ=255)
                                                                                                                                      tctTCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAg      >  5:127382/1‑149 (MQ=255)
                                                                                                                                        ttCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCg    >  6:286432/1‑149 (MQ=255)
                                                                                                                                        ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCAtt                                             <  4:155249/108‑1 (MQ=255)
                                                                                                                                        ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCAtt                                             >  3:155249/1‑108 (MQ=255)
                                                                                                                                        ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCg    >  3:135903/1‑149 (MQ=255)
                                                                                                                                          ggcgGCCTGGCCCGCTAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCggg  <  2:333251/148‑1 (MQ=255)
                                                                                                                                          cgcgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCgg   <  8:238480/148‑1 (MQ=255)
                                                                                                                                           gcgGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCggg  <  6:127382/148‑1 (MQ=255)
                                                                                                                                               gCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTg                                         <  8:199621/105‑1 (MQ=255)
                                                                                                                                               gCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTg                                         >  7:199621/1‑105 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
TGGCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGG  >  NC_002947/6080843‑6081129

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