Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,439,030 A→G 5.9% intergenic (‑2/‑61) PP_3933 ← / → PP_3934 membrane protein/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,439,0300AG 5.9% 199.2 / 2.5 68intergenic (‑2/‑61)PP_3933/PP_3934membrane protein/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base G (2/2);  total (31/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 = 6.61e-01

TGCTGCCGAGCAAGGCGGCCGGTTGCGCGGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACATGCTGCGGGAAATTCAGGCGTTCGTCAGCGTGG  >  NC_002947/4438887‑4439172
                                                                                                                                               |                                                                                                                                              
tgctgcCGAGCAAGGCGGCCGGTTGCGCGGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCgg                                                                                                                                           >  7:115029/1‑149 (MQ=255)
 gctgcCGAGCAAGGCGGCCGGTTGCGCGGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCgg                                                                                                                                           <  6:253378/148‑1 (MQ=255)
   tgcCGAGCAAGGCGGCCGGTTGCGCGGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTc                                                                                                                                        >  5:278650/1‑149 (MQ=255)
   tgcCGAGCAAGGCGGCCGGTTGCGCGGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTc                                                                                                                                        >  2:94752/1‑149 (MQ=255)
   tgcCGAGCAAGGCGGCCGGTTGCGCGGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTc                                                                                                                                        >  6:426572/1‑149 (MQ=255)
     ccGAGCAAGGCGGCCGGTTGCGCGGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTcc                                                                                                                                       <  5:120983/148‑1 (MQ=255)
          cAAGGCGGCCGGTTGCGCGGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACAGCGGCGTCGATCAACCCAGCGGGGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAg                                                                                                                                 >  1:334278/1‑149 (MQ=255)
           aaGGCGGCCGGTTGCGCGGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGt                                                                                                                                <  6:369466/149‑1 (MQ=255)
               cgggcgggTGCGCGGTTGGCAAAAAGTTGAAACACGCCGGGCTCTGGAACATGCCGCCGCCAACCACCGCCGCGTTGAATAAACCATGGGGGAATGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGaa                                                                                                                             <  6:237006/148‑1 (MQ=255)
                ggcCGGTTGCGCGGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGa                                                                                                                           <  6:278650/149‑1 (MQ=255)
                           cgGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAgg                                                                                                                <  7:25918/149‑1 (MQ=255)
                            ggAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGc                                                                                                               <  8:236826/149‑1 (MQ=255)
                               gggCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTaa                                                                                                            <  8:228696/149‑1 (MQ=255)
                                  cAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCAAAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATc                                                                                                          >  6:420208/1‑148 (MQ=255)
                                   aaCACGTTGAACAAGGCCGGGACCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAgg                                                                                                                                    >  7:258010/1‑121 (MQ=255)
                                   aaCACGTTGAACAAGGCCGGGACCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAgg                                                                                                                                    <  8:258010/121‑1 (MQ=255)
                                    acacGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCAAAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCaaa                                                                                                       <  5:420208/149‑1 (MQ=255)
                                     cacGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAg                                                                                                      <  6:408209/149‑1 (MQ=255)
                                      acGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGt                                                                                                     <  1:94752/149‑1 (MQ=255)
                                        gTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACTCCTGCAAGGCTAATCAAAGTgg                                                                                                   <  2:334278/149‑1 (MQ=255)
                                           gAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGca                                                                                                 >  3:205090/1‑148 (MQ=255)
                                           gAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGAGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGca                                                                                                 >  5:276289/1‑148 (MQ=255)
                                             aCAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCacac                                                                                              >  6:160030/1‑149 (MQ=255)
                                                 ggCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTtgt                                                                                          <  6:294899/149‑1 (MQ=255)
                                                    cGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGtt                                                                                       <  7:93140/149‑1 (MQ=255)
                                                       gATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGaa                                                                                    <  5:426572/149‑1 (MQ=255)
                                                             gATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGaaa                                                                                   <  2:187543/144‑1 (MQ=255)
                                                               tCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGca                                                                                                 >  2:94258/1‑128 (MQ=255)
                                                               tCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGca                                                                                                 <  1:94258/128‑1 (MQ=255)
                                                                   cccGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGt                                                                                                     <  3:217012/120‑1 (MQ=255)
                                                                   cccGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGt                                                                                                     >  4:217012/1‑120 (MQ=255)
                                                                   cccGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGc                                                                        >  5:364434/1‑149 (MQ=255)
                                                                      gccgccACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGa                                                                      >  6:203519/1‑148 (MQ=255)
                                                                           cacccaccGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGGTGCGGTTc                                                                                                                                        <  2:377786/77‑1 (MQ=255)
                                                                           cacccaccGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGGTGCGGTTc                                                                                                                                        >  1:377786/1‑77 (MQ=255)
                                                                             cccaccGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGc                                                               >  8:138016/1‑148 (MQ=255)
                                                                               caccGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAAc                                                            <  8:390727/149‑1 (MQ=255)
                                                                                 ccGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACga                                                          <  6:364434/149‑1 (MQ=255)
                                                                                 aggctgcgTCGATCATCCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGGTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACga                                                          <  7:350515/144‑1 (MQ=255)
                                                                                    cggcgTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACgagga                                                       <  7:138016/149‑1 (MQ=255)
                                                                                        gTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCaa                                                             >  7:327382/1‑139 (MQ=255)
                                                                                        gTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCaa                                                             <  8:327382/139‑1 (MQ=255)
                                                                                        gTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAg                                                   <  4:205090/149‑1 (MQ=255)
                                                                                        gTCGATCAACCCAGCGGAGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAg                                                   <  6:276289/149‑1 (MQ=255)
                                                                                            aTCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCa                                               >  2:19764/1‑149 (MQ=255)
                                                                                                  ccAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGGTGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAgg                                                  >  8:310686/1‑140 (MQ=255)
                                                                                                  ccAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGGTGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAgg                                                  <  7:310686/140‑1 (MQ=255)
                                                                                                  ccAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATg                                             <  8:123086/145‑1 (MQ=255)
                                                                                                  ccAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATg                                             >  7:123086/1‑145 (MQ=255)
                                                                                                    aGCGGCGAAGGCGAACACAGAAAGCTCAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTg                                       <  2:210077/149‑1 (MQ=255)
                                                                                                         cgaaggcgaaCACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACAt                                  <  5:23265/149‑1 (MQ=255)
                                                                                                              gcgaaCACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACAtgctgc                             <  5:160030/149‑1 (MQ=255)
                                                                                                                      aGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTgg                                      >  3:81952/1‑132 (MQ=255)
                                                                                                                      aGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTgg                                      <  4:81952/132‑1 (MQ=255)
                                                                                                                        aaaGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACATGCTGCGGGAAAtt                     >  3:2277/1‑147 (MQ=255)
                                                                                                                        aaaGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGAACTGGACATGCTGCGGGAAAtt                     >  7:212051/1‑147 (MQ=255)
                                                                                                                         aaGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTTGACATGCTGCGGGAAATTCAg                  <  6:204783/149‑1 (MQ=255)
                                                                                                                         aaGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTAGACATGCTGCGGGAAATTCAg                  >  5:204783/1‑149 (MQ=255)
                                                                                                                           gCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACATGCTGCGGGAAATTCAGGc                >  4:296917/1‑149 (MQ=255)
                                                                                                                             acaGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACATGCTGCGGGAAATTCAGGCGt              >  2:23042/1‑149 (MQ=255)
                                                                                                                                 cAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACATGCTGCGGGAAATTCAGGCGTTCGt          >  5:229440/1‑149 (MQ=255)
                                                                                                                                       ccATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACATGCTGCGGGAAATTCAGGCGTTc            >  4:94826/1‑141 (MQ=255)
                                                                                                                                       ccATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACATGCTGCGGGAAATTCAGGCGTTc            <  3:94826/141‑1 (MQ=255)
                                                                                                                                         atcatGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTgg                                      >  1:311225/1‑113 (MQ=255)
                                                                                                                                         atcatGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTgg                                      <  2:311225/113‑1 (MQ=255)
                                                                                                                                         atcatGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACATGCTGCGGGaaa                       >  5:200329/1‑128 (MQ=255)
                                                                                                                                         atcatGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACATGCTGCGGGaaa                       <  6:200329/128‑1 (MQ=255)
                                                                                                                                         atcatGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACATGCTGCGGGAAATACAGGCGTTCGTCAGCGTgg  >  2:90258/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TGCTGCCGAGCAAGGCGGCCGGTTGCGCGGAGGGCAACACGTTGAACAAGGCCGGGATCTGGATCAGCCCGCCGCCACCCACCGCGGCGTCGATCAACCCAGCGGCGAAGGCGAACACAGAAAGCACAGCAATATCCATCATGATGCGGTTCCAGGCAGTGAAGACGCCTGCAAGGCTAATCAAAGTGGCACACTTGTGTTGAAATGCCATATTGCGATAACTGCAACGAGGAGCAGGCCATGGCACTGGACATGCTGCGGGAAATTCAGGCGTTCGTCAGCGTGG  >  NC_002947/4438887‑4439172

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