Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,584,198 G→A 5.5% P477L (CCG→CTG)  glgE ← alpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,584,1980GA 5.5% 210.4 / 2.3 73P477L (CCG→CTG) glgEalpha‑1,4‑glucan:maltose‑1‑phosphate maltosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (33/35);  new base A (2/2);  total (35/38)
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.89e-01

TACTCCTCCTTGCCCGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATC  >  NC_002947/4584051‑4584347
                                                                                                                                                    |                                                                                                                                                       
tgctcctcCTTGCCCGGCAACGGCTGGCCCGCGCACAGTTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGTTCGCCAGCGCGGCACGGATCAGTAAACCTTCGCGGCCTTAGTTGTGCAGGAAAAACTGGTTTATGTCCTgt                                                                                                                                                       <  2:2232322/148‑1 (MQ=255)
tACTCCTCCTTGCCCGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGgt                                                                                                                                                       <  2:677481/150‑1 (MQ=255)
tACTCCTCCTTGCCCGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGgt                                                                                                                                                       <  1:1217079/150‑1 (MQ=255)
 aCTCCTCCTTACCCGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAAGACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCAAGGAACAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGgtg                                                                                                                                                      >  1:2176174/1‑150 (MQ=255)
   tcctccTTGCCCGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACCTGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTa                                                                                                                                                    >  2:51279/1‑150 (MQ=255)
   tcctacTTGCCAGGAACAGGCGTGACCTCGCAAAGCTAGAAGCCAGAATACATGCCCCACAAGCCTGAGCCAATGGTCGCAAGCGCGGCACGGATAAGGAAAAATACGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTa                                                                                                                                                    >  2:2231072/1‑150 (MQ=255)
     ctcctTGCCCGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTAtt                                                                                                                                                  >  1:2035346/1‑150 (MQ=255)
     ctcctTGCCCGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTAtt                                                                                                                                                  >  2:622486/1‑150 (MQ=255)
         ttGCCCGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACg                                                                                                                                              >  2:2147898/1‑150 (MQ=255)
           gCCCGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCAGTGTATTGACgaa                                                                                                                                            >  1:883644/1‑150 (MQ=255)
             ccGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACgaaga                                                                                                                                          >  1:1025332/1‑150 (MQ=255)
             ccGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCACACAAGGCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGGGTGCAGGAAAAACGGGTTGATGGCCGGTGTCTTGAcgacga                                                                                                                                          >  2:13530/1‑150 (MQ=255)
               ggCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCACCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGaagaag                                                                                                                                        <  2:51829/150‑1 (MQ=255)
                gCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGt                                                                                                                                       >  1:1142831/1‑150 (MQ=255)
                   aCGGCGTGCCCTCGCACAGCTCGAAGCCCGCATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTgg                                                                                                                                    >  1:1022160/1‑150 (MQ=255)
                      gCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGcc                                                                                                                                 >  1:710312/1‑150 (MQ=255)
                            ccTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGc                                                                                                                           >  1:78588/1‑150 (MQ=255)
                               cGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCACATGTGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACt                                                                                                                        <  1:1778882/150‑1 (MQ=255)
                                    aGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTc                                                                                                                   <  2:2035346/150‑1 (MQ=255)
                                         gAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCAcgg                                                                                                              >  1:533398/1‑150 (MQ=255)
                                           aGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGACGGTAGCACTGGCTCCAcggcg                                                                                                            <  1:622486/150‑1 (MQ=255)
                                            gCCCGAATACATGCCCCACAAGAATGAGCGCATGGGCGCCATCGCGGCACGGATCAGTAAACCTGCGCGGCCTGAGGCGCCAGCCAGTAACCGGCGGATGTCCGGTGTATTGACGAAGAAGTTGGGCAGGTCGCACTGTCTCCAcggcgg                                                                                                           >  1:1695480/1‑150 (MQ=255)
                                              ccGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTgg                                                                                                                                    >  1:748583/1‑123 (MQ=255)
                                              ccGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTgg                                                                                                                                    <  2:748583/123‑1 (MQ=255)
                                              ccGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCAcggcgg                                                                                                           <  2:1613339/148‑1 (MQ=255)
                                              ccGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCAcggcgg                                                                                                           >  1:1613339/1‑148 (MQ=255)
                                                gAATACATGCCCCAAAAGCCTGAGCCTATGGTCGCCAGCGCGGCACGGATCAGCAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGa                                                                                                       <  1:2231072/150‑1 (MQ=255)
                                                  aTACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGAtt                                                                                                     <  2:648550/150‑1 (MQ=255)
                                                     cctGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAg                                                                                                  <  2:626917/148‑1 (MQ=255)
                                                     cATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAg                                                                                                  <  2:395444/150‑1 (MQ=255)
                                                     cATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAg                                                                                                  <  2:1912400/150‑1 (MQ=255)
                                                        gCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAgctg                                                                                               >  1:1456466/1‑150 (MQ=255)
                                                               aaGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGaaa                                                                                       <  2:1025332/150‑1 (MQ=255)
                                                                    tGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTc                                                                                 <  2:2197367/150‑1 (MQ=255)
                                                                      aGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgc                                                                               <  1:717902/150‑1 (MQ=255)
                                                                      aGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgc                                                                               <  2:2170963/150‑1 (MQ=255)
                                                                       gCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgcg                                                                              <  2:78588/150‑1 (MQ=255)
                                                                       gCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCAGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgag                                                                              <  2:883644/150‑3 (MQ=255)
                                                                        cccATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCgcgc                                                                             >  1:170441/1‑150 (MQ=255)
                                                                         ccATGGTCGCCAGCGCGGCACCGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCGGATTGAGCTGCTCGAAATATTCGCGCa                                                                            <  1:1473406/150‑1 (MQ=255)
                                                                             ggTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTc                                                                        <  2:1142831/150‑1 (MQ=255)
                                                                              gTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTcc                                                                       <  2:1695480/150‑1 (MQ=255)
                                                                              gTCGCCAGCGCGGCACGGATCAGGAAACCTCAGCGGCCTGAGGAGTGCAGGAAAAACGGGTTGATGTCGGTTGGATTGAAGAAGAAGCTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAAAATTCGCGCAGCTcc                                                                       >  1:1332899/1‑150 (MQ=255)
                                                                                cGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCGCGAATATTTCGAGCAGCTcaat                                                                     <  2:238089/150‑4 (MQ=255)
                                                                                     gcgcGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGGAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTtgt                                                                <  1:51279/150‑1 (MQ=255)
                                                                                        cgGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGt                                                             <  2:170441/150‑1 (MQ=255)
                                                                                          gCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCAc                                                                                                                >  1:1928489/1‑99 (MQ=255)
                                                                                          gCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCAc                                                                                                                <  2:1928489/99‑1 (MQ=255)
                                                                                               gATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTcc                                                                       >  1:1890837/1‑133 (MQ=255)
                                                                                               gATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTcc                                                                       <  2:1890837/133‑1 (MQ=255)
                                                                                                  cAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGt                                                   <  1:157275/150‑1 (MQ=255)
                                                                                                    ggAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCCCGCCGCTCCTGCGTGTGGTTGCGCCAGGTGa                                                 >  2:1774750/1‑150 (MQ=255)
                                                                                                     gAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGaa                                                >  2:1522409/1‑150 (MQ=255)
                                                                                                         ccTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAg                                            <  2:1456466/150‑1 (MQ=255)
                                                                                                                 gCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTAACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGct                                    >  1:1142769/1‑150 (MQ=255)
                                                                                                                           tgCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTACTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGt                                          <  1:1595130/134‑1 (MQ=255)
                                                                                                                           tgCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTACTCGAAATATTCGCGCAGCTCCTGCGTGTGGTTGCGCCAGGTGAAATAGGt                                          >  2:1595130/1‑134 (MQ=255)
                                                                                                                             cAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGcc                        <  2:1332899/150‑1 (MQ=255)
                                                                                                                              aGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAg                                            <  1:1372075/129‑1 (MQ=255)
                                                                                                                              aGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAg                                            >  2:1372075/1‑129 (MQ=255)
                                                                                                                                 gAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgcg                                                                              <  1:1169327/93‑1 (MQ=255)
                                                                                                                                 gAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTcgcg                                                                              >  2:1169327/1‑93 (MQ=255)
                                                                                                                                 gAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGAc                     >  1:156241/1‑150 (MQ=255)
                                                                                                                                   aaaaCGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCt                   >  2:2100695/1‑150 (MQ=255)
                                                                                                                                        ctgTTGATTTCCGGTGTATTGATGAAGAAGTTGGGCCTGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGccc              <  2:1048425/148‑1 (MQ=255)
                                                                                                                                         ggTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCa             <  1:1689403/150‑1 (MQ=255)
                                                                                                                                         ggTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCa             <  2:710312/150‑1 (MQ=255)
                                                                                                                                           ttGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCaag           >  1:1048425/1‑148 (MQ=255)
                                                                                                                                               tGTCCAGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCAcggcgg                                                                                                           >  1:1402186/1‑52 (MQ=255)
                                                                                                                                               tGTCCAGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCAcggcgg                                                                                                           <  2:1402186/52‑1 (MQ=255)
                                                                                                                                                 tCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGcc     >  1:753472/1‑150 (MQ=255)
                                                                                                                                                 tCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGcc     >  2:1553387/1‑150 (MQ=255)
                                                                                                                                                   cGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGATAGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCcat   >  2:1727525/1‑150 (MQ=255)
                                                                                                                                                    ggTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATCAACAGCGTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCcatc  <  1:1194996/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
TACTCCTCCTTGCCCGGCAACGGCGTGCCCTCGCACAGCTCGAAGCCCGAATACATGCCCCACAAGCCTGAGCCCATGGTCGCCAGCGCGGCACGGATCAGGAAACCTGCGCGGCCTGAGGTGTGCAGGAAAAACGGGTTGATGTCCGGTGTATTGACGAAGAAGTTGGGCCGGTAGCACTGGCTCCACGGCGGCTGATTGAGCTGCTCGAAATATTCGCGCAGCTCCTGCTTGTGGTTGCGCCAGGTGAAATAGGTATAGCTCTGTGCGTAGCCGACCTTGCCCAGGCGCGCCATC  >  NC_002947/4584051‑4584347

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