Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,596,373 A→T 7.3% L213Q (CTG→CAG)  dctD‑III ← C4‑dicarboxylate transport transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,596,3730AT 7.3% 199.4 / 3.6 67L213Q (CTG→CAG) dctD‑IIIC4‑dicarboxylate transport transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (28/34);  new base T (2/3);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAACACGTTC  >  NC_002947/1596231‑1596513
                                                                                                                                                |                                                                                                                                             
tATCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCg                                                                                                                                          <  8:135589/148‑1 (MQ=255)
tATCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTAATGGCCGAACAGTTCgg                                                                                                                                         >  6:306004/1‑149 (MQ=255)
  tCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACTGTTCGGc                                                                                                                                        >  2:176989/1‑148 (MQ=255)
    aCCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACTGTTCGGccat                                                                                                                                     <  1:176989/149‑3 (MQ=255)
    aCCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACTGTTCGGccat                                                                                                                                     <  2:293468/149‑3 (MQ=255)
     cccGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCg                                                                                                                                    <  5:294025/149‑1 (MQ=255)
      ccGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCg                                                                                                                                    >  3:198506/1‑148 (MQ=255)
           aCAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAATGCGCCGCTCTCATGGCCGAACTGTTCGGccat                                                                                                                                     <  5:58932/142‑3 (MQ=255)
           aCAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGACACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACTGTTCGGccat                                                                                                                                     >  6:58932/1‑140 (MQ=255)
              gCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGAc                                                                                                                           >  7:209955/1‑149 (MQ=255)
               cTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACt                                                                                                                          >  6:174912/1‑149 (MQ=255)
                ttGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTc                                                                                                                         >  4:183350/1‑149 (MQ=255)
                    ccTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCgggg                                                                                                                     <  1:290297/149‑1 (MQ=255)
                     cTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGa                                                                                                                    <  3:90912/149‑1 (MQ=255)
                            gCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCcgc                                                                                                             >  5:260407/1‑149 (MQ=255)
                             ccAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATAATCCAGATTGCGGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTAGGGCGTCGAACAGCGACACGGGGATGGCcgcg                                                                                                            >  7:338188/1‑149 (MQ=38)
                                 gggCCTGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGt                                                                                                       <  8:144115/149‑1 (MQ=255)
                                   gCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGtt                                                                                                      >  8:322314/1‑148 (MQ=255)
                                    cATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTaa                                                                                                    <  6:260407/149‑1 (MQ=255)
                                      tGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAgg                                                                                                  >  1:201359/1‑149 (MQ=255)
                                        ctctCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCCACTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAgg                                                                                                  <  1:103005/148‑1 (MQ=255)
                                        ctctCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGc                                                                                                >  1:175057/1‑149 (MQ=255)
                                         tctcGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGc                                                                                                >  7:30455/1‑148 (MQ=255)
                                          ctcGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTa                                                                                              >  8:57441/1‑149 (MQ=255)
                                           tcGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCt                                                                                               <  5:306004/147‑1 (MQ=255)
                                                   gtcgAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCt                                                                                     >  5:228198/1‑149 (MQ=255)
                                                    tcgAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCtt                                                                                    <  3:146026/149‑1 (MQ=255)
                                                      gatgAACTCGGTGGCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTgg                                                                                  <  8:338188/146‑1 (MQ=255)
                                                        ggAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTggcg                                                                                >  6:329114/1‑149 (MQ=255)
                                                         gAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCgggg                                                                                                                     >  6:220532/1‑112 (MQ=255)
                                                         gAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCgggg                                                                                                                     <  5:220532/112‑1 (MQ=255)
                                                          aaCACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTggcggc                                                                              <  3:183350/149‑1 (MQ=255)
                                                                     aCCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCg                                                                    >  7:95773/1‑148 (MQ=255)
                                                                     aCCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCg                                                                    >  5:230495/1‑148 (MQ=255)
                                                                      ccGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGaa                                                                  <  7:57441/149‑1 (MQ=255)
                                                                      ccGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGaa                                                                  <  7:322314/149‑1 (MQ=255)
                                                                            gCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTg                                                            <  2:64852/149‑1 (MQ=255)
                                                                             cATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTgc                                                           <  5:174912/149‑1 (MQ=255)
                                                                             cATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTg                                                            >  6:140690/1‑148 (MQ=255)
                                                                             cATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCGCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTgc                                                           <  2:69968/149‑1 (MQ=255)
                                                                               tACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACTAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTgcag                                                         <  2:201359/149‑1 (MQ=255)
                                                                               tACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTgcag                                                         <  5:329114/149‑1 (MQ=255)
                                                                                  tCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGcc                                                                     >  1:313468/1‑134 (MQ=255)
                                                                                  tCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGcc                                                                     <  2:313468/134‑1 (MQ=255)
                                                                                  tCCAGCTTGCCGATGCGTTTACACTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAAAAGTTAGTCCTCGAACATCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGTCACGGTGGCCGAAGTCGTgcagg                                                        >  3:343172/1‑147 (MQ=255)
                                                                                           ccGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCAc                                             <  8:209955/149‑1 (MQ=255)
                                                                                             gATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCAccc                                           <  3:20299/149‑1 (MQ=255)
                                                                                                  gTTTACCCTGGGCGCCCGTGAAGGCGCCGGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTc                                       >  1:209224/1‑149 (MQ=255)
                                                                                                    ttACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCtt                                    >  7:196744/1‑149 (MQ=255)
                                                                                                       cccTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGc                                  <  8:95773/148‑1 (MQ=255)
                                                                                                       cccTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGc                                  >  1:159519/1‑148 (MQ=255)
                                                                                                       cccTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGc                                  <  2:175057/148‑1 (MQ=255)
                                                                                                       cccTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGc                                  <  2:159519/148‑1 (MQ=255)
                                                                                                          tGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACt                              >  7:61413/1‑149 (MQ=255)
                                                                                                             gcgcCAGTGAATGCGCCGCTCTCATGGCAGAACAGTTCGGCCTCGAACAGCGACTCGGTGATGGCCTCGCAGTTAAGGGCTACAAAGGTCTTGGCGGCACGGGTGCCGAAGTCGTGCAGGCAGCGGGACACCCGCTCCTTGCCACTGcc                           <  4:343172/149‑1 (MQ=255)
                                                                                                                 ccGTGAAGGCGCCGCTCTCATGGCCGATCAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCa                                <  4:348917/140‑1 (MQ=255)
                                                                                                                 ccGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCa                                >  3:348917/1‑140 (MQ=255)
                                                                                                                   gTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGt                                                             <  1:255544/110‑1 (MQ=255)
                                                                                                                   gTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGt                                                             >  2:255544/1‑109 (MQ=255)
                                                                                                                    ggAAGGCGCCGCTCTCATGGCCGAACAGTGCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCa                    <  8:231919/148‑1 (MQ=255)
                                                                                                                     gAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGt                        >  7:174877/1‑144 (MQ=255)
                                                                                                                     gAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGt                        <  8:174877/144‑1 (MQ=255)
                                                                                                                         gcgcCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACcgcg               <  5:332349/149‑1 (MQ=255)
                                                                                                                             cgcTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATc           <  6:230495/149‑1 (MQ=255)
                                                                                                                                 ctctCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCaa         <  8:209360/149‑1 (MQ=255)
                                                                                                                                 ctctCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCaa         <  2:209224/149‑1 (MQ=255)
                                                                                                                                        ggCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAACACGTTc  >  2:158613/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                             
TATCACCCGCAACAGCTTGGCCTGTTGGGCCAGGGGCATGCTCTCGATTTCGTCGAGGAACACGGTGCCACCGTTGGCATACTCCAGCTTGCCGATGCGTTTACCCTGGGCGCCCGTGAAGGCGCCGCTCTCATGGCCGAACAGTTCGGCCTCGAACAGCGACTCGGGGATGGCCGCGCAGTTAAGGGCTACAAAGGGCTTGGCGGCACGGGGGCCGAAGTCGTGCAGGCAGCGGGCCACCCGCTCCTTGCCACTGCCGGTTTCACCGCGGATCAACACGTTC  >  NC_002947/1596231‑1596513

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: