Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,520,440 T→C 7.8% intergenic (+31/+41) 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,4400TC 7.8% 175.2 / 4.8 64intergenic (+31/+41)cspD/icdDNA replication inhibitor/NADP(+)‑specific isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base C (3/2);  total (32/32)
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.72e-01

AACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCG  >  NC_002947/4520293‑4520586
                                                                                                                                                   |                                                                                                                                                  
aaCATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATc                                                                                                                                                   <  1:329387/149‑1 (MQ=255)
aaCATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATc                                                                                                                                                   >  2:53721/1‑149 (MQ=255)
 aCATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCTAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCg                                                                                                                                                  <  1:53721/149‑1 (MQ=255)
  cATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCgg                                                                                                                                                 >  6:187104/1‑149 (MQ=255)
       tCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAACCGGTCGct                                                                                                                                            >  3:186293/1‑149 (MQ=255)
          aGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAACCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGcttgt                                                                                                                                         >  3:72471/1‑147 (MQ=255)
                ccAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGcttct                                                                                                                                         >  8:205104/1‑143 (MQ=255)
                ccAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGcttct                                                                                                                                         <  7:205104/143‑1 (MQ=255)
                ccAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTaa                                                                                                                                     >  4:256270/1‑147 (MQ=255)
                ccAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTaa                                                                                                                                     >  7:131340/1‑147 (MQ=255)
                   aGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaa                                                                                                                                <  1:4799/149‑1 (MQ=255)
                   aGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaa                                                                                                                                <  7:229867/149‑1 (MQ=255)
                    gggCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAACCGGTCGCTTCTTTTAaagaag                                                                                                                               >  3:19958/1‑149 (MQ=255)
                          gCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAACCGGTCGCTTCTTTTAAAGAAGCggc                                                                                                                           <  4:19958/147‑3 (MQ=255)
                             cGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAACGCCTGCCTGCCGGCGCGGTGCCGATAGACAAAAAACTCGGCCGCTTCTTTAAAAGAAGCGACCgggt                                                                                                                      <  2:287755/149‑3 (MQ=255)
                              gCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGACACTGCAAACGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGcttcct                                                                                                                                        >  2:46350/1‑128 (MQ=255)
                                  ttGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAACCGGTCGCTTCTTTTAAAGAAGCGGCCGAttttttt                                                                                                                 <  4:186293/149‑4 (MQ=38)
                                     aaaaccaaGATGCCGCGGCCGCCGCGGCAGAGGCTGCCAGCGCCTTTGGCGAACCGCAAAGCGAAACCCCAAAAGCCTGCCTGCCCGCGCGGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGtt                                                                                                               <  4:72471/142‑1 (MQ=255)
                                     aaaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGtt                                                                                                               <  5:273314/148‑1 (MQ=255)
                                     aaaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGtt                                                                                                               <  3:256270/148‑1 (MQ=255)
                                       aTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCGAGAGCCTTGGGCGAGCCGCTAAGCGAAACTGCAACAGCCTGACTGCCTGCGCTGAGCCGATTGACAAAAAAATCGGCCGCTTCTTTAAATGAAGCGACCTGGTTTCGGTGGcc                                                                                                            >  4:188164/1‑149 (MQ=255)
                                          aaaGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTa                                                                                                         >  6:352464/1‑149 (MQ=255)
                                             gATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACAt                                                                                                      <  3:81812/149‑1 (MQ=255)
                                                gCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACAtgtg                                                                                                   >  6:139114/1‑149 (MQ=255)
                                                  cgcgGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGtt                                                                                                 <  7:329390/149‑1 (MQ=255)
                                                   gcgGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGttt                                                                                                <  5:352464/149‑1 (MQ=255)
                                                   gaggCCGCGGCGGCAGAGGCTGCCAGAGCCTTGGGCGATCCGCAAAGCGAAACTGCAAAAGCGTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGttt                                                                                                <  6:337003/147‑1 (MQ=255)
                                                    cgGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTg                                                                                               <  5:99125/149‑1 (MQ=255)
                                                     ggCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACTAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGa                                                                                              <  5:9918/149‑1 (MQ=255)
                                                          cTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCAtt                                                                                         <  8:131340/149‑1 (MQ=255)
                                                           tGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCg                                                                                                                         <  5:181661/116‑1 (MQ=255)
                                                           tGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCg                                                                                                                         >  6:181661/1‑116 (MQ=255)
                                                              ggCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGACGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCAt                                                                                     >  7:253064/1‑149 (MQ=255)
                                                               gCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATc                                                                                    >  6:49134/1‑149 (MQ=255)
                                                                cAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCg                                                                                   >  5:297128/1‑149 (MQ=255)
                                                                   agGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGcc                                                                                 >  7:307206/1‑148 (MQ=255)
                                                                     gCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCaaa                                                                              <  6:297128/149‑1 (MQ=255)
                                                                          cAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCt                                                                          >  4:103886/1‑148 (MQ=255)
                                                                           agagCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGa                                                                        <  5:139114/149‑1 (MQ=255)
                                                                           agagCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGa                                                                        <  3:103886/149‑1 (MQ=255)
                                                                                cTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAgctgc                                                                   >  5:262872/1‑149 (MQ=255)
                                                                                     gCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTAcc                                                              <  6:262872/149‑1 (MQ=255)
                                                                                      cGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCa                                                             <  4:125956/149‑1 (MQ=255)
                                                                                       gAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAg                                                            >  6:152371/1‑149 (MQ=255)
                                                                                         gCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCt                                                          >  6:337215/1‑149 (MQ=255)
                                                                                                    aaaCTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTcc                                               <  3:265876/149‑1 (MQ=255)
                                                                                                           aaaaGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCAt                                             <  4:341903/144‑1 (MQ=255)
                                                                                                           aaaaGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCAt                                             >  3:341903/1‑144 (MQ=255)
                                                                                                             aaGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGc                                      <  8:307206/149‑1 (MQ=255)
                                                                                                               gCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTc                                    <  5:337215/149‑1 (MQ=255)
                                                                                                                  tGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGa                                  >  5:89989/1‑148 (MQ=255)
                                                                                                                   gACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGaaa                                >  3:96497/1‑149 (MQ=255)
                                                                                                                   gACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGaaa                                >  8:294513/1‑149 (MQ=255)
                                                                                                                                      aTAGACAAAAAAATCGGCAGCTTCTTAAAAAGAAGCGACCGGTTTTCTGTTGCCATACATGTGTTTGATCATATCATCGCCAAAGCCTGAGCTGCTTACCAGCTAGCCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGCTCTTgg             >  4:303622/1‑149 (MQ=255)
                                                                                                                                          aCAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCg          >  2:69380/1‑148 (MQ=255)
                                                                                                                                           acagaaaaCCGGTCGCTTCTTTTAAAGAAGCGGCCGATTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGc         >  5:238676/5‑148 (MQ=38)
                                                                                                                                            aaaaaaaTCGGCCGCTTCTTTAAAATAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTa                           <  2:307809/129‑1 (MQ=255)
                                                                                                                                            aaaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCa                      <  7:165925/134‑1 (MQ=255)
                                                                                                                                            aaaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCa                      >  8:165925/1‑134 (MQ=255)
                                                                                                                                            aaaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCaa       <  4:96497/149‑1 (MQ=255)
                                                                                                                                            aaaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCaa       >  8:215315/1‑149 (MQ=255)
                                                                                                                                             aaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAAt      <  1:69380/149‑1 (MQ=255)
                                                                                                                                               aaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATgg    >  8:201578/1‑149 (MQ=255)
                                                                                                                                                 aaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGgcg  >  7:27704/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCG  >  NC_002947/4520293‑4520586

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: ‑