Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,681,667 A→T 8.4% E62D (GAA→GAT recJ → ssDNA‑specific exonuclease RecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,681,6670AT 8.4% 209.9 / 7.6 71E62D (GAA→GATrecJssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base T (3/3);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTACG  >  NC_002947/1681532‑1681815
                                                                                                                                       |                                                                                                                                                      
ccTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagc                                                                                                                                           >  1:125309/1‑149 (MQ=255)
 cTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGAccagc                                                                                                                                           <  8:35104/148‑1 (MQ=255)
 cTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGAccagc                                                                                                                                           >  7:35104/1‑148 (MQ=255)
   gCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgcc                                                                                                                                        >  4:219086/1‑149 (MQ=255)
    ccACCGTTGCTGACCCGCCTTTACGCCGCGGGGGGGGTGAAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCATACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgcca                                                                                                                                       <  2:381590/149‑1 (MQ=255)
      tcccTTGCTGCCCCGCCTTTACGCCGCGCGCGGCGTGCAGCGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGAccagcaggtcc                                                                                                                                     <  1:366434/148‑7 (MQ=255)
          ttGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGaa                                                                                                                                                        >  4:218556/1‑126 (MQ=255)
          ttGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGaa                                                                                                                                                        <  3:218556/126‑1 (MQ=255)
             cTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCt                                                                                                                             <  6:294418/149‑1 (MQ=255)
                  ccGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGCccagcaggt                                                                                                                                       >  8:292344/1‑131 (MQ=255)
                  ccGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGAccagcaggt                                                                                                                                       <  7:292344/135‑5 (MQ=255)
                  ccGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtcg                                                                                                                        >  2:400168/1‑149 (MQ=255)
                              cgcgcgcgGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTc                                                                                                            >  2:311067/1‑148 (MQ=255)
                               gcgcgcgGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGa                                                                                                          <  3:25499/149‑1 (MQ=255)
                                          gCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTg                                                                                               >  1:3331/1‑149 (MQ=255)
                                            agagCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGcc                                                                                             <  3:219086/149‑1 (MQ=255)
                                            agagCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGc                                                                                              >  7:256224/1‑148 (MQ=255)
                                            agagCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGc                                                                                              >  7:79040/1‑148 (MQ=255)
                                               gCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCAcc                                                                                          >  2:92056/1‑149 (MQ=255)
                                                 gAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGc                                                                                        <  5:404257/149‑1 (MQ=255)
                                                 gAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGAATTCGATGCCGACGGTGCCACCGc                                                                                        >  4:219507/1‑149 (MQ=255)
                                                      cGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCa                                                                                   >  5:346543/1‑149 (MQ=255)
                                                        aGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCAcc                                                                                 <  6:346543/149‑1 (MQ=255)
                                                        aGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCAcc                                                                                 >  1:320243/1‑149 (MQ=255)
                                                            ggACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCg                                                                             >  1:142444/1‑149 (MQ=255)
                                                              aCAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCAcc                                                                                          <  7:20268/134‑1 (MQ=255)
                                                              aCAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCAcc                                                                                          >  8:20268/1‑134 (MQ=255)
                                                              aCAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCggg                                                                           <  2:125309/149‑1 (MQ=255)
                                                              aCAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCggg                                                                           <  8:79040/149‑1 (MQ=255)
                                                               cAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGAAGGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCggg                                                                           >  1:137024/1‑149 (MQ=255)
                                                               cAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCggg                                                                           <  8:256224/148‑1 (MQ=255)
                                                                aaGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGAt                                                                                                         <  5:373742/117‑1 (MQ=255)
                                                                aaGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGAt                                                                                                         >  6:373742/1‑117 (MQ=255)
                                                                    gCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTg                                                                     <  2:320243/149‑1 (MQ=255)
                                                                     ccTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTgg                                                                    <  2:142444/149‑1 (MQ=255)
                                                                      cTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCAc                                                                                           >  1:188706/1‑125 (MQ=255)
                                                                      cTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCAc                                                                                           <  2:188706/125‑1 (MQ=255)
                                                                        ggCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGcc                                                                                       >  6:237116/1‑127 (MQ=255)
                                                                        ggCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGcc                                                                                       <  5:237116/127‑1 (MQ=255)
                                                                             gCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcctgcgc                                                            <  1:311067/149‑1 (MQ=255)
                                                                                tgctgcCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTgg                                                                    >  5:258823/1‑138 (MQ=255)
                                                                                tgctgcCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTgg                                                                    <  6:258823/138‑1 (MQ=255)
                                                                                  ctgcCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACAGTCGGGGTGCTGGgcctgcgcctg                                                         <  7:253363/147‑1 (MQ=255)
                                                                                     cccTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCctgctg                                                      <  1:191037/147‑1 (MQ=255)
                                                                                        tACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGc                                                 >  3:226414/1‑149 (MQ=255)
                                                                                               agcTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCgg                                                                            <  4:3575/115‑1 (MQ=255)
                                                                                               agcTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCgg                                                                            >  3:3575/1‑115 (MQ=255)
                                                                                               agcTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGt                                                   <  4:244990/140‑1 (MQ=255)
                                                                                               agcTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGt                                                   >  3:244990/1‑140 (MQ=255)
                                                                                               agcTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCAt                                          <  1:92056/149‑1 (MQ=255)
                                                                                                gcTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATg                                         >  5:421789/1‑149 (MQ=255)
                                                                                                  tCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCCCCGCCAGCACCGTCggg                                                                           >  3:260809/1‑113 (MQ=255)
                                                                                                  tCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCggg                                                                           <  4:260809/113‑1 (MQ=255)
                                                                                                    aaGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTc                                       <  1:400168/147‑1 (MQ=255)
                                                                                                    aaGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGa                                     >  7:5955/1‑149 (MQ=255)
                                                                                                        gCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTAc                                 >  5:30398/1‑149 (MQ=255)
                                                                                                         cATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTAcc                                >  4:183406/1‑149 (MQ=255)
                                                                                                                gCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGcc                         >  1:176650/1‑149 (MQ=255)
                                                                                                                gCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGcc                         >  1:256061/1‑149 (MQ=255)
                                                                                                                  cGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGt                                        <  4:259140/132‑1 (MQ=255)
                                                                                                                  cGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGt                                        >  3:259140/1‑132 (MQ=255)
                                                                                                                       tgGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGc                  >  7:280200/1‑149 (MQ=255)
                                                                                                                        ggACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCt                 >  6:170168/1‑149 (MQ=255)
                                                                                                                         gACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCtt                >  6:72553/1‑149 (MQ=255)
                                                                                                                          aCCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCCCCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTc               >  2:254041/1‑149 (MQ=255)
                                                                                                                          aCCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTc               <  2:3331/149‑1 (MQ=255)
                                                                                                                            ctgctgGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGa             <  3:183406/149‑1 (MQ=255)
                                                                                                                               ctgGTCGATGGCTTCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTa          >  8:271967/1‑149 (MQ=255)
                                                                                                                                 ggTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATg        >  8:302393/1‑149 (MQ=255)
                                                                                                                                  gTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATg        <  5:170168/148‑1 (MQ=255)
                                                                                                                                  gTCGAAGCCATCCACCAGCGCCAGCGCATCTCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCAATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATg        >  1:245390/1‑149 (MQ=255)
                                                                                                                                     gAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTa    >  7:181230/1‑149 (MQ=255)
                                                                                                                                     gAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTa    >  8:328798/1‑149 (MQ=255)
                                                                                                                                       aGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTACg  <  5:72553/149‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                      
CCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTACG  >  NC_002947/1681532‑1681815

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: