Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,520,454 A→T 6.7% 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 6.7% 215.0 / 4.1 74intergenic (+45/+27)cspD/icdDNA replication inhibitor/NADP(+)‑specific isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (37/32);  new base T (3/2);  total (40/34)
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.47e-01

CCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGA  >  NC_002947/4520307‑4520602
                                                                                                                                                        |                                                                                                                                                     
ccccAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTa                                                                                                                                                       >  5:233106/1‑148 (MQ=255)
     aGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGATGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaa                                                                                                                                                 <  7:144566/149‑1 (MQ=255)
      gggCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaag                                                                                                                                                <  4:84966/149‑1 (MQ=255)
          ttGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGTTTATTAAGATGAAAGTACAGACAGCCTGCGCTGTGCCGATAGAAAAAAAAATAGGCCGCATATTTAAaag                                                                                                                                                   >  1:130567/1‑147 (MQ=255)
             cACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCgg                                                                                                                                        >  5:61672/1‑149 (MQ=255)
               cGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGtt                                                                                                                                      <  3:20681/149‑1 (MQ=255)
               cGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGtt                                                                                                                                      >  7:27003/1‑149 (MQ=255)
                gCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGttt                                                                                                                                     >  4:163825/1‑149 (MQ=255)
                    ttGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaag                                                                                                                                                >  4:259225/1‑135 (MQ=255)
                    ttGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaag                                                                                                                                                <  3:259225/135‑1 (MQ=255)
                     tGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGt                                                                                                                                <  1:349672/149‑1 (MQ=255)
                       aaaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCAGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGtt                                                                                                                               <  2:258317/148‑1 (MQ=255)
                        aaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGc                                                                                                                             <  8:27003/149‑1 (MQ=255)
                          tCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCt                                                                                                                           >  5:76926/1‑149 (MQ=255)
                           cAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCt                                                                                                                           >  1:67756/1‑148 (MQ=255)
                            aaaGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCtga                                                                                                                         >  5:343351/1‑147 (MQ=255)
                              aGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACa                                                                                                                       >  6:335201/1‑149 (MQ=255)
                               gATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACAt                                                                                                                      >  6:198846/1‑149 (MQ=255)
                                    cgcgGCCGCGGGGGATGAGGGTGCCTGCGCCTTGGGCTAGCCCAAACCCT‑‑‑A‑‑AACTGAAAACGCGTCAGTCCCTGCGCTGTGCCTTAACACAAAAAAATCGGCCGCTTCTTTAAAACAAGCGACCGGTTTTCTGTTGCCTTACATGTGtt                                                                                                                 <  6:61672/149‑1 (MQ=255)
                                           gCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCAt                                                                                                          <  4:181185/149‑1 (MQ=255)
                                             tGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCAttt                                                                                                        >  6:2910/1‑149 (MQ=255)
                                               cGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTc                                                                                                       <  6:76926/148‑1 (MQ=255)
                                                 gCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATc                                                                                                    >  3:319629/1‑149 (MQ=255)
                                                   tctggCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGGCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGa                                                                                                              <  8:56049/134‑1 (MQ=255)
                                                   agagGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGc                                                                                                  >  7:206097/1‑149 (MQ=255)
                                                   agagGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAAACGGCCGCATCTTTAAAACAAGCGAACGGTTTTCTGTTGCCTTACATGTGtttta                                                                                                              >  7:56049/1‑135 (MQ=255)
                                                       gCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTCTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCaaa                                                                                              >  6:409164/1‑149 (MQ=255)
                                                            cAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGTCGCTTCTTTTAAAGAAGCGGCCGAttttttt                                                                                                                                 <  1:375212/109‑4 (MQ=37)
                                                            cAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGTCGCTTCTTTTAAAGAAGCGGCCGAttttttt                                                                                                                                 >  2:375212/1‑106 (MQ=37)
                                                            cAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCt                                                                                          >  3:290789/1‑148 (MQ=255)
                                                              gagCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCTAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGa                                                                                        <  1:105182/149‑1 (MQ=255)
                                                              gagCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAg                                                                                       >  2:54226/1‑149 (MQ=255)
                                                                  cTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAgctgc                                                                                   >  3:163871/1‑149 (MQ=255)
                                                                   ttGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCaaa                                                                                              >  7:378447/1‑137 (MQ=255)
                                                                   ttGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCaaa                                                                                              <  8:378447/137‑1 (MQ=255)
                                                                     gggCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCtt                                                                                 <  6:233106/148‑1 (MQ=255)
                                                                         gacccGCCAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGCGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAg                                                                            <  4:319629/146‑1 (MQ=255)
                                                                             cGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAGAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTg                                                                        >  8:331144/1‑149 (MQ=255)
                                                                                aaaGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCGTTACATGTGTTTGATAATTTCATCGCCAAAGCCTGAGCTGCTTACACGCTTGgcg                                                                     >  8:355936/1‑149 (MQ=255)
                                                                                aaaGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGgcg                                                                     <  5:335201/149‑1 (MQ=255)
                                                                                aaaGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGgcg                                                                     <  4:163871/149‑1 (MQ=255)
                                                                                    cG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCt                                                                 >  7:202040/1‑149 (MQ=255)
                                                                                     g‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCt                                                                 >  3:242738/1‑148 (MQ=255)
                                                                                             aCTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCAt                                                             <  4:62871/149‑1 (MQ=255)
                                                                                             aCTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCAt                                                             <  4:242738/149‑1 (MQ=255)
                                                                                                    aaGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGc                                                      >  7:154309/1‑149 (MQ=255)
                                                                                                      gccggcctgcctgcTCTGTGCCGATAGACAACAAAAGCGGCCGCTTCTTTAAACGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTc                                                    <  6:343351/143‑1 (MQ=255)
                                                                                                       ccTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCg                                                   <  4:290789/149‑1 (MQ=255)
                                                                                                        cTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTg                                                                        >  2:144154/1‑127 (MQ=255)
                                                                                                               cctgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCg                                             >  1:2552/1‑147 (MQ=255)
                                                                                                               cctgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCg                                             >  3:398029/1‑147 (MQ=255)
                                                                                                                ctgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGt                                            <  2:308059/147‑1 (MQ=255)
                                                                                                                ctgcGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGt                                            >  1:308059/1‑147 (MQ=255)
                                                                                                                      tgtgCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACg                                    >  2:160623/1‑149 (MQ=255)
                                                                                                                            gATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTg                              <  5:183021/149‑1 (MQ=255)
                                                                                                                             aTAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTgg                             >  4:101271/1‑149 (MQ=255)
                                                                                                                              tAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGc                            <  4:79623/149‑1 (MQ=255)
                                                                                                                                 aCAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCg                          >  5:318029/1‑148 (MQ=255)
                                                                                                                                  cAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGc                         >  8:12443/1‑148 (MQ=255)
                                                                                                                                    aaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAAt                      <  1:54226/149‑1 (MQ=255)
                                                                                                                                     aaaaaTCGGTCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATg                     >  8:276896/1‑149 (MQ=255)
                                                                                                                                         aTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGgcgc                 <  8:154309/149‑1 (MQ=255)
                                                                                                                                          tCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGcc                <  6:410708/149‑1 (MQ=255)
                                                                                                                                            ggtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCAt              >  1:30046/4‑149 (MQ=255)
                                                                                                                                             gtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAAt                      >  6:289271/3‑140 (MQ=255)
                                                                                                                                             gtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAAt                      <  5:289271/138‑1 (MQ=255)
                                                                                                                                             gCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCAtt             <  2:67756/149‑1 (MQ=255)
                                                                                                                                              ccGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTg            <  8:202040/149‑1 (MQ=255)
                                                                                                                                              ccGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTg            <  1:160623/149‑1 (MQ=255)
                                                                                                                                                gCTTCTTTAAAAGCAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGt          <  3:373330/149‑1 (MQ=255)
                                                                                                                                                gCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGt          >  8:229975/1‑149 (MQ=255)
                                                                                                                                                 cttcttTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGTCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTgg                             >  6:347846/1‑129 (MQ=255)
                                                                                                                                                 cttcttTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTgg                             <  5:347846/129‑1 (MQ=255)
                                                                                                                                                 cttcttTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCACCATCAATCGCTCGAAATCGTAGGTCAAGGTCTTGGCCGCAATGCAGCCAATGGTg         >  8:253593/1‑149 (MQ=255)
                                                                                                                                                     tttAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCt     >  6:156905/1‑149 (MQ=255)
                                                                                                                                                        aaaaGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTtga  >  7:119438/1‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
CCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCG‑‑‑A‑‑AACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGA  >  NC_002947/4520307‑4520602

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: