Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,252,246 G→C 5.1% intergenic (‑170/+70) PP_3726 ← / ← rocE enoyl‑CoA hydratase/isomerase family protein/amino acid permease RocE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,252,2460GC 5.1% 255.5 / 3.1 79intergenic (‑170/+70)PP_3726/rocEenoyl‑CoA hydratase/isomerase family protein/amino acid permease RocE
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base C (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ATGGTGTTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGGCAA  >  NC_002947/4252099‑4252395
                                                                                                                                                    |                                                                                                                                                     
aTGGTGTTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGGAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTGTATGGCTTAAGCGATCGGCCACCGCGTa                                                                                                                                                     >  2:94562/1‑150 (MQ=255)
aTGGTGTTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTa                                                                                                                                                     >  2:1186041/1‑150 (MQ=255)
           ggaCTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACCCAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCg                                                                                                                                          <  2:1457591/150‑1 (MQ=255)
             aCTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGAc                                                                                                                                        <  2:460016/150‑1 (MQ=255)
             aCTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTACGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGAc                                                                                                                                        <  2:948691/150‑1 (MQ=255)
               ttGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGAcc                                                                                                                                       >  1:460732/1‑149 (MQ=255)
               ttGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGAcc                                                                                                                                       <  2:460732/149‑1 (MQ=255)
                             aCTAAACGAGAGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGaaa                                                                                                                        >  2:1637479/1‑150 (MQ=255)
                                   cGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGCCGGCTTTATGGCATGTGGATACACGCGCCGCTGATGCTGCGCTCGACCTCACAGGCGCTGAAAATGCAg                                                                                                                  <  1:1518713/150‑1 (MQ=255)
                                   cGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAg                                                                                                                  <  2:209358/150‑1 (MQ=255)
                                    gAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTAGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGc                                                                                                                 <  1:1592521/150‑1 (MQ=255)
                                      gtgtTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGa                                                                                                               >  1:671719/1‑150 (MQ=255)
                                       tgtTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGAc                                                                                                              >  1:1166663/1‑150 (MQ=255)
                                       tgtTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGAc                                                                                                              >  2:1423122/1‑150 (MQ=255)
                                         ttCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGgcgc                                                                                                                             >  2:96190/1‑133 (MQ=255)
                                         ttCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGgcgc                                                                                                                             <  1:96190/133‑1 (MQ=255)
                                          tCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAACGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCaa                                                                                                           <  1:1090351/150‑1 (MQ=255)
                                            aGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTc                                                                                                                                           >  2:995671/1‑116 (MQ=255)
                                            aGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTc                                                                                                                                           <  1:995671/116‑1 (MQ=255)
                                            aGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGCAAATGCAGCGACCAaca                                                                                                         >  2:700453/1‑150 (MQ=255)
                                            aGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAaca                                                                                                         >  2:375310/1‑150 (MQ=255)
                                                          ggCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGgcgc                                                                                                                             >  2:1129659/1‑116 (MQ=255)
                                                          ggCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGgcgc                                                                                                                             <  1:1129659/116‑1 (MQ=255)
                                                               ttACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTgc                                                                                      >  1:1385335/1‑150 (MQ=255)
                                                                tACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTgcg                                                                                     >  1:645859/1‑150 (MQ=255)
                                                                  cAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTgcggc                                                                                   <  1:1637479/150‑1 (MQ=255)
                                                                   aaCTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTgcggcg                                                                                  >  1:625345/1‑150 (MQ=255)
                                                                   aaCTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGAGCGCCGCCTACGCGGTGCGCGATCTCACAGgccat                                                                                                                            <  1:557467/108‑4 (MQ=21)
                                                                   aaCTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGAGCGCCGCCTACGCGGTGCGCGATCTCACAGgccat                                                                                                                            >  2:557467/1‑105 (MQ=18)
                                                                     cTTCGCACAAGTCGGTGTCTGGGGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGGGAGATCGCGCACCGCGTAGGCGGCGCTCGCACTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTc                                                                                <  1:700453/150‑1 (MQ=255)
                                                                        cGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCaaa                                                                             <  2:1166663/150‑1 (MQ=255)
                                                                          cacaACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAAcc                                                                           >  2:1597985/1‑150 (MQ=255)
                                                                           acaACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAAccc                                                                          >  1:244820/1‑150 (MQ=255)
                                                                           acaACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGAGCGCCGCCTAcgcg                                                                                                                                                 <  1:469529/79‑1 (MQ=38)
                                                                           acaACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGAGCGCCGCCTAcgcg                                                                                                                                                 >  2:469529/1‑79 (MQ=38)
                                                                            caACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGAGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAAcccc                                                                         <  1:1186041/150‑1 (MQ=255)
                                                                             aaCTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGgcgc                                                                                                                                             >  2:688177/1‑81 (MQ=255)
                                                                             aaCTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGgcgc                                                                                                                                             <  1:688177/81‑1 (MQ=255)
                                                                             aaCTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGAcc                                                                                                             >  2:522317/1‑113 (MQ=255)
                                                                             aaCTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGAcc                                                                                                             <  1:522317/113‑1 (MQ=255)
                                                                              aCTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTg                                                                       <  1:1423122/150‑1 (MQ=255)
                                                                              aCTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTg                                                                       <  1:94562/150‑1 (MQ=255)
                                                                              aCTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTg                                                                       <  2:1460477/150‑1 (MQ=255)
                                                                                tCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGgt                                                                     <  2:625345/150‑1 (MQ=255)
                                                                                        tGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCt                                                                                                      <  1:896355/109‑1 (MQ=255)
                                                                                        tGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCt                                                                                                      >  2:896355/1‑109 (MQ=255)
                                                                                          gCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGgt                                                                     >  1:1287022/1‑140 (MQ=255)
                                                                                          gCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGgt                                                                     <  2:1287022/140‑1 (MQ=255)
                                                                                               acacTAGGTCTTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAg                                                      <  2:1235510/150‑1 (MQ=255)
                                                                                                 acaAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGgtgt                                                                   <  1:249590/135‑1 (MQ=255)
                                                                                                 acaAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGgtgt                                                                   >  2:249590/1‑135 (MQ=255)
                                                                                                   aaGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGccc                                                  >  2:622477/1‑150 (MQ=255)
                                                                                                     gggCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAacac                                                                                                        >  1:416541/1‑94 (MQ=255)
                                                                                                     gggCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAacac                                                                                                        <  2:416541/94‑1 (MQ=255)
                                                                                                     gggCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGaa                                                       >  2:696512/1‑143 (MQ=255)
                                                                                                     gggCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGaa                                                       <  1:696512/143‑1 (MQ=255)
                                                                                                          ttCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACGCCTGCCGACGGCGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGc                                                    <  1:656821/141‑1 (MQ=255)
                                                                                                          ttCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGc                                                    >  2:656821/1‑141 (MQ=255)
                                                                                                                  gACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGt                                   >  1:906839/1‑150 (MQ=255)
                                                                                                                  gACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGt                                   >  1:295743/1‑150 (MQ=255)
                                                                                                                   aCGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGtt                                  <  2:295743/150‑1 (MQ=255)
                                                                                                                     ggCTTTATGGCCTGTGAGATCGCGCACCGCGTCCGCTGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGa                                >  1:1621183/1‑150 (MQ=255)
                                                                                                                     ggCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGTCCCAGGGCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGa                                <  1:522905/150‑1 (MQ=255)
                                                                                                                        tttATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGTAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACga                             >  1:995055/1‑150 (MQ=255)
                                                                                                                          tATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACgaag                           <  2:1385335/150‑1 (MQ=255)
                                                                                                                          tATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACgaag                           <  1:1639689/150‑1 (MQ=255)
                                                                                                                            tGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACgaagaa                         <  2:244820/150‑1 (MQ=255)
                                                                                                                            tGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACgaagaa                         <  2:645859/150‑1 (MQ=255)
                                                                                                                            tGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACgaagaa                         >  1:1527731/1‑150 (MQ=255)
                                                                                                                             ggCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGaagaag                        >  2:326397/1‑150 (MQ=255)
                                                                                                                               ccTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGaaaa                                                                                                                       <  2:432549/53‑1 (MQ=255)
                                                                                                                               ccTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGaaaa                                                                                                                       >  1:432549/1‑53 (MQ=255)
                                                                                                                                cTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGTAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAc                     <  2:995055/150‑1 (MQ=255)
                                                                                                                                cTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAc                     >  1:1399300/1‑150 (MQ=255)
                                                                                                                                  gtgAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAcgc                   <  2:1527731/150‑1 (MQ=255)
                                                                                                                                       aTCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCAt              >  2:1120445/1‑150 (MQ=255)
                                                                                                                                             gcACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAAccc                                                                          <  2:257267/85‑1 (MQ=255)
                                                                                                                                             gcACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAAccc                                                                          >  1:257267/1‑85 (MQ=255)
                                                                                                                                              cACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGcaca        >  1:337011/1‑150 (MQ=255)
                                                                                                                                                 cgcgTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAggg     <  2:906839/150‑1 (MQ=255)
                                                                                                                                                    gTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGGCaa  <  2:337011/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGGTGTTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGGCAA  >  NC_002947/4252099‑4252395

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