Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,865,939 Δ1 bp 9.8% coding (983/1020 nt) PP_3415 ← LacI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,865,9390T. 9.8% 174.5 / 14.3 61coding (983/1020 nt)PP_3415LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (32/23);  new base . (3/3);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

CATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGAC  >  NC_002947/3865802‑3866084
                                                                                                                                         |                                                                                                                                                               
cATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGATAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCAt                                                                                                                                                     >  1:1405795/1‑150 (MQ=255)
              ccTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTc                                                                                                                                       >  1:1176369/1‑150 (MQ=255)
                     cGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCa                                                                                                                       >  1:762182/1‑150 (MQ=255)
                       cAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCtt                                                                                                                                        <  1:1190904/140‑1 (MQ=255)
                       cAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCtt                                                                                                                                        >  2:1190904/1‑140 (MQ=255)
                                  gCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCCCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCtt                                                                                                                                        >  2:756874/1‑129 (MQ=255)
                                  gCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCCCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCtt                                                                                                                                        <  1:756874/129‑1 (MQ=255)
                                   cTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGACCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAgcg                                                                                                        <  2:210676/150‑1 (MQ=255)
                                      cACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGc‑‑‑‑‑c                                                                                                 >  2:106988/1‑150 (MQ=255)
                                          ttCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑Ctgtg                                                                                             >  1:720178/1‑150 (MQ=255)
                                                 cGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGcgcc                                                                                      >  2:478551/1‑150 (MQ=255)
                                                    gCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGACCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGcgccccgc                                                                                  >  1:623567/1‑150 (MQ=255)
                                                     ccAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGcgccccgc                                                                                  >  2:878581/1‑150 (MQ=255)
                                                       aaTCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGcgccccgccc                                                                                >  1:929954/1‑150 (MQ=255)
                                                        aTCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCCTCAGCTCAAAGCCCAGATCCTCCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAACAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCgccccgcccc                                                                               >  1:192683/1‑150 (MQ=255)
                                                         tCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCca                                                                                                              >  1:70647/1‑123 (MQ=255)
                                                         tCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCca                                                                                                              <  2:70647/123‑1 (MQ=255)
                                                         tCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAACAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCa                                                                              >  1:637691/1‑150 (MQ=255)
                                                             tCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGc                                                                          >  1:225965/1‑150 (MQ=255)
                                                               cTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCctgc                                                                                                                                                            <  2:94747/80‑1 (MQ=255)
                                                               cTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCctgc                                                                                                                                                            >  1:94747/1‑80 (MQ=255)
                                                                tgtgGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGc                                                                       >  1:1091672/1‑150 (MQ=255)
                                                                     gAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAggaggtc                                                                  <  2:225965/150‑7 (MQ=255)
                                                                          gcgcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCAccc                                                                                                                                            >  1:343345/1‑85 (MQ=255)
                                                                          gcgcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCAccc                                                                                                                                            <  2:343345/85‑1 (MQ=255)
                                                                            gcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGAt                                                                                                                                                                 >  2:941834/1‑62 (MQ=255)
                                                                            gcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGAt                                                                                                                                                                 <  1:941834/62‑1 (MQ=255)
                                                                              gagaAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGaa                                                         <  2:71219/150‑1 (MQ=255)
                                                                                  aGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGAACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGa                                                     <  2:1366968/150‑1 (MQ=255)
                                                                                   ggCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAg                                                    <  1:106988/150‑1 (MQ=255)
                                                                                    gCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGcc                                                                                           <  1:171869/110‑1 (MQ=255)
                                                                                    gCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGcc                                                                                           >  2:171869/1‑110 (MQ=255)
                                                                                     ccTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCcgc                                                                                         <  1:238076/111‑1 (MQ=255)
                                                                                     ccTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCcgc                                                                                         >  2:238076/1‑111 (MQ=255)
                                                                                      cTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGa                                                          >  1:1140522/1‑150 (MQ=255)
                                                                                       ttGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCa                                                >  1:131466/1‑150 (MQ=255)
                                                                                          gccggCATTCTCAGCCGCGCTCGCGCACCATCAGCTCAAAGTCACACTACTGAGAAGCACGTGACTGGCGTTTCCCAC‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCCCCGCAGCCCGCTGTGTGCGAATGGAGGTCAGcc                                             <  1:878581/147‑1 (MQ=255)
                                                                                                    tcAGGTGCTCTCGCGCACAATCAGCTGACAGAAAAGATGATGCTGCGCATTTGTCACCGGGTTTCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGCCTGTGCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGc                                        >  2:1231707/1‑150 (MQ=255)
                                                                                                        gTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATgc                               <  2:720178/150‑1 (MQ=255)
                                                                                                          gCTCTCGCGCTCCTTCAGCGCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATgcgc                             >  2:1308283/1‑150 (MQ=255)
                                                                                                          gCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAc                                      <  2:1140522/150‑1 (MQ=255)
                                                                                                             ctcGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGa                          >  1:818417/1‑150 (MQ=255)
                                                                                                                   cACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                    <  2:1405795/150‑1 (MQ=255)
                                                                                                                   cACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAACAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                    <  2:637691/150‑1 (MQ=255)
                                                                                                                    aCCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCggg                   >  2:743438/1‑150 (MQ=255)
                                                                                                                      cctCAGCTCAAAGCCCGGATCCTGCGGCGCATTCGGCACGCGCTTTCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCAGCGACCCGACCCACCACACCACTGAAGACGCGATACAATGGCAGCTGGTGCACCATATGCGCCGAGGCCGGgag                 >  2:312364/3‑150 (MQ=255)
                                                                                                                        tCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGAt               >  2:480032/1‑150 (MQ=255)
                                                                                                                        tCAGCTCAAAGCCCAGACCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                    >  2:789230/1‑144 (MQ=255)
                                                                                                                        tCAGCTCAAAGCCCAGACCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                    <  1:789230/144‑1 (MQ=255)
                                                                                                                        ccAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGACTGGACGACAGCCGTGGCACCATATGCGCCGAGGCCGGGAGAt               <  1:743438/149‑1 (MQ=255)
                                                                                                                            cTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCCGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt           <  1:178402/150‑1 (MQ=255)
                                                                                                                            cTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCCCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCTCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt           <  2:192683/150‑1 (MQ=255)
                                                                                                                             tCAAAGCCCAGACCTGCTGCGCATTCGTCACCCGCTTCTCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCg             >  2:484513/1‑146 (MQ=255)
                                                                                                                             tCAAAGCCCAGACCTGCTGCGCATTCGTCACCCGCTTCTCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCg             <  1:484513/146‑1 (MQ=255)
                                                                                                                              cAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGa         >  1:966583/1‑150 (MQ=255)
                                                                                                                               aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGCCTGTGCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCg             >  1:742737/1‑150 (MQ=255)
                                                                                                                                aaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGaaa       <  1:478551/150‑1 (MQ=255)
                                                                                                                                  gCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt           >  2:22870/1‑144 (MQ=255)
                                                                                                                                  gCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt           <  1:22870/144‑1 (MQ=255)
                                                                                                                                   cccAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGc                                                                       <  2:196372/83‑1 (MQ=255)
                                                                                                                                   cccAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGc                                                                       >  1:196372/1‑83 (MQ=255)
                                                                                                                                    ccAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGa   <  2:1091672/150‑1 (MQ=255)
                                                                                                                                     cAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGAc  >  1:1004971/1‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                                               
CATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCAT‑‑‑‑‑‑‑‑‑CCAGCAACGCCAGCAGCGC‑‑‑‑‑CTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGAC  >  NC_002947/3865802‑3866084

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