Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,313,825 C→A 5.3% intergenic (+30/‑129) hemB → / → proP 5‑aminolevulinate dehydratase/osmosensory proline/betaine/H+ permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,313,8250CA 5.3% 309.6 / 4.4 94intergenic (+30/‑129)hemB/proP5‑aminolevulinate dehydratase/osmosensory proline/betaine/H+ permease
Reads supporting (aligned to +/- strand):  ref base C (47/42);  new base A (2/3);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GCCGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAG  >  NC_002947/3313679‑3313973
                                                                                                                                                  |                                                                                                                                                      
gccGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCa                                                                                                                                                       <  5:155168/148‑1 (MQ=255)
   ggTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTg                                                                                                                                                    >  2:317006/1‑148 (MQ=255)
       ccATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCg                                                                                                                                              <  8:100815/149‑1 (MQ=255)
        cATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCgg                                                                                                                                             >  7:51446/1‑149 (MQ=255)
         aTCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCCGCGGc                                                                                                                                            >  4:323336/1‑149 (MQ=255)
         aTCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGc                                                                                                                                            <  1:197327/149‑1 (MQ=255)
         aTCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGc                                                                                                                                            >  4:2844/1‑149 (MQ=255)
          tCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCt                                                                                                                                           >  1:314135/1‑149 (MQ=255)
             acgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCtttt                                                                                                                                       <  1:317006/149‑1 (MQ=255)
             acgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCtttt                                                                                                                                       <  4:222812/149‑1 (MQ=255)
                 aGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCttt                                                                                                                                   <  8:100676/149‑1 (MQ=255)
                 aGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCttt                                                                                                                                   <  8:90651/149‑1 (MQ=255)
                  ggCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTtg                                                                                                                                  <  8:37440/149‑1 (MQ=255)
                   gCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTtgt                                                                                                                                 <  8:331472/149‑1 (MQ=255)
                   gCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTtgt                                                                                                                                 <  8:51446/149‑1 (MQ=255)
                     cGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTgc                                                                                                                               <  1:272544/149‑1 (MQ=255)
                     cGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTgc                                                                                                                               >  5:201174/1‑149 (MQ=255)
                       ttgattgCGCCAATGCTTTCGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTgcgc                                                                                                                             >  6:330004/6‑149 (MQ=255)
                               gcgAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCg                                                                                                                                              >  1:189394/1‑125 (MQ=255)
                               gcgAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCg                                                                                                                                              <  2:189394/125‑1 (MQ=255)
                                 gAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTg                                                                                                                        <  7:321708/144‑1 (MQ=255)
                                 gAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTg                                                                                                                        >  8:321708/1‑144 (MQ=255)
                                 gAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGggtgg                                                                                                                   >  7:187694/1‑149 (MQ=255)
                                       aTTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTgg                                                                                                                       >  7:176503/1‑139 (MQ=255)
                                       aTTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTgg                                                                                                                       <  8:176503/139‑1 (MQ=255)
                                        ttGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGc                                                                                                             >  3:78439/1‑148 (MQ=255)
                                          ggCGCAATCAAGCGTGCTGGTGAACACCTCATCCAGACCTACTTCCACAAGGACCACGCCCGCGAAGCTATCTGAAGCCACTCCATGGTCATGAAAGCAGCTGCCAAAGCAGCGGCTATTCTTATAGCCCCGTAGGGTGCAGCCGCaaa                                                                                                          >  8:84454/1‑149 (MQ=255)
                                                 tCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGcaca                                                                                                    >  3:204657/1‑148 (MQ=255)
                                                 tCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGcaca                                                                                                    >  3:326973/1‑148 (MQ=255)
                                                    aGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACaaaa                                                                                                 >  5:165220/1‑148 (MQ=255)
                                                        tGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACa                                                                                            <  4:78439/149‑1 (MQ=255)
                                                             gTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGc                                                                                                                                            >  6:208254/1‑97 (MQ=255)
                                                             gTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGc                                                                                                                                            <  5:208254/97‑1 (MQ=255)
                                                             gTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTaa                                                                                        <  6:217530/148‑1 (MQ=255)
                                                             gTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAg                                                                                       >  4:20024/1‑149 (MQ=255)
                                                             gTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAg                                                                                       >  4:88811/1‑149 (MQ=255)
                                                             gTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAg                                                                                       >  5:194514/1‑149 (MQ=255)
                                                              tGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAgc                                                                                      >  7:143929/1‑149 (MQ=255)
                                                                    ccTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGaa                                                                                >  5:149921/1‑149 (MQ=255)
                                                                     cTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAat                                                                               >  5:123479/1‑149 (MQ=255)
                                                                     cTGATCCTGACCTACTTCGCCATGGACCTGGCGAGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGACGTTGGGTGGTGCAGCAAACGCACAAAAAAACAGTAAGAGCTGAat                                                                               >  6:261791/1‑149 (MQ=255)
                                                                         tCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGGCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAattatt                                                                           <  5:366598/149‑1 (MQ=255)
                                                                           cTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTAttt                                                                          <  2:314135/148‑1 (MQ=255)
                                                                           cTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTa                                                                         <  4:305301/149‑1 (MQ=255)
                                                                           cTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTa                                                                         <  3:20024/149‑1 (MQ=255)
                                                                           cTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTa                                                                         <  8:143929/149‑1 (MQ=255)
                                                                           cTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAACAGTAAGCGCTGAATTATTTa                                                                         <  6:165220/148‑1 (MQ=255)
                                                                            tGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAg                                                                        >  4:266332/1‑149 (MQ=255)
                                                                             gACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGaa                                                                                <  2:288568/140‑1 (MQ=255)
                                                                             gACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGaa                                                                                >  1:288568/1‑140 (MQ=255)
                                                                              aCCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAgt                                                                       <  4:326973/148‑1 (MQ=255)
                                                                                cTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTc                                                                    >  1:361100/1‑149 (MQ=255)
                                                                                     tCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgtt                                                               >  1:249110/1‑149 (MQ=255)
                                                                                      cGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgtt                                                               <  6:201174/148‑1 (MQ=255)
                                                                                      cGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgtt                                                               <  5:348802/148‑1 (MQ=255)
                                                                                       gCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTCTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGcac                                                                                                     <  1:60299/110‑1 (MQ=255)
                                                                                       gCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGcac                                                                                                     >  2:60299/1‑109 (MQ=255)
                                                                                         cATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCGTTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGaa                                                                                >  3:353146/1‑128 (MQ=255)
                                                                                         cATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCGTTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGaa                                                                                <  4:353146/128‑1 (MQ=255)
                                                                                         cATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGc                                                                                                                                                >  1:52638/1‑66 (MQ=255)
                                                                                         cATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGc                                                                                                                                                <  2:52638/66‑1 (MQ=255)
                                                                                         cATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTa                                                            >  1:133522/1‑148 (MQ=255)
                                                                                          aTGGACCTGGCGAGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTATGTGAGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGATGAATTATTTAGTGTCCTGTTGTAta                                                          >  7:324783/1‑148 (MQ=255)
                                                                                          aTGCACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGGTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATc                                                          <  3:2844/149‑1 (MQ=255)
                                                                                              aCCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGAt                                                      >  7:71479/1‑149 (MQ=255)
                                                                                                  ggCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATATGATCCCt                                                  >  3:253306/1‑149 (MQ=255)
                                                                                                    cgcgcgAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTa                                                <  8:71479/149‑1 (MQ=255)
                                                                                                     gcgcgAAGGTATCTGAAGCCAGGCCTTGGTCATTAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTag                                               <  5:330004/149‑1 (MQ=255)
                                                                                                      cgcgAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTaga                                              <  2:290973/149‑1 (MQ=255)
                                                                                                       gcgAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTagag                                             >  2:68565/1‑149 (MQ=255)
                                                                                                         gAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAACCGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCa                                           >  4:210217/1‑149 (MQ=255)
                                                                                                         gAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCa                                           >  1:217728/1‑149 (MQ=255)
                                                                                                                tCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATc                                     >  6:275378/1‑148 (MQ=255)
                                                                                                                  tGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGGATTTACTTTGTGGGCGGTTGGGTGGTGACGCAAACGAAAAAAAAAAAATTAAGCGATTAATTATTTAGTGTCATGTTGTAGATGATCCCTTAGAGAATTAATCg                                    >  1:255764/1‑148 (MQ=255)
                                                                                                                  tGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGa                                   >  8:209332/1‑148 (MQ=255)
                                                                                                                     aGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAg                                <  5:195943/148‑1 (MQ=255)
                                                                                                                        cAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGa                             >  7:310378/1‑148 (MQ=255)
                                                                                                                         aGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATa                           >  5:40699/1‑149 (MQ=255)
                                                                                                                         aGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGCTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTCTCCTGTTGTATCTGATCCCTTAGAGCATTCAACGACAGGGAt                            >  2:119213/1‑148 (MQ=255)
                                                                                                                            ccTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACgg                        <  5:303679/149‑1 (MQ=255)
                                                                                                                                ggTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTa                     <  6:149921/148‑1 (MQ=255)
                                                                                                                                   cATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTTATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGa                  >  1:308365/1‑148 (MQ=255)
                                                                                                                                    agaaaaGCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATc                                     >  4:77890/4‑128 (MQ=255)
                                                                                                                                    agaaaaGCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATc                                     <  3:77890/125‑1 (MQ=255)
                                                                                                                                    aTGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGaaa                <  3:266332/149‑1 (MQ=255)
                                                                                                                                      gAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACgg                        >  7:70776/1‑139 (MQ=255)
                                                                                                                                      gAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACgg                        <  8:70776/139‑1 (MQ=255)
                                                                                                                                      gAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAAcc              >  6:133434/1‑149 (MQ=255)
                                                                                                                                       aaaGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAAcc              <  1:68565/148‑1 (MQ=255)
                                                                                                                                        aaGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTc            <  3:323336/149‑1 (MQ=255)
                                                                                                                                          gCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTTTCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCg           >  4:269866/1‑148 (MQ=255)
                                                                                                                                           ccGCTGCATTTGGAGCGGCTTTTGTTTGTGCTCCGTTTGGTGGTGCCGCAAACTCACAAAAAAACAGTAATCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCACTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGc          <  8:324783/148‑1 (MQ=255)
                                                                                                                                            cGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCaa        <  5:275378/149‑1 (MQ=255)
                                                                                                                                             gCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAg       <  2:308365/149‑1 (MQ=255)
                                                                                                                                                 ccATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGaaaa   <  2:361100/149‑1 (MQ=255)
                                                                                                                                                  cATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAg  <  6:194514/149‑1 (MQ=255)
                                                                                                                                                  cATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAg  >  5:167083/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GCCGGTGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAG  >  NC_002947/3313679‑3313973

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: