Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,520,454 A→T 7.4% intergenic (+45/+27) 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,4540AT 7.4% 214.4 / 5.4 77intergenic (+45/+27)cspD/icdDNA replication inhibitor/NADP(+)‑specific isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (34/37);  new base T (2/4);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCT  >  NC_002947/4520308‑4520599
                                                                                                                                                     |                                                                                                                                                 
cccAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTaaa                                                                                                                                                 <  8:17685/149‑1 (MQ=255)
 ccAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGGTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAGACTGCAAGAGCCTGGCTGCCTGAGCTGTGCCGCTAGACAA‑‑TAAAATCGGCCGCTTCTTTaa                                                                                                                                                  >  8:379784/1‑147 (MQ=255)
 ccAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTaa                                                                                                                                                  >  7:20373/1‑147 (MQ=255)
 ccAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTaa                                                                                                                                                  >  2:386591/1‑147 (MQ=255)
  cAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGAGCATGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATAGGCCGCTTCTTTaaa                                                                                                                                                 >  1:131368/1‑148 (MQ=255)
  cAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTaa                                                                                                                                                  >  7:259752/1‑146 (MQ=255)
   aaGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGTTAGACAA‑‑AAAAATCGGCCGCTTCTTTAAaaga                                                                                                                                              <  1:52103/149‑1 (MQ=255)
      ggCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGc                                                                                                                                           <  7:251823/149‑1 (MQ=255)
         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                                                                                                                            <  7:193558/148‑1 (MQ=255)
                      aaaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGtt                                                                                                                            <  2:74946/148‑1 (MQ=255)
                         tCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCt                                                                                                                        >  5:217131/1‑149 (MQ=255)
                         tCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCt                                                                                                                        >  7:250281/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                                                                     <  7:229265/149‑1 (MQ=255)
                                                                             cGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTg                                                                     >  6:341016/1‑149 (MQ=255)
                                                                                aaaGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGgcg                                                                  >  1:30293/1‑149 (MQ=255)
                                                                                aaaGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGgcg                                                                  >  3:214942/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)
                                                                                                                                                    tAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGccc   <  1:126862/146‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAA‑‑AAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCT  >  NC_002947/4520308‑4520599

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: