Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,546,229 G→A 7.3% intergenic (+30/‑330) nfuA → / → gntT iron‑sulfur cluster carrier protein NfuA/high‑affinity gluconate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,546,2290GA 7.3% 292.2 / 8.5 95intergenic (+30/‑330)nfuA/gntTiron‑sulfur cluster carrier protein NfuA/high‑affinity gluconate transporter
Reads supporting (aligned to +/- strand):  ref base G (50/38);  new base A (2/5);  total (52/43)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCATGGTCGATGTGACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACG  >  NC_000913/3546080‑3546378
                                                                                                                                                     |                                                                                                                                                     
tCCATGGTCGATGTGACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATg                                                                                                                                                       <  1:1107144/150‑1 (MQ=255)
    tGGTCGATGTGACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGAc                                                                                                                                                   >  1:452912/1‑150 (MQ=255)
      gTCGATGTGACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGc                                                                                                                                                 <  2:160254/150‑1 (MQ=255)
         gATGTGACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTaa                                                                                                                                              <  1:666647/150‑1 (MQ=255)
             tgACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCg                                                                                                                                          >  1:1638848/1‑150 (MQ=255)
             ggACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCAAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCg                                                                                                                                          <  2:1765661/149‑1 (MQ=255)
              gACGCTGAAAGAAGGGATCGAGAATCAGCAGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGAAATGACGCTTTAGCGt                                                                                                                                         >  2:302741/1‑150 (MQ=255)
              gACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTa                                                                                                                                               <  1:1828169/144‑1 (MQ=255)
              gACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTa                                                                                                                                               >  2:1828169/1‑144 (MQ=255)
                      aagaagGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCgg                                                                                                                                 <  2:1043165/150‑1 (MQ=255)
                      aagaagGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCgg                                                                                                                                 >  1:1720437/1‑150 (MQ=255)
                      aagaagGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCgg                                                                                                                                 <  1:1318998/150‑1 (MQ=255)
                       tgaagGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCggg                                                                                                                                <  1:58217/149‑1 (MQ=255)
                          agGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGAAATGACGCTTTAGCGTCGCATCGGGcaa                                                                                                                             <  1:302741/150‑3 (MQ=255)
                          agGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGAAATGACGCTTTAGCGTCGCATCGGGcaa                                                                                                                             <  1:1654145/150‑3 (MQ=255)
                            ggATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTAc                                                                                                                           >  1:1002784/1‑150 (MQ=255)
                            ggATCGAGAAGCAGCTGCTGAACGAATTCACGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCGCTCTTTTGTAGATTGCCCGAAATGACGCTTTAGCGTCGCATCGGGcaatc                                                                                                                           <  2:414673/150‑5 (MQ=38)
                               tCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGAAATGACGCTTTAGCGTcgcatcgg                                                                                                                                 <  2:1800544/141‑8 (MQ=255)
                               tCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGAAATGACGCTTTAGCGTcgcatcgg                                                                                                                                 >  1:1800544/1‑134 (MQ=255)
                               tCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGAAATGACGCTTTAGCGTCGCATCGGGcaatctac                                                                                                                        <  2:527283/150‑8 (MQ=38)
                                cGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTg                                                                                                                       >  1:1612710/1‑150 (MQ=255)
                                 gagaAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGa                                                                                                                      >  2:209470/1‑150 (MQ=255)
                                  agaAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGaa                                                                                                                     >  1:214339/1‑150 (MQ=255)
                                  agaAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGaa                                                                                                                     <  1:1412557/150‑1 (MQ=255)
                                    aaGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCACGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACa                                                                                                                   <  1:125516/150‑1 (MQ=255)
                                     agcagcTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACAt                                                                                                                  >  2:1312765/1‑150 (MQ=255)
                                     agcagcTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACAt                                                                                                                  >  2:684055/1‑150 (MQ=255)
                                       cagcTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACAtgt                                                                                                                >  2:561826/1‑150 (MQ=255)
                                           tgctgAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCa                                                                                                            >  2:1741259/1‑150 (MQ=255)
                                              tgAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGAc                                                                                                         >  1:166637/1‑150 (MQ=255)
                                               gAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAAc                                                                                                                    <  1:1147910/138‑1 (MQ=255)
                                               gAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAAc                                                                                                                    >  2:1147910/1‑138 (MQ=255)
                                               gAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACg                                                                                                        >  1:294393/1‑150 (MQ=255)
                                               gAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACg                                                                                                        >  1:53252/1‑150 (MQ=255)
                                                    aaTTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGgc                                                                                                   <  1:1224165/150‑1 (MQ=255)
                                                     aTTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGgcg                                                                                                  >  2:470161/1‑150 (MQ=255)
                                                      ttCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGgcgc                                                                                                 <  2:697084/150‑1 (MQ=255)
                                                           ggAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCATACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCAt                                                                                            <  1:1313287/150‑1 (MQ=255)
                                                              gCTGAAAGGTGTGTGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGc                                                                                         <  2:1481739/150‑1 (MQ=255)
                                                              gCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTgg                                                                                         >  2:1936621/1‑149 (MQ=255)
                                                              gCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGc                                                                                         >  1:1945885/1‑150 (MQ=255)
                                                                tGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGc                                                                                                                                                      >  1:388868/1‑87 (MQ=255)
                                                                tGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGc                                                                                                                                                      <  2:388868/87‑1 (MQ=255)
                                                                tGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCtt                                                                                       <  2:1694964/150‑1 (MQ=255)
                                                                 gAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCttt                                                                                      >  2:751853/1‑150 (MQ=255)
                                                                   aaGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATATCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGt                                                                                    >  2:1604788/1‑150 (MQ=255)
                                                                          gcgcGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATcc                                                                             >  2:1337866/1‑150 (MQ=255)
                                                                            gcgATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTc                                                                           >  1:13251/1‑150 (MQ=255)
                                                                              gATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCgg                                                                         <  1:684055/150‑1 (MQ=255)
                                                                              gATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCgg                                                                         <  2:915852/150‑1 (MQ=255)
                                                                                tctcACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGtt                                                                       <  1:1741259/150‑1 (MQ=255)
                                                                                tctcACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGtt                                                                       <  2:1785191/150‑1 (MQ=255)
                                                                                 ctcACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTc                                                                      >  2:1332161/1‑150 (MQ=255)
                                                                                  tcACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCg                                                                     >  2:1351240/1‑150 (MQ=255)
                                                                                   cACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGt                                                                    >  2:678490/1‑150 (MQ=255)
                                                                                     ccGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGttt                                                                  >  2:74861/1‑150 (MQ=255)
                                                                                        aaCACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTccc                                                               >  1:603317/1‑150 (MQ=255)
                                                                                         acacCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTcccc                                                              <  1:1312765/150‑1 (MQ=255)
                                                                                          cacCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCcg                                                             >  2:470963/1‑150 (MQ=255)
                                                                                            ccAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCcgcg                                                           >  1:1206161/1‑150 (MQ=255)
                                                                                            ccAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCcgcg                                                           >  1:1818227/1‑150 (MQ=255)
                                                                                                 gcggcgAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGataat                                                      <  1:1019813/150‑1 (MQ=255)
                                                                                                 gcggcgAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGataat                                                      >  2:344894/1‑150 (MQ=255)
                                                                                                    gcgAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATg                                                   >  2:1889536/1‑150 (MQ=255)
                                                                                                    gcgAACACTCCTACTACTAAGTAAACCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATg                                                   <  2:1930998/150‑1 (MQ=255)
                                                                                                      gAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGAc                                                 <  1:470161/150‑1 (MQ=255)
                                                                                                         cacTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCaa                                              >  2:1633479/1‑150 (MQ=255)
                                                                                                          acTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAAc                                             <  1:561826/150‑1 (MQ=255)
                                                                                                                 ctactaAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCAt                                      >  1:351993/1‑150 (MQ=255)
                                                                                                                   actaAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaa                                    >  1:1559314/1‑150 (MQ=255)
                                                                                                                    ctaAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaat                                   >  1:1741199/1‑150 (MQ=255)
                                                                                                                    ctaAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGAAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaat                                   >  2:1653887/1‑150 (MQ=255)
                                                                                                                    ctaAGTTATCCCCTCTTTTGTAGATTGCCAGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaat                                   >  2:1660927/1‑150 (MQ=255)
                                                                                                                    ctaAGTTATCCCCTCTTTTGTAGATTGACCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaat                                   >  1:1795663/1‑150 (MQ=255)
                                                                                                                        gTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaatttaa                               <  1:1332161/150‑1 (MQ=255)
                                                                                                                         ttATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaatttaaa                              >  2:340302/1‑150 (MQ=255)
                                                                                                                          tATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaatttaaat                             <  2:166637/150‑1 (MQ=255)
                                                                                                                          tATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaatttaaat                             <  1:678490/150‑1 (MQ=255)
                                                                                                                           aTCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCACCCTCTCATaatttaaatt                            >  2:1576534/1‑150 (MQ=255)
                                                                                                                            tCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaatttaaattt                           >  2:1245074/1‑150 (MQ=255)
                                                                                                                             ccccTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaatttaaattta                          >  1:411344/1‑150 (MQ=255)
                                                                                                                             ccccTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTACCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATaatttaaattta                          <  1:209470/150‑1 (MQ=255)
                                                                                                                                 tctTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTAcccc                      <  2:13251/150‑1 (MQ=255)
                                                                                                                                   ttttGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGc                    >  2:817602/1‑150 (MQ=255)
                                                                                                                                        tAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTgg               <  1:1633479/150‑1 (MQ=255)
                                                                                                                                        tAGATTGCCCGATGCGACGCTAAAGCGTCATATCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTgg               <  1:1604788/150‑1 (MQ=255)
                                                                                                                                         aGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGt              >  2:1326391/1‑150 (MQ=255)
                                                                                                                                          gATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTtgtg             <  2:53252/150‑1 (MQ=255)
                                                                                                                                           aTTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGa            <  1:313996/150‑1 (MQ=255)
                                                                                                                                             tGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTggcgttt          >  1:617827/1‑145 (MQ=255)
                                                                                                                                             tGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGAtt          <  1:751853/150‑1 (MQ=255)
                                                                                                                                               cccGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATtct        <  1:1245074/150‑1 (MQ=255)
                                                                                                                                               cccGATGCGACACTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTAccc                       <  2:1128586/135‑1 (MQ=255)
                                                                                                                                               cccGATGCGACACTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTAccc                       >  1:1128586/1‑135 (MQ=255)
                                                                                                                                                 cGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCa      <  1:1337866/150‑1 (MQ=255)
                                                                                                                                                 cGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCa      >  2:130691/1‑150 (MQ=255)
                                                                                                                                                     gCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACg  <  2:1638848/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCATGGTCGATGTGACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACG  >  NC_000913/3546080‑3546378

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: