Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,520,464 A→G 11.1% intergenic (+55/+17) cspD → / ← icd DNA replication inhibitor/NADP(+)‑specific isocitrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,520,4640AG11.1% 205.6 / 14.6 82intergenic (+55/+17)cspD/icdDNA replication inhibitor/NADP(+)‑specific isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (34/38);  new base G (3/6);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

TTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATCAGGTC  >  NC_002947/4520317‑4520612
                                                                                                                                                     |                                                                                                                                                    
ttGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAACGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGAc                                                                                                                                                     >  5:261075/1‑149 (MQ=255)
   cACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCg                                                                                                                                                   <  1:386591/148‑1 (MQ=255)
    aCGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGt                                                                                                                                                 >  2:131005/1‑149 (MQ=255)
      gCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGttt                                                                                                                                               >  4:205900/1‑149 (MQ=255)
        cGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAACCGGTCGCTTCTTTTAAAGAAGCGGCCGAttttt                                                                                                                                             <  5:217558/149‑2 (MQ=39)
        cGTTGAAATCAAAGATGCCGCGGCCGCGGCGGCAGAGGCTGACAGAGCCTTGTGCGATCCGCAAAGCGAAACTGCAAAAGCCTGACTTCCTGCGCTGTTCCGAGAGACAA‑‑AAAAATCGGTCGCTTCTTTTAAAGAAGCGGCCGAttttt                                                                                                                                             <  7:371834/149‑2 (MQ=38)
         gTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCt                                                                                                                                            <  7:76246/149‑1 (MQ=255)
          ttGAAATCAAAGATGCCGCCGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAACCGGTCGCTTCTTTTAAAGAAGCGGCCGAttttttt                                                                                                                                           <  5:442833/149‑4 (MQ=38)
             aaaTCAAATATGCAGCGGCCGCTGCGGTTTAGGCTGCAAGAGCCTTGGGCGTGCCGCAAAGCGAAGCTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATTTACAA‑‑AAAAATGGTCCGTTTCTTTAAAAGAAGCGACCGGTTTTCtttt                                                                                                                                         <  4:410951/148‑1 (MQ=255)
             aaaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGtt                                                                                                                                         <  2:74946/148‑1 (MQ=255)
             aaaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGtt                                                                                                                                         <  7:193558/148‑1 (MQ=255)
                tCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCt                                                                                                                                     >  7:250281/1‑149 (MQ=255)
                tCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCt                                                                                                                                     >  5:217131/1‑149 (MQ=255)
                  aaaGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCtt                                                                                                                                    <  3:205900/148‑1 (MQ=255)
                     tctGCCGCAGCGGCCGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACAt                                                                                                                                >  3:379479/3‑149 (MQ=255)
                     gATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTAGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACAt                                                                                                                                <  8:318236/149‑1 (MQ=255)
                          cgcgGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGtt                                                                                                                           >  5:440151/1‑149 (MQ=255)
                             ggCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGa                                                                                                                        <  1:268722/149‑1 (MQ=255)
                                cGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGGTTGATCa                                                                                                                     >  8:380466/1‑149 (MQ=255)
                                    gCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTc                                                                                                                 <  2:308097/149‑1 (MQ=255)
                                     cGGCATAGTCTGCCAGAGCCTTGGGCGAGCCGCAAATCGAAACAGCAAAAGCCTGACTTGCTGCGCTGTGCCGTTATACAACCTAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCAtt                                                                                                                   <  1:131005/148‑1 (MQ=255)
                                       gCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATc                                                                                                              <  8:259752/149‑1 (MQ=255)
                                       gCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATc                                                                                                              >  7:362769/1‑149 (MQ=255)
                                             gCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCaaa                                                                                                        <  8:250281/149‑1 (MQ=255)
                                                 ccAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGc                                                                                                      <  7:306050/147‑1 (MQ=255)
                                                  cAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCAttt                                                                                                                  <  7:113382/134‑1 (MQ=255)
                                                  cAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCAttt                                                                                                                  >  8:113382/1‑134 (MQ=255)
                                                        cTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAgctgc                                                                                             <  8:20373/149‑1 (MQ=255)
                                                           gggCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTa                                                                                          >  5:8186/1‑149 (MQ=255)
                                                                   cGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTg                                                                                  >  6:341016/1‑149 (MQ=255)
                                                                   cGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTg                                                                                  <  7:229265/149‑1 (MQ=255)
                                                                      aaaGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGgcg                                                                               >  3:214942/1‑149 (MQ=255)
                                                                      aaaGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGgcg                                                                               >  1:30293/1‑149 (MQ=255)
                                                                       aaGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTAcc                                                                                        >  8:332694/1‑139 (MQ=255)
                                                                       aaGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTAcc                                                                                        <  7:332694/139‑1 (MQ=255)
                                                                          cGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCt                                                                           <  7:342639/149‑1 (MQ=255)
                                                                           gAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTc                                                                          >  1:186009/1‑149 (MQ=255)
                                                                            gagatccAAAAGCGTGACTGCGTGCGGTGGGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAACCCTGAGCTGCTTACCAGCTTGGCGCCCTcc                                                                         <  3:247522/143‑1 (MQ=255)
                                                                            aaaCTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTc                                                                          >  2:12479/1‑148 (MQ=255)
                                                                             aaCTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCa                                                                        <  4:214942/149‑1 (MQ=255)
                                                                             aaCTGCAAAAGACTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCa                                                                        <  8:362769/149‑1 (MQ=255)
                                                                                tGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCa                                                                        >  1:69466/1‑146 (MQ=255)
                                                                                tGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATc                                                                      >  7:342301/1‑148 (MQ=255)
                                                                                tGCAAAAGCCTGACTGCCTGCGCTGTCCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCa                                                                        <  2:69466/146‑1 (MQ=255)
                                                                                  cAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAAt                                                                   >  1:430981/1‑149 (MQ=255)
                                                                                  cAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCCCCAAAGCCTGAGCTGATTACCAGCTTGGCGGCCTCCaagagc                                                                   >  4:247522/1‑144 (MQ=255)
                                                                                    aaaGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATTGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCg                                                                 <  4:379479/149‑1 (MQ=255)
                                                                                     aaGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAg                                                                                                       >  8:180845/1‑110 (MQ=255)
                                                                                     aaGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAg                                                                                                       <  7:180845/110‑1 (MQ=255)
                                                                                     aaGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCATTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAAt                                                                   >  5:86045/1‑146 (MQ=255)
                                                                                     aaGCCGTACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAAt                                                                   <  6:86045/146‑1 (MQ=255)
                                                                                      aGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCt                                                               >  7:281670/1‑149 (MQ=255)
                                                                                      aGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCt                                                               <  2:30293/149‑1 (MQ=255)
                                                                                        ccTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCg                                                             <  1:75539/149‑1 (MQ=255)
                                                                                           gACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGCCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGATCGaaa                                                          >  3:159798/1‑149 (MQ=255)
                                                                                                 ctgcGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAg                                                    <  1:346457/149‑1 (MQ=255)
                                                                                                    cgcTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTAc                                                                                         >  1:408859/1‑109 (MQ=255)
                                                                                                    cgcTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGt                                                  >  2:99838/1‑148 (MQ=255)
                                                                                                       tgtgCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACg                                              >  4:239825/1‑149 (MQ=255)
                                                                                                             gATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTg                                        <  6:217131/149‑1 (MQ=255)
                                                                                                             gATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTg                                        >  2:353829/1‑149 (MQ=255)
                                                                                                              aTAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTgg                                       >  5:3891/1‑149 (MQ=255)
                                                                                                              aTAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTgg                                       <  1:99838/149‑1 (MQ=255)
                                                                                                                 taCAA‑‑AAAAATGGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGGCACGGTCTTGGcc                                     <  6:261075/147‑1 (MQ=255)
                                                                                                                   cAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGc                                   >  2:126862/1‑148 (MQ=255)
                                                                                                                    aa‑‑aaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCaa                                 >  5:133268/1‑149 (MQ=255)
                                                                                                                    aa‑‑aaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCaa                                 >  4:258086/1‑149 (MQ=255)
                                                                                                                         aaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATgg                              <  1:12479/149‑1 (MQ=255)
                                                                                                                          aaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGgc                             <  2:430981/149‑1 (MQ=255)
                                                                                                                             ccGGTCGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGgcgc                           >  3:62834/2‑148 (MQ=255)
                                                                                                                               ggtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCCAAATAGTAGGTCACGGTCTTGGCCGCAATGGCGACAt                        >  8:150077/4‑149 (MQ=255)
                                                                                                                                  cGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTgg                     <  3:239825/149‑1 (MQ=255)
                                                                                                                                       cttTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGcc                 <  4:62834/148‑1 (MQ=255)
                                                                                                                                        tttAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCt               >  6:124673/1‑149 (MQ=255)
                                                                                                                                          ttaaaGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAg                                                    <  6:427687/108‑1 (MQ=255)
                                                                                                                                          ttaaaGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAg                                                    >  5:427687/3‑110 (MQ=255)
                                                                                                                                          tAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGccc                <  1:126862/146‑1 (MQ=255)
                                                                                                                                            aaaGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTtgat           <  2:131368/149‑1 (MQ=255)
                                                                                                                                            aaaGAAGCGACCGGTTTTCTGTTGCCATACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTtga            <  7:379784/148‑1 (MQ=255)
                                                                                                                                             aagaagCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTtgat           >  8:42313/1‑148 (MQ=255)
                                                                                                                                               gaagCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTgatgat        <  6:8186/149‑1 (MQ=255)
                                                                                                                                                 agCGGCCGGTTTTCTGTTGCCTTACTTGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATc       <  7:150077/148‑1 (MQ=255)
                                                                                                                                                     aCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTAACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATCAGGTc  >  7:36738/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATCAGGTC  >  NC_002947/4520317‑4520612

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: