Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 87,493 T→A 5.9% N98I (AAT→ATT)  PP_0076 ← choline betaine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294787,4930TA 5.9% 345.2 / 8.6 113N98I (AAT→ATT) PP_0076choline betaine‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (48/58);  new base A (5/2);  total (53/60)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 114 total aligned reads displayed.

AGCCGATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAACAG  >  NC_002947/87344‑87641
                                                                                                                                                      |                                                                                                                                                    
aGCCGATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAt                                                                                                                                                      <  2:1742226/150‑1 (MQ=255)
  ccGATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTg                                                                                                                                                    <  1:1546030/150‑1 (MQ=255)
  ccGATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGGCGGGCACTGCGAGGGTGAACTCGGTTTCCTCTATGTTGCGCGACAGCTGCTGCACATCGCCATTg                                                                                                                                                    <  1:1045095/150‑1 (MQ=255)
    gATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTATTGTGTCAATTCAGCAAAATAATTCAACCCGGCATTGAATAATTAAGTAAAGCATTCGAAGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCACCATTGGc                                                                                                                                                  >  1:1751793/1‑150 (MQ=255)
      ttCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGa                                                                                                                                                <  2:946387/150‑1 (MQ=255)
      ttCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGa                                                                                                                                                <  2:2304996/150‑1 (MQ=255)
       tCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGAt                                                                                                                                               <  1:2363570/150‑1 (MQ=255)
       tCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGGTTCGCGACAGCTGCTGCACATCGCCATTGGCGAt                                                                                                                                               <  1:789487/150‑1 (MQ=255)
                ttttttGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTc                                                                                                                                      <  1:2290373/150‑1 (MQ=255)
                  ttttGTCGAACTGCCCGGCGTGTTTTTGCAGGTCATCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCAt                                                                                                                                    >  1:2435713/1‑150 (MQ=255)
                          aaCTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCt                                                                                                                            >  1:1419512/1‑150 (MQ=255)
                          aaCTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCt                                                                                                                            >  2:1231109/1‑150 (MQ=255)
                          aaCTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGAAAGCCt                                                                                                                            >  2:539933/1‑150 (MQ=255)
                          aaCTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCAGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCt                                                                                                                            >  1:2102542/1‑150 (MQ=255)
                                ccGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGa                                                                                                                                                <  1:1889060/124‑1 (MQ=255)
                                   gCGTGTTTTTGCAGGTCAGCGAAGTTCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGcc                                                                                                                   <  2:116733/150‑1 (MQ=255)
                                     gtgtTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCgg                                                                                                                 >  2:88559/1‑150 (MQ=255)
                                         ttttGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGGAGCCTTGCTGCGCCGGCATc                                                                                                             <  2:1551298/150‑1 (MQ=255)
                                                ggTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTc                                                                                                      >  1:196474/1‑150 (MQ=255)
                                                ggTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTc                                                                                                      <  2:2112560/150‑1 (MQ=255)
                                                ggTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAATGGCGAt                                                                                                                                               <  2:668411/109‑1 (MQ=255)
                                                ggTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAATGGCGAt                                                                                                                                               >  1:668411/1‑109 (MQ=255)
                                                  tCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGcc                                                                                                                   >  1:1229963/1‑135 (MQ=255)
                                                  tCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGcc                                                                                                                   <  2:1229963/135‑1 (MQ=255)
                                                  tCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGGTGCGCGACAGCGGCGGCACGTCGCCATGGGCGATGAACTTGTCATGGAAGCCTAACAGCGGCGGCATCCAGTTTCCg                                                                                                    >  1:1306345/1‑150 (MQ=255)
                                                   cAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGa                                                                                                   <  2:1736339/150‑1 (MQ=255)
                                                   cAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGa                                                                                                   <  1:252300/150‑1 (MQ=255)
                                                   aaGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGa                                                                                                   <  2:1705040/149‑1 (MQ=255)
                                                        aaGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAg                                                                                              <  1:1295772/150‑1 (MQ=255)
                                                        aaGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAg                                                                                              <  2:671663/150‑1 (MQ=255)
                                                         aGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAgg                                                                                             <  2:1751793/150‑1 (MQ=255)
                                                          gTCCTTCACCCCGGCATTCCACCCCAGGCCAGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGc                                                                                            <  1:1181600/150‑1 (MQ=255)
                                                          gTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGc                                                                                            >  2:1467854/1‑150 (MQ=255)
                                                            ccTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGt                                                                                          >  2:1914461/1‑150 (MQ=255)
                                                             cTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTg                                                                                                                                                    >  1:234876/1‑91 (MQ=255)
                                                             cTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTc                                                                                         <  2:1419512/150‑1 (MQ=255)
                                                             cTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTc                                                                                         >  1:2212833/1‑150 (MQ=255)
                                                                 aCCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGAcc                                                                                     >  1:246151/1‑150 (MQ=255)
                                                                 aCCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGAcc                                                                                     >  2:844242/1‑150 (MQ=255)
                                                                  ccccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCt                                                                                    <  2:1680683/150‑1 (MQ=255)
                                                                   cccNGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTg                                                                                   >  2:1510633/1‑150 (MQ=255)
                                                                   cccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTg                                                                                   <  2:196474/150‑1 (MQ=255)
                                                                   cccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTg                                                                                   >  2:905189/1‑150 (MQ=255)
                                                                   cccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTg                                                                                   >  1:1531936/1‑150 (MQ=255)
                                                                    ccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGCACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTGCCGAGGAAGGCGGCGCCCTgg                                                                                  >  2:516037/1‑150 (MQ=255)
                                                                         aTTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTc                                                                             >  1:699204/1‑150 (MQ=255)
                                                                          ttCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGAt                                                                                                                                               >  1:1136582/1‑83 (MQ=255)
                                                                          ttCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGAt                                                                                                                                               <  2:1136582/83‑1 (MQ=255)
                                                                          ttCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTcc                                                                            >  1:569694/1‑150 (MQ=255)
                                                                           tCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCt                                                                           <  1:1510633/150‑1 (MQ=255)
                                                                              acacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTg                                                                                                                                                    >  2:1105557/1‑74 (MQ=255)
                                                                              acacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTg                                                                                                                                                    <  1:1105557/74‑1 (MQ=255)
                                                                              acacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGa                                                                        <  1:1914461/150‑1 (MQ=255)
                                                                              acacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGa                                                                        >  1:560080/1‑150 (MQ=255)
                                                                               cacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAg                                                                       <  2:875679/150‑1 (MQ=255)
                                                                                acaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAgg                                                                      <  2:329956/150‑1 (MQ=255)
                                                                                  aTAGTCGGGCAATGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGgcc                                                                    >  1:2434373/1‑150 (MQ=255)
                                                                                         gggCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTa                                                              <  1:1231109/150‑1 (MQ=255)
                                                                                           gCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGTGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAga                                                            <  2:1212929/150‑1 (MQ=255)
                                                                                           gCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAga                                                            <  1:1467854/150‑1 (MQ=255)
                                                                                               tGCGAGGGTGAAGTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAgatgat                                                        <  1:5033/150‑1 (MQ=255)
                                                                                               tGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAgatgat                                                        <  2:1531936/150‑1 (MQ=255)
                                                                                                     ggTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGa                                                                                                   >  2:158253/1‑101 (MQ=255)
                                                                                                     ggTGAACTCGGTTCCCTCCAGGTTGCCCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGa                                                                                                   <  1:158253/101‑1 (MQ=255)
                                                                                                       tGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCa                                                                                                                                     >  1:833254/1‑65 (MQ=255)
                                                                                                       tGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCa                                                                                                                                     <  2:833254/65‑1 (MQ=255)
                                                                                                           cTCGGTTCCGTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGAGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAg                                            <  1:816285/150‑1 (MQ=255)
                                                                                                                ttCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGt                                       <  1:905189/150‑1 (MQ=255)
                                                                                                                   ccTCCAGGTTGCGCGCCAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGa                                    <  1:783526/150‑1 (MQ=255)
                                                                                                                    cgccAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCa                                                                                                           <  1:643879/76‑1 (MQ=255)
                                                                                                                    cTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCa                                                                                                           >  2:643879/1‑78 (MQ=255)
                                                                                                                    cTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGAt                                   >  1:2406768/1‑150 (MQ=255)
                                                                                                                     tCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATc                                  >  2:1741621/1‑150 (MQ=255)
                                                                                                                       cAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTcatatcg                                >  1:1646800/1‑144 (MQ=255)
                                                                                                                          gcggcgcTACAGCTGCTGCACAGCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCAc                             <  2:2350842/147‑1 (MQ=255)
                                                                                                                          gTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCAc                             >  1:1673058/1‑150 (MQ=255)
                                                                                                                               gcgACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGCCCTTGAGGCCGCCGTAgatgat                                                        <  2:1557718/118‑1 (MQ=255)
                                                                                                                               gcgACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGCCCTTGAGGCCGCCGTAgatgat                                                        >  1:1557718/1‑118 (MQ=255)
                                                                                                                                  aCAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTg                                                                                   >  1:2098885/1‑88 (MQ=255)
                                                                                                                                  aCAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTg                                                                                   <  2:2098885/88‑1 (MQ=255)
                                                                                                                                    aGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAAccc                   <  2:2406768/150‑1 (MQ=255)
                                                                                                                                     gctgctGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCa                  >  1:763010/1‑150 (MQ=255)
                                                                                                                                       tgctgcATATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCAGTTGTATATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAgg                >  1:743887/1‑150 (MQ=255)
                                                                                                                                           gttcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTcc                                                                            >  2:1322316/4‑86 (MQ=255)
                                                                                                                                           gttcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCttt            >  2:636125/4‑150 (MQ=255)
                                                                                                                                            ttcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTc           >  2:2348626/3‑150 (MQ=255)
                                                                                                                                             acaTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTcc          <  2:246151/150‑1 (MQ=255)
                                                                                                                                              caTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCa         <  1:353935/150‑1 (MQ=255)
                                                                                                                                              caTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCa         <  1:88559/150‑1 (MQ=255)
                                                                                                                                               aTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCCTCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAg        >  2:806350/1‑150 (MQ=255)
                                                                                                                                               aTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCtt                                                                          >  2:2426285/1‑84 (MQ=255)
                                                                                                                                               aTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCa                              >  2:353687/1‑128 (MQ=255)
                                                                                                                                               aTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCa                              <  1:353687/128‑1 (MQ=255)
                                                                                                                                               aTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAg        >  2:352560/1‑150 (MQ=255)
                                                                                                                                               aTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACGCAGGCTTTCCAg        >  1:1698275/1‑150 (MQ=255)
                                                                                                                                                tCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGt       >  2:376197/1‑150 (MQ=255)
                                                                                                                                                  gCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTaa     <  2:1861078/150‑1 (MQ=255)
                                                                                                                                                   ccATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAAc    >  2:2125314/1‑150 (MQ=255)
                                                                                                                                                    cATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGACGTCAATCTTCACCTGGTAACCCAGGCTTTCCAGTAACa   >  2:1145692/1‑150 (MQ=255)
                                                                                                                                                     aTTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAACAg  >  1:903962/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AGCCGATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAACAG  >  NC_002947/87344‑87641

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: