Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,994 T→G 5.0% intergenic (‑56/‑37) crcB ← / → ybeM F(‑) channel/deamidated glutathione amidase, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,9940TG 5.0% 288.8 / 3.2 93intergenic (‑56/‑37)crcB/ybeMF(‑) channel/deamidated glutathione amidase, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (49/39);  new base G (3/2);  total (52/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCC  >  NC_000913/657846‑658142
                                                                                                                                                    |                                                                                                                                                    
cTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAAtcgtc                                                                                                                                                     >  2:870419/1‑150 (MQ=255)
  gATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATcgtcgt                                                                                                                                                   <  1:65016/150‑1 (MQ=255)
   aTGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACGCTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGtt                                                                                                                                                  >  1:751552/1‑150 (MQ=255)
        gTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGa                                                                                                                                             >  1:524646/1‑150 (MQ=255)
        gTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGa                                                                                                                                             >  1:1608230/1‑150 (MQ=255)
            gTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAAttt                                                                                                                                         >  1:307748/1‑150 (MQ=255)
                aaTCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACg                                                                                                                                     <  2:383618/150‑1 (MQ=255)
                 aTCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGa                                                                                                                                    >  1:1217499/1‑150 (MQ=255)
                  tCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGAc                                                                                                                                   >  2:39146/1‑150 (MQ=255)
                  tCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGAc                                                                                                                                   >  2:1542242/1‑150 (MQ=255)
                    gCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGActc                                                                                                                                 >  2:1123009/1‑150 (MQ=255)
                     cATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGActct                                                                                                                                >  2:880785/1‑150 (MQ=255)
                      aTACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTg                                                                                                                               <  1:1073444/150‑1 (MQ=255)
                            aaCAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACagga                                                                                                                         >  2:975484/1‑150 (MQ=255)
                              cAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACaggagg                                                                                                                       >  2:1449826/1‑150 (MQ=255)
                              cAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACaggagg                                                                                                                       <  2:524646/150‑1 (MQ=255)
                               aGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGt                                                                                                                      <  1:985198/150‑1 (MQ=255)
                                 ccATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTgg                                                                                                                    >  2:1066448/1‑150 (MQ=255)
                                 ccATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTgg                                                                                                                    >  2:809581/1‑150 (MQ=255)
                                  cATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAACCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGc                                                                                                                   >  2:728050/1‑150 (MQ=255)
                                     ctcGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAAt                                                                                                                >  1:802063/1‑150 (MQ=255)
                                     ctcGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAAt                                                                                                                <  1:11763/150‑1 (MQ=255)
                                       cGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGc                                                                                                              >  1:950011/1‑150 (MQ=255)
                                          cACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTgg                                                                                                           >  2:1519869/1‑150 (MQ=255)
                                           aCGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGt                                                                                                          <  1:446368/150‑1 (MQ=255)
                                             gCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCGAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCAAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTg                                                                                                        <  1:1517104/150‑1 (MQ=255)
                                               ttCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGTATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGcc                                                                                                      <  1:1049338/150‑1 (MQ=255)
                                                 cccGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGc                                                                                                    >  2:586627/1‑150 (MQ=255)
                                                 cccGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACGCTGACAGGAGGTGGCAATGCTGGTTGCCGc                                                                                                    >  1:1256892/1‑150 (MQ=255)
                                                 cccGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAAGCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGc                                                                                                    <  2:148557/150‑1 (MQ=255)
                                                  ccGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTGTTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCa                                                                                                   <  2:1389478/150‑1 (MQ=255)
                                                   cGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAg                                                                                                  <  1:661349/150‑1 (MQ=255)
                                                   cGTACCACCGCCAATAAAAACTCCTAAAAGAACCGGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGCCACCCCGAATAATCGTCGTTGGCGAAGTTTACGACTCTGACAGGAGGTGGCAATGCTGGGTGCCGCAg                                                                                                  <  2:1558718/150‑1 (MQ=255)
                                                     tACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGa                                                                                                >  2:514876/1‑150 (MQ=255)
                                                       ccaccGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTCCAAAATCCTGCTAATTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACa                                                                                              >  2:1318801/1‑150 (MQ=255)
                                                          ccGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGtt                                                                                           >  1:1033813/1‑150 (MQ=255)
                                                           cGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGttt                                                                                          <  1:39146/150‑1 (MQ=255)
                                                             ccAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGc                                                                                        <  2:1608230/150‑1 (MQ=255)
                                                              cAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTgg                                                                                                           <  2:264259/130‑1 (MQ=255)
                                                              cAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTgg                                                                                                           >  1:264259/1‑130 (MQ=255)
                                                                 tAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGtt                                                                                    <  1:586627/150‑1 (MQ=255)
                                                                   aaaaCTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTAc                                                                                  >  2:1404204/1‑150 (MQ=255)
                                                                   aaaaCTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTAc                                                                                  <  1:654719/150‑1 (MQ=255)
                                                                      aCTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATc                                                                               <  1:1519869/150‑1 (MQ=255)
                                                                          cTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCtgtg                                                                           <  1:92419/150‑1 (MQ=255)
                                                                          cTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCtgtg                                                                           >  1:32398/1‑150 (MQ=255)
                                                                          cTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCtgtg                                                                           <  1:738373/150‑1 (MQ=255)
                                                                             aaaGAAGTTGTAACCATGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGATGGTTGCCGCAGGCCAGTTTGCTGTTACATCTGTGTgg                                                                        >  2:328769/1‑150 (MQ=255)
                                                                             aaaGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTgg                                                                        >  1:1234637/1‑150 (MQ=255)
                                                                                  agTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTTATAATTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGaaaa                                                                   <  2:868346/150‑1 (MQ=255)
                                                                                  agTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGaaaa                                                                   <  1:1542242/150‑1 (MQ=255)
                                                                                   gTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACATTTTGCTGTTACATCTGTGTGGGAAAAg                                                                  <  1:1443080/150‑1 (MQ=255)
                                                                                   cttGTAAAACCTCAAAAACCTTCCAATTTATTTGTTTTTCTGAAAACTGTTAAGGCGAAAAATCCGCGTTGGCCTAATTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAg                                                                  <  1:990066/149‑1 (MQ=255)
                                                                                     tGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGaa                                                                >  1:1458976/1‑150 (MQ=255)
                                                                                      gTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAAc                                                               <  1:514876/150‑1 (MQ=255)
                                                                                       tAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACg                                                              >  2:690394/1‑150 (MQ=255)
                                                                                        aaCACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGc                                                             <  1:1190195/150‑1 (MQ=255)
                                                                                        aaCAATGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATAGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGc                                                             >  2:1602294/1‑150 (MQ=255)
                                                                                         acacTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCt                                                            <  1:842006/150‑1 (MQ=255)
                                                                                              gCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGAt                                                       >  1:953737/1‑150 (MQ=255)
                                                                                              gCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGAt                                                       >  1:699915/1‑150 (MQ=255)
                                                                                              gCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGAt                                                       >  2:726675/1‑150 (MQ=255)
                                                                                              gCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGAt                                                       >  1:774837/1‑150 (MQ=255)
                                                                                              gCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGAt                                                       >  2:1487137/1‑150 (MQ=255)
                                                                                                 aaaTCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTtg                                                    >  2:32044/1‑150 (MQ=255)
                                                                                                 aaaTCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTtg                                                    <  1:870419/150‑1 (MQ=255)
                                                                                                   aTCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTtgtg                                                  <  1:915190/150‑1 (MQ=255)
                                                                                                   aTCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTtgtg                                                  <  1:1404204/150‑1 (MQ=255)
                                                                                                     ccTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGcc                                                >  2:390190/1‑150 (MQ=255)
                                                                                                      cTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCtgtgtg                                                                         >  1:377364/1‑124 (MQ=255)
                                                                                                      cTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCtgtgtg                                                                         <  2:377364/124‑1 (MQ=255)
                                                                                                      cTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCt                                               >  2:1342868/1‑150 (MQ=255)
                                                                                                       tGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTCTGTGGGAAAAGAACGCTGAGATTTGTGCCTc                                              >  2:668479/1‑150 (MQ=255)
                                                                                                         cTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGt                                            <  1:690394/150‑1 (MQ=255)
                                                                                                          tATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGtt                                           <  1:728050/150‑1 (MQ=255)
                                                                                                              tgatttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATg                                       <  2:1217499/150‑1 (MQ=255)
                                                                                                               gatttgTATGAGTGATAAGTGTAACGCCGAATAATCGGCGTTGGCGAATTTTACGACGCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATgg                                      <  2:1256892/150‑1 (MQ=255)
                                                                                                                 tttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGTCTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGgcg                                    >  2:260646/1‑150 (MQ=255)
                                                                                                                  ttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCGTCGTTGATGgcgc                                   <  2:642981/150‑1 (MQ=255)
                                                                                                                  ttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGgcgc                                   <  1:975484/150‑1 (MQ=255)
                                                                                                                   tgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCa                                  <  2:445654/150‑1 (MQ=255)
                                                                                                                    gTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAg                                 >  2:56753/1‑150 (MQ=255)
                                                                                                                    gTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAg                                 <  2:1394731/150‑1 (MQ=255)
                                                                                                                      aTGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAggc                               >  2:270642/1‑150 (MQ=255)
                                                                                                                       tGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAggcg                              >  1:1370498/1‑150 (MQ=255)
                                                                                                                           tGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGcggcgg                          <  2:1522381/150‑1 (MQ=255)
                                                                                                                            gATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGa                         >  1:1547857/1‑150 (MQ=255)
                                                                                                                               aaGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGaaaa                      <  1:582767/150‑1 (MQ=255)
                                                                                                                                 gtgtAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAacg                    >  1:799834/1‑150 (MQ=255)
                                                                                                                                 cgccaaCGACGATTATTCGGCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAacg                    >  2:214276/5‑150 (MQ=255)
                                                                                                                                    caaCGACGATTATTCGGCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAacgacg                 >  1:1321076/2‑150 (MQ=255)
                                                                                                                                    caaCGACGATTATTCGGCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAacgacg                 >  2:1011113/2‑150 (MQ=255)
                                                                                                                                      aCGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCa               >  2:1223132/1‑150 (MQ=255)
                                                                                                                                        gCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAGACGACGCATc             >  1:521642/1‑150 (MQ=255)
                                                                                                                                                atCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTg     >  2:126100/1‑150 (MQ=255)
                                                                                                                                                 tcgtcgTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGc    <  1:1318801/150‑1 (MQ=255)
                                                                                                                                                 tcgtcgTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGc    <  2:950011/150‑1 (MQ=255)
                                                                                                                                                   gtcgtTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGccc  <  1:1066448/150‑1 (MQ=255)
                                                                                                                                                   gtcgtTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGccc  <  1:809581/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCC  >  NC_000913/657846‑658142

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