Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,030,732 A→G 8.2% intergenic (‑64/+27) 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,7320AG 8.2% 319.0 / 15.0 110intergenic (‑64/+27)tusE/yccAsulfur transfer protein TusE/modulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base A (54/47);  new base G (3/6);  total (57/53)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 112 total aligned reads displayed.

AGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGTTGCTGGT  >  NC_000913/1030585‑1030881
                                                                                                                                                   |                                                                                                                                                     
agccagcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAgg                                                                                                                                                     <  1:804645/150‑1 (MQ=255)
   cagcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAATAAAAAACAAAGCCTCGCTTTATGCGGGGCtt                                                                                                                                                  <  2:464130/150‑2 (MQ=255)
     gcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAgg                                                                                                                                                     <  1:1611753/145‑1 (MQ=255)
     gcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAgg                                                                                                                                                     >  2:1611753/1‑145 (MQ=255)
      cGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGa                                                                                                                                               <  2:883290/150‑1 (MQ=255)
       ggCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGaa                                                                                                                                              >  2:1367685/1‑150 (MQ=255)
        gCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCTCGCTTTATGCGGGGCTTTgttt                                                                                                                                             <  2:56336/150‑4 (MQ=255)
         cTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAg                                                                                                                                            >  2:835240/1‑150 (MQ=255)
          tCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAgt                                                                                                                                           <  2:180159/150‑1 (MQ=255)
           cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCTCGCTTTATGCGGGGCTTTgtttttt                                                                                                                                          <  1:110212/150‑7 (MQ=255)
           cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCTCGCTTTATGCGGGGCTTTgtttttt                                                                                                                                          <  1:6065/150‑7 (MQ=255)
           cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAgtg                                                                                                                                          <  1:1048889/150‑1 (MQ=255)
            aCTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAgtgt                                                                                                                                         >  1:342211/1‑150 (MQ=255)
            aCTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAgtgt                                                                                                                                         >  1:1518780/1‑150 (MQ=255)
            aCTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAgtgt                                                                                                                                         <  1:885992/150‑1 (MQ=255)
                 aCTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCg                                                                                                                                    <  1:673226/150‑1 (MQ=255)
                     gctgctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTg                                                                                                                                <  1:634734/150‑1 (MQ=255)
                      ctgctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTga                                                                                                                               >  2:1157821/1‑150 (MQ=255)
                      ctgctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTga                                                                                                                               >  2:1163707/1‑150 (MQ=255)
                       tgctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTgag                                                                                                                              <  1:762945/150‑1 (MQ=255)
                            tctttGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTaa                                                                                                                         >  1:1464546/1‑150 (MQ=255)
                                tGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTg                                                                                                                                >  2:754818/1‑139 (MQ=255)
                                tGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTg                                                                                                                                <  1:754818/139‑1 (MQ=255)
                                   gaTAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGc                                                                                                                  >  2:638436/1‑150 (MQ=255)
                                       gCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGc                                                                                                              >  1:996376/1‑150 (MQ=255)
                                       gCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGc                                                                                                              >  1:413733/1‑150 (MQ=255)
                                       gCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGc                                                                                                              >  1:235344/1‑150 (MQ=255)
                                         cTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGa                                                                                                            <  2:1513922/150‑1 (MQ=255)
                                              gTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGccc                                                                                                       <  2:1070291/150‑1 (MQ=255)
                                               tATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGAGAGATTAATCGCGGCTAGCGAAGCCCa                                                                                                      >  2:1562853/1‑150 (MQ=255)
                                                aTCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAg                                                                                                     >  2:1378766/1‑150 (MQ=255)
                                                 tCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGa                                                                                                    >  2:909784/1‑150 (MQ=255)
                                                        cTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCt                                                                                                                 <  2:1193373/130‑1 (MQ=255)
                                                        cTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCt                                                                                                                 >  1:1193373/1‑130 (MQ=255)
                                                         tATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCTCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTg                                                                                                                                <  2:412176/114‑1 (MQ=255)
                                                         tATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGAGGGTg                                                                                                                                >  1:412176/1‑114 (MQ=255)
                                                              ctTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTca                                                                                             >  2:25264/1‑144 (MQ=255)
                                                              ctTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGc                                                                                       >  2:880853/1‑150 (MQ=255)
                                                              ctTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGc                                                                                       >  2:983285/1‑150 (MQ=255)
                                                                ttACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTg                                                                                     <  1:983285/150‑1 (MQ=255)
                                                                   ccTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACg                                                                                  <  2:1259565/150‑1 (MQ=255)
                                                                   ccTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACg                                                                                  <  2:1511849/150‑1 (MQ=255)
                                                                    cTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGa                                                                                 >  1:568981/1‑150 (MQ=255)
                                                                     ttCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGaa                                                                                >  1:1028247/1‑150 (MQ=255)
                                                                         aaGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAtg                                                                            >  2:32042/1‑150 (MQ=255)
                                                                          aGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAtgt                                                                           <  2:413733/150‑1 (MQ=255)
                                                                           gATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAtgtt                                                                          >  1:402375/1‑150 (MQ=255)
                                                                             tCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTgag                                                                                                                              >  2:1251108/1‑96 (MQ=255)
                                                                             tCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTgag                                                                                                                              <  1:1251108/96‑1 (MQ=255)
                                                                             tCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAtgttgt                                                                        <  1:345424/150‑1 (MQ=255)
                                                                              cAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTa                                                                       >  2:1448133/1‑150 (MQ=255)
                                                                                gCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACa                                                                     <  1:802385/150‑1 (MQ=255)
                                                                                gCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACa                                                                     >  1:1042559/1‑150 (MQ=255)
                                                                                  aTAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGc                                                                   >  2:567445/1‑150 (MQ=255)
                                                                                    aaTTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGa                                                                 <  2:627294/150‑1 (MQ=255)
                                                                                     attattACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGc                                                                   >  2:1632499/1‑147 (MQ=255)
                                                                                     attattACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGc                                                                   <  1:1632499/147‑1 (MQ=255)
                                                                                      ttattaCTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGaaa                                                               <  1:116394/150‑1 (MQ=255)
                                                                                       tattaCTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAAc                                                              <  2:1103508/150‑1 (MQ=255)
                                                                                          tactacCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATa                                                           >  1:1515361/1‑150 (MQ=255)
                                                                                            ctacCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACa                                                         <  1:1383121/150‑1 (MQ=255)
                                                                                            ctacATTAATCAGACTCGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAACTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACa                                                         >  2:1551143/1‑150 (MQ=255)
                                                                                                cTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGc                                                                   >  1:537540/1‑136 (MQ=255)
                                                                                                cTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGc                                                                   <  2:537540/136‑1 (MQ=255)
                                                                                                cTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCt                                                     >  1:1134008/1‑150 (MQ=255)
                                                                                                cTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCt                                                     >  1:395409/1‑150 (MQ=255)
                                                                                                cTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCt                                                     <  2:213467/150‑1 (MQ=255)
                                                                                                  tAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTa                                                    <  2:386788/149‑1 (MQ=255)
                                                                                                  tAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTa                                                    >  1:386788/1‑149 (MQ=255)
                                                                                                       aGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTg                                              >  1:197634/1‑150 (MQ=255)
                                                                                                            ggcggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACg                                         >  1:36437/1‑150 (MQ=255)
                                                                                                            ggcggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACg                                         >  1:1187603/1‑150 (MQ=255)
                                                                                                             gcggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGa                                        >  1:678166/1‑150 (MQ=255)
                                                                                                              cggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGaa                                       <  2:719590/150‑1 (MQ=255)
                                                                                                                 cAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAatat                                    <  2:857451/150‑1 (MQ=255)
                                                                                                                       aaCAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCg                              >  2:870904/1‑150 (MQ=255)
                                                                                                                       aaCAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCg                              <  1:1562853/150‑1 (MQ=255)
                                                                                                                       aaCAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCg                              <  1:835240/150‑1 (MQ=255)
                                                                                                                        aCAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGt                             >  1:1422590/1‑150 (MQ=255)
                                                                                                                        aCAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGt                             >  1:834761/1‑150 (MQ=255)
                                                                                                                         cAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGCGGTTGTACAGCGAAACATAGCGGCTAACCGTGGCACGAATATAGTTCGtc                            >  2:781421/1‑150 (MQ=255)
                                                                                                                         cAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGCCTAACCGTGGCACGAATATAGTTCGtc                            >  2:579973/1‑150 (MQ=255)
                                                                                                                           aaaaaCAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGtctc                          >  2:650594/1‑150 (MQ=255)
                                                                                                                           aaaaaCAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGGGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTACCCGTGGCACGAATATAGTTCGtctc                          >  1:1040673/1‑150 (MQ=255)
                                                                                                                              aaCAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCAcc                       <  1:1238650/150‑1 (MQ=255)
                                                                                                                               tcataGCCTCGCTTTATGCGGGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCg                      >  2:167465/5‑150 (MQ=255)
                                                                                                                                 ataGCCTCGCTTTATGCGGGGATATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGcc                    >  2:811443/3‑150 (MQ=255)
                                                                                                                                   aGCCTCGCTTTATGCGGGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAAtata                                   <  2:1399539/133‑1 (MQ=255)
                                                                                                                                   aGCCTCGCTTTATGCGGGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAAtata                                   >  1:1399539/1‑133 (MQ=255)
                                                                                                                                    gCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTGCAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCatg                 <  1:909784/150‑1 (MQ=255)
                                                                                                                                          cATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAatgat           <  1:1367685/150‑1 (MQ=255)
                                                                                                                                            tAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGt         >  1:840278/1‑150 (MQ=255)
                                                                                                                                            tAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGt         <  1:32042/150‑1 (MQ=255)
                                                                                                                                            tAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGt         <  1:1163707/150‑1 (MQ=255)
                                                                                                                                            tAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGt         <  1:1157821/150‑1 (MQ=255)
                                                                                                                                              aaGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGTTg       >  2:636088/1‑150 (MQ=255)
                                                                                                                                              aaGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGTTg       >  1:584225/1‑150 (MQ=255)
                                                                                                                                               aGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGTTGc      <  2:1134008/150‑1 (MQ=255)
                                                                                                                                                   aGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGTTGCTGGt  >  2:6886/1‑150 (MQ=255)
                                                                                                                                                   aGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGTTGCTGGt  <  1:638436/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGTTGCTGGT  >  NC_000913/1030585‑1030881

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: