Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,577,815 G→C 5.9% intergenic (+1/‑44) PP_4906 → / → rhdA hypothetical protein/thiosulfate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,577,8150GC 5.9% 223.2 / 3.7 68intergenic (+1/‑44)PP_4906/rhdAhypothetical protein/thiosulfate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (27/37);  new base C (2/2);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

TCGCTGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGC  >  NC_002947/5577668‑5577956
                                                                                                                                                      |                                                                                                                                             
tcgctGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGc                                                                                                                                              >  1:298769/1‑149 (MQ=255)
   ctGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCtt                                                                                                                                           <  8:69343/149‑1 (MQ=255)
      gCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTg                                                                                                                                          <  2:298769/147‑1 (MQ=255)
       cAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGCTCAACCCTTgcgg                                                                                                                                       <  5:207103/149‑2 (MQ=255)
          ccAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTa                                                                                                                                     >  1:218820/1‑148 (MQ=255)
          ccAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTa                                                                                                                                     >  3:80666/1‑148 (MQ=255)
           cAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTAc                                                                                                                                    >  3:34929/1‑148 (MQ=255)
           cAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTAc                                                                                                                                    >  6:79847/1‑148 (MQ=255)
             aCGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAATCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAAt                                                                                                                                 <  2:304273/149‑1 (MQ=255)
                  cAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGccc                                                                                                                            <  6:148465/149‑1 (MQ=255)
                   aGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGccc                                                                                                                            <  1:88883/149‑1 (MQ=255)
                     gTGCTGATGCTGCTGGTGGCCGACCATGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCCGCAAGGGTTGAGCCTTGCGCTACAATCGc                                                                                                                              <  1:196634/145‑1 (MQ=255)
                     gTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGc                                                                                                                              >  2:196634/1‑144 (MQ=255)
                        cTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCttt                                                                                                                      <  7:188139/149‑1 (MQ=255)
                           aTGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTc                                                                                                                     <  6:311843/147‑1 (MQ=255)
                           aTGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTc                                                                                                                     <  1:252480/147‑1 (MQ=255)
                                    gtggCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGga                                                                                                            <  1:246174/148‑1 (MQ=255)
                                    gtggCCGACCAGGGCTGCCGGCCTTTGGACGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGTTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCTCCCCTCTTTCCAAACTGgaga                                                                                                          >  6:78947/1‑149 (MQ=255)
                                     tggCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGAc                                                                                                         >  5:328002/1‑149 (MQ=255)
                                       gCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGCTCAACCCTTgcggttggt                                                                                                                                  <  2:274503/122‑7 (MQ=255)
                                       gCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGCTCAACCCTTgcggttggt                                                                                                                                  >  1:274503/1‑116 (MQ=255)
                                         cGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTg                                                                                                     <  3:303875/149‑1 (MQ=255)
                                         cGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTg                                                                                                     >  8:215857/1‑149 (MQ=255)
                                         cGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTg                                                                                                     >  4:59440/1‑149 (MQ=255)
                                          gACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTgg                                                                                                    <  6:328002/149‑1 (MQ=255)
                                          gACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTgg                                                                                                    <  6:295722/149‑1 (MQ=255)
                                          gACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTgg                                                                                                    <  4:70982/149‑1 (MQ=255)
                                               gggCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGAc                                                                                               >  4:265164/1‑149 (MQ=255)
                                                ggCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGAcc                                                                                              <  8:73232/149‑1 (MQ=255)
                                                ggCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGAcc                                                                                              <  4:34929/149‑1 (MQ=255)
                                                   tGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGAc                                                                                           <  4:217385/149‑1 (MQ=255)
                                                     ccGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGAc                                                                                           <  2:218820/147‑1 (MQ=255)
                                                       ggCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACtttt                                                                                       >  3:267989/1‑149 (MQ=255)
                                                        gCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTc                                                                                      <  5:79847/149‑1 (MQ=255)
                                                         ccTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTcc                                                                                     <  5:141877/149‑1 (MQ=255)
                                                         ccTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTcc                                                                                     <  3:280847/149‑1 (MQ=255)
                                                                ccGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGAATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCAGACTCGAGACCGTGGATGACCGACTATTCCGGCCTGc                                                                              >  8:251234/1‑149 (MQ=255)
                                                                      tCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTgg                                                                        <  7:105667/149‑1 (MQ=255)
                                                                      tCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTgg                                                                        <  3:59440/149‑1 (MQ=255)
                                                                        aGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTg                                                                      <  5:1752/149‑1 (MQ=255)
                                                                               aaGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAg                                                               >  1:172273/1‑149 (MQ=255)
                                                                                  ccTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTg                                                            <  4:267989/149‑1 (MQ=255)
                                                                                   cTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGAt                                                                    >  6:166559/1‑140 (MQ=255)
                                                                                   cTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGAt                                                                    <  5:166559/140‑1 (MQ=255)
                                                                                           aGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACctgc                                                   <  2:172273/149‑1 (MQ=255)
                                                                                            gACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACctgct                                                  >  7:16005/1‑149 (MQ=255)
                                                                                              ccGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCtgctgc                                                <  7:256091/149‑1 (MQ=255)
                                                                                               cgcgCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGc                                                                                  >  7:266074/1‑114 (MQ=255)
                                                                                               cgcgCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGc                                                                                  <  8:266074/114‑1 (MQ=255)
                                                                                                     aTGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGCCCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCt                                         >  7:206283/1‑149 (MQ=255)
                                                                                                         aTCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGAt                                     >  5:164546/1‑149 (MQ=255)
                                                                                                             aGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTcgc                                 >  5:36451/1‑149 (MQ=255)
                                                                                                                  ggCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTcaccgc                              >  8:244567/1‑148 (MQ=255)
                                                                                                                         tCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCACCGCAACTGATc                      <  7:244567/149‑1 (MQ=255)
                                                                                                                           ggTGTTTACCAACCGCAATGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATcc                     <  4:80666/148‑1 (MQ=255)
                                                                                                                            gtgtTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTg                   >  5:192299/1‑149 (MQ=255)
                                                                                                                               tttACCAACCGGAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGAGCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTc                <  7:251234/149‑1 (MQ=255)
                                                                                                                               tttACCAACCGCAAGGGTTGAGCCTTTCGCTACAATCGCCCGTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGt                 <  5:78947/148‑1 (MQ=255)
                                                                                                                               tttACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGt                 <  6:192299/148‑1 (MQ=255)
                                                                                                                                 tACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGa              >  7:174829/1‑149 (MQ=255)
                                                                                                                                  aCCAACCGCAAGGGTTGCGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGAc             >  8:236142/1‑149 (MQ=255)
                                                                                                                                  aCCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGAt                                                                    <  3:43668/94‑1 (MQ=255)
                                                                                                                                  aCCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGAt                                                                    >  4:43668/1‑94 (MQ=255)
                                                                                                                                     cAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTCGAGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGACGACCt           >  3:108334/1‑149 (MQ=255)
                                                                                                                                       tagcGCAAGGCTCAACCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCt           >  5:24475/4‑148 (MQ=255)
                                                                                                                                            cAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTcgc                                 <  6:10068/121‑1 (MQ=255)
                                                                                                                                            cAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTcgc                                 >  5:10068/1‑121 (MQ=255)
                                                                                                                                               gggTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAgcgc  <  7:215857/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
TCGCTGGCAGCCAACGACCAGGTGCTGATGCTGCTGGTGGCCGACCAGGGCTGCCGGCCTTTGGCCGATATCAGCGTGCAAGCCTTCGGCAAGACCGCGCAATGCATCGAGCGGGCCTTGTCGGTGTTTACCAACCGCAAGGGTTGAGCCTTGCGCTACAATCGCCCCTCTTTCCACACTGGAGACCGTGGATGACCGACTTTTCCGGCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGC  >  NC_002947/5577668‑5577956

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: