Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,030,724 A→T 5.4% intergenic (‑56/+35) tusE ← / ← yccA sulfur transfer protein TusE/modulator of FtsH protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,030,7240AT 5.4% 358.2 / 5.7 109intergenic (‑56/+35)tusE/yccAsulfur transfer protein TusE/modulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base A (54/49);  new base T (3/3);  total (57/52)
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.90e-01

Only 100 of 109 total aligned reads displayed.

CACCACAGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCAT  >  NC_000913/1030579‑1030865
                                                                                                                                                 |                                                                                                                                             
caccacAGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATaaa                                                                                                                                           >  1:422690/1‑150 (MQ=255)
caccacAGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATaaa                                                                                                                                           >  1:1000271/1‑150 (MQ=255)
caccacAGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATaaa                                                                                                                                           >  2:1638012/1‑150 (MQ=255)
 accacAGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAg                                                                                                                                          >  1:1425779/1‑150 (MQ=255)
 accacAGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTAACTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAg                                                                                                                                          >  2:1061962/1‑150 (MQ=255)
  ccacAGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGc                                                                                                                                         >  1:1004279/1‑150 (MQ=255)
    acaGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGa                                                                                                                                       <  1:247753/150‑1 (MQ=255)
     cagccagcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAg                                                                                                                                      >  1:178873/1‑150 (MQ=255)
       gccagcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGc                                                                                                                                    >  2:232953/1‑150 (MQ=255)
         cagcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTa                                                                                                                                  <  2:422690/150‑1 (MQ=255)
          agcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCTCGCTTTATGCGGGGCttt                                                                                                                                 <  1:164405/150‑3 (MQ=255)
          agcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTAt                                                                                                                                 >  1:1065203/1‑150 (MQ=255)
             ggCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGaa                                                                                                                              >  2:1202696/1‑150 (MQ=255)
             ggCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGaa                                                                                                                              >  1:945603/1‑150 (MQ=255)
              gCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGc                                                                                                                                         >  2:312994/1‑138 (MQ=255)
              gCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGc                                                                                                                                         <  1:312994/138‑1 (MQ=255)
               cTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAg                                                                                                                            <  1:913551/150‑1 (MQ=255)
                 cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTTTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAgtg                                                                                                                          <  2:524627/150‑1 (MQ=255)
                 cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAgtg                                                                                                                          >  2:1301626/1‑150 (MQ=255)
                 cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAACCCTCGCTTTATGCGGGGCTTTgtttttt                                                                                                                          <  1:910024/150‑7 (MQ=255)
                  aCTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCATAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCTCGCTTTATGCGGGGCTTTgttttttg                                                                                                                         <  2:1109500/150‑8 (MQ=255)
                   cTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGtt                                                                                                                        <  1:534315/150‑1 (MQ=255)
                   cTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGtt                                                                                                                        >  1:1664856/1‑150 (MQ=255)
                   cTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGtt                                                                                                                        <  2:1663949/150‑1 (MQ=255)
                     ccACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAg                                                                                                                      <  1:232953/150‑1 (MQ=255)
                       aCTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCg                                                                                                                    <  2:46149/150‑1 (MQ=255)
                        cTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCgg                                                                                                                   >  2:1011540/1‑150 (MQ=255)
                         tGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCggg                                                                                                                  <  1:1668123/150‑1 (MQ=255)
                         tGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCggg                                                                                                                  <  2:1568902/150‑1 (MQ=255)
                             tgctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTgag                                                                                                              <  2:945603/150‑1 (MQ=255)
                              gctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGtt                                                                                                                        <  2:1609677/139‑1 (MQ=255)
                              gctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGtt                                                                                                                        >  1:1609677/1‑139 (MQ=255)
                              gctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTgaga                                                                                                             >  2:1616197/1‑150 (MQ=255)
                                tttctttGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGAtt                                                                                                           <  2:215510/150‑1 (MQ=255)
                                  tctttGAGATAGCCTTCGCTATCCGTTTATATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGAGGCAGTGTAACAAAAAACAAAGCCCCGAATAAACCGAGGCTATGAAAGTGATAGCGGGTGAGATTaa                                                                                                         >  2:530466/1‑150 (MQ=255)
                                   ctttGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGGAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAAt                                                                                                        >  1:622243/1‑150 (MQ=255)
                                         gaTAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGc                                                                                                  <  1:538566/150‑1 (MQ=255)
                                            aGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAg                                                                                               <  1:449075/150‑1 (MQ=255)
                                            aGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAg                                                                                               <  1:530466/150‑1 (MQ=255)
                                             gCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGc                                                                                              >  2:123411/1‑150 (MQ=255)
                                              ccTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAg                                                                                               <  2:217189/148‑1 (MQ=255)
                                              ccTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAg                                                                                               >  1:217189/1‑148 (MQ=255)
                                              ccTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCg                                                                                             <  1:137494/150‑1 (MQ=255)
                                               cTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGa                                                                                            <  1:290722/150‑1 (MQ=255)
                                               cTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGa                                                                                            <  2:1000271/150‑1 (MQ=255)
                                                ttCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGaa                                                                                           <  2:1375861/150‑1 (MQ=255)
                                                  cGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGc                                                                                         <  2:216840/150‑1 (MQ=255)
                                                    gTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGccc                                                                                       <  2:939313/150‑1 (MQ=255)
                                                               tATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTcag                                                                            >  1:823657/1‑150 (MQ=255)
                                                                 tctcTTTACCCTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTcagca                                                                          <  1:174391/150‑1 (MQ=255)
                                                                    ctTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGc                                                                       >  1:999554/1‑150 (MQ=255)
                                                                    ctTTACCTTCGAAGATCAGCATAATTATTACTACCCTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGCTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGc                                                                       <  2:312921/150‑1 (MQ=255)
                                                                     tttACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCt                                                                      <  2:178873/150‑1 (MQ=255)
                                                                       tACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGa                                                                    <  2:1535204/150‑1 (MQ=255)
                                                                         ccTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACg                                                                  >  1:616498/1‑150 (MQ=255)
                                                                            tCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAg                                                               <  2:1425779/150‑1 (MQ=255)
                                                                            tCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAg                                                               <  1:1204392/150‑1 (MQ=255)
                                                                            tCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAg                                                               >  2:1595693/1‑150 (MQ=255)
                                                                             cGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGa                                                              >  1:177796/1‑150 (MQ=255)
                                                                               aaGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAtg                                                            <  2:1065203/150‑1 (MQ=255)
                                                                                  aTCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAtgttg                                                         <  2:1215874/150‑1 (MQ=255)
                                                                                   tCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACGAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAtgttgt                                                        >  2:1127138/1‑150 (MQ=255)
                                                                                    cAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTa                                                       >  1:1391042/1‑150 (MQ=255)
                                                                                     aGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTAc                                                      >  1:1618359/1‑150 (MQ=255)
                                                                                            ttattaCTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGaaa                                               <  1:1431410/150‑1 (MQ=255)
                                                                                                  ctacCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACa                                         >  1:33720/1‑150 (MQ=255)
                                                                                                   tacCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTg                                                                                                                >  1:1470453/1‑78 (MQ=255)
                                                                                                   tacCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTg                                                                                                                <  2:1470453/78‑1 (MQ=255)
                                                                                                    acCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAgg                                       >  2:517455/1‑150 (MQ=255)
                                                                                                    acCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAgg                                       >  2:523106/1‑150 (MQ=255)
                                                                                                    acCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAgg                                       >  1:1051786/1‑150 (MQ=255)
                                                                                                     ccTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGc                                      <  1:78087/150‑1 (MQ=255)
                                                                                                       ttAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTa                                    >  1:1199563/1‑150 (MQ=255)
                                                                                                             aGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACAt                                            <  2:314235/136‑1 (MQ=255)
                                                                                                             aGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTg                              <  1:1638012/150‑1 (MQ=255)
                                                                                                             aGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTg                              <  2:874706/150‑1 (MQ=255)
                                                                                                                cTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCa                           <  1:1202696/150‑1 (MQ=255)
                                                                                                                  ggcggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAGCATACa                                         >  1:880112/1‑134 (MQ=255)
                                                                                                                  ggcggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACa                                         <  2:880112/134‑1 (MQ=255)
                                                                                                                  ggcggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACg                         >  1:252850/1‑150 (MQ=255)
                                                                                                                   gcggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACa                                             <  2:1333651/129‑1 (MQ=255)
                                                                                                                   gcggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACa                                             >  1:1333651/1‑129 (MQ=255)
                                                                                                                   gcggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGa                        <  1:123411/150‑1 (MQ=255)
                                                                                                                    cggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGaa                       <  2:1277495/150‑1 (MQ=255)
                                                                                                                      gcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGAGAAACATACAGGCTAACCGTGGCACGAata                     >  1:367282/1‑150 (MQ=255)
                                                                                                                       cAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTg                              <  1:894258/140‑1 (MQ=255)
                                                                                                                       cAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTg                              >  2:894258/1‑140 (MQ=255)
                                                                                                                              aCAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAtg                                                            >  1:1202338/1‑103 (MQ=255)
                                                                                                                              aCAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAtg                                                            <  2:1202338/103‑1 (MQ=255)
                                                                                                                               cAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGtc            <  2:252850/150‑1 (MQ=255)
                                                                                                                               cAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGtc            >  1:720776/1‑150 (MQ=255)
                                                                                                                               cAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCACCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGtc            >  1:1264771/1‑150 (MQ=255)
                                                                                                                                aaaaaaCAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCa                           <  2:699001/134‑1 (MQ=255)
                                                                                                                                 aaaaaCAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGtctc          >  1:790504/1‑150 (MQ=255)
                                                                                                                                 aaaaaCAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGCGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCTGCAGCCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGtctc          >  1:1065850/1‑150 (MQ=255)
                                                                                                                                  ctttcataGCCTCGCTTTATGCGGGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCa         >  2:135925/8‑150 (MQ=255)
                                                                                                                                       ataGCCTCGCTTTATGCGGGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGcc    >  1:1534501/3‑150 (MQ=255)
                                                                                                                                       ataGCCTCGCTTTATGCGGGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGcc    >  1:662823/3‑150 (MQ=255)
                                                                                                                                         aGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCat  >  1:552023/1‑150 (MQ=255)
                                                                                                                                         aGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCat  >  1:1027433/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CACCACAGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCAT  >  NC_000913/1030579‑1030865

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: