Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,585 C→T 7.6% D108D (GAC→GAT) ‡ mcbR → DNA‑binding transcriptional dual regulator McbR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,520,5850CT 7.6% 241.4 / 6.4 79D108D (GAC→GAT) ‡mcbRDNA‑binding transcriptional dual regulator McbR
Reads supporting (aligned to +/- strand):  ref base C (32/41);  new base T (4/2);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGCA  >  NC_000913/1520440‑1520731
                                                                                                                                                 |                                                                                                                                                  
cGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTg                                                                                                                                                <  1:951175/150‑1 (MQ=255)
cGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTg                                                                                                                                                <  1:851462/150‑1 (MQ=255)
 ggTGAATGCGCTTTCTGTCGCACCTGCAGAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGATATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGc                                                                                                                                               <  2:762508/150‑1 (MQ=255)
 ggTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGc                                                                                                                                               >  1:1078899/1‑150 (MQ=255)
 agagAATGAGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGc                                                                                                                                               >  2:408066/4‑150 (MQ=255)
  gTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCg                                                                                                                                              <  2:1026169/150‑1 (MQ=255)
  gTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCg                                                                                                                                              <  1:43486/150‑1 (MQ=255)
   tGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCgg                                                                                                                                             >  2:441673/1‑150 (MQ=255)
     aaTGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGGTCTTGCgggg                                                                                                                                           <  1:1429148/150‑3 (MQ=255)
      aTGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAAc                                                                                                                                          <  2:1294871/150‑1 (MQ=255)
      aTGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAAc                                                                                                                                          <  1:709592/150‑1 (MQ=255)
        gcgcTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTc                                                                                                                                        >  2:42762/1‑150 (MQ=255)
        gcgcTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTc                                                                                                                                        >  1:944758/1‑150 (MQ=255)
            tttCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAgg                                                                                                                                    >  1:767472/1‑150 (MQ=255)
                 gTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTg                                                                                                                               >  2:1138080/1‑150 (MQ=255)
                   ctcACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCt                                                                                                                             <  2:624541/148‑1 (MQ=255)
                    gCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTc                                                                                                                            <  2:759426/150‑1 (MQ=255)
                    gCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTc                                                                                                                            >  2:608894/1‑150 (MQ=255)
                     cACCTTCACAATCATTTACAGTTCCTGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCg                                                                                                                           <  2:610177/150‑1 (MQ=255)
                     cACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCg                                                                                                                           <  2:1396336/150‑1 (MQ=255)
                     cACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCg                                                                                                                           <  1:207002/150‑1 (MQ=255)
                       ccTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCgag                                                                                                                         <  1:1365049/150‑1 (MQ=255)
                              aaGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAc                                                                                                                  <  2:159628/150‑1 (MQ=255)
                               aGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAca                                                                                                                 >  1:219825/1‑150 (MQ=255)
                               aGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAca                                                                                                                 >  1:186246/1‑150 (MQ=255)
                                 cATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAcagc                                                                                                               <  1:408066/150‑1 (MQ=255)
                                                  gTGGGGAAACGTCAATTCGACGAACTCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGc                                                                                              <  1:254165/150‑1 (MQ=255)
                                                  gTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGc                                                                                              <  1:249075/150‑1 (MQ=255)
                                                  gTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGc                                                                                              <  1:1277826/150‑1 (MQ=255)
                                                    ggggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTg                                                                                                                               <  2:144534/115‑1 (MQ=255)
                                                    ggggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTg                                                                                                                               >  1:144534/1‑115 (MQ=255)
                                                    ggggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGa                                                                                            <  2:1405560/150‑1 (MQ=255)
                                                    ggggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGa                                                                                            >  2:508524/1‑150 (MQ=255)
                                                     gggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGAc                                                                                           >  1:248688/1‑150 (MQ=255)
                                                           cGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGaa                                                                                     >  1:635099/1‑150 (MQ=255)
                                                            gTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAAc                                                                                    >  1:366570/1‑150 (MQ=255)
                                                            gTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAAc                                                                                    >  1:1188548/1‑150 (MQ=255)
                                                            gTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAAc                                                                                    <  1:1227090/150‑1 (MQ=255)
                                                               aaTTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACaaa                                                                                 >  2:841108/1‑150 (MQ=255)
                                                                   aGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCAt                                                                             >  1:119725/1‑150 (MQ=255)
                                                                      cGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTaa                                                                          <  1:441673/150‑1 (MQ=255)
                                                                         aatcaatcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGt                                                                       <  2:1368686/150‑1 (MQ=255)
                                                                            caatcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGa                                                                                            <  1:1216403/126‑1 (MQ=255)
                                                                            caatcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGa                                                                                            >  2:1216403/1‑126 (MQ=255)
                                                                             aatcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAAc                                                                   <  2:219825/150‑1 (MQ=255)
                                                                              atcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACa                                                                  >  2:906207/1‑150 (MQ=255)
                                                                              atcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACa                                                                  <  2:782754/150‑1 (MQ=255)
                                                                                cGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAgg                                                                <  1:608894/150‑1 (MQ=255)
                                                                                 ggATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGc                                                               <  2:896044/150‑1 (MQ=255)
                                                                                 ggATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGc                                                               <  2:747182/150‑1 (MQ=255)
                                                                                  gATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCt                                                              >  1:1263057/1‑150 (MQ=255)
                                                                                             aaCTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCCTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAg                                                   >  1:906726/1‑150 (MQ=255)
                                                                                              aCTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGc                                                  <  2:944758/150‑1 (MQ=255)
                                                                                                     ttAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTat                                           <  2:1263057/150‑1 (MQ=255)
                                                                                                      tAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatc                                          <  2:635099/150‑1 (MQ=255)
                                                                                                      tAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatc                                          >  1:542135/1‑150 (MQ=255)
                                                                                                       aaTGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatca                                         >  2:323226/1‑150 (MQ=255)
                                                                                                          ggCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACCGGCTATTTCGCTTAGCGATTTATCATCg                                      >  1:595303/1‑150 (MQ=255)
                                                                                                           gCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGc                                     <  1:252559/150‑1 (MQ=255)
                                                                                                              gTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCa                                  <  2:767472/150‑1 (MQ=255)
                                                                                                                 gCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGc                                                                                              >  1:109841/1‑87 (MQ=255)
                                                                                                                 gCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGc                                                                                              <  2:109841/87‑1 (MQ=255)
                                                                                                                  ctctGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAata                              >  1:539575/1‑150 (MQ=255)
                                                                                                                      ggCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGcc                          <  2:186246/150‑1 (MQ=255)
                                                                                                                      ggCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGcc                          <  1:508524/150‑1 (MQ=255)
                                                                                                                               aaaCCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCtgtg                 >  1:887770/1‑150 (MQ=255)
                                                                                                                               aaaCCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCtgtg                 <  2:542135/150‑1 (MQ=255)
                                                                                                                                  ccTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTga              >  1:544670/1‑150 (MQ=255)
                                                                                                                                   cTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTgag             <  1:1138080/150‑1 (MQ=255)
                                                                                                                                      gttccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgatg          >  1:1214189/4‑150 (MQ=255)
                                                                                                                                      gttccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgatg          >  1:633650/4‑150 (MQ=255)
                                                                                                                                      gttccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgatg          >  1:126274/4‑150 (MQ=255)
                                                                                                                                       ccccGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgatga         >  1:1046590/1‑150 (MQ=255)
                                                                                                                                        cccGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgatgat        >  1:1374061/1‑150 (MQ=255)
                                                                                                                                          cGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTg      <  2:1214189/150‑1 (MQ=255)
                                                                                                                                           gCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGa     >  1:405204/1‑150 (MQ=255)
                                                                                                                                           gCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGATTAGCGATTTATCATCACTCAAATATGCCCATCCTGTGTGAGATGATTGa     >  1:188265/1‑150 (MQ=255)
                                                                                                                                             aaGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGc   <  1:405478/150‑1 (MQ=255)
                                                                                                                                              aGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGCa  <  2:1338268/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGCA  >  NC_000913/1520440‑1520731

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