Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,546,237 A→T 9.3% intergenic (+38/‑322) 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,2370AT 9.3% 287.1 / 14.3 96intergenic (+38/‑322)nfuA/gntTiron‑sulfur cluster carrier protein NfuA/high‑affinity gluconate transporter
Reads supporting (aligned to +/- strand):  ref base A (50/37);  new base T (4/5);  total (54/42)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GATGTGACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATG  >  NC_000913/3546089‑3546385
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                                                      >  2:1828169/1‑144 (MQ=255)
     gACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTa                                                                                                                                                      <  1:1828169/144‑1 (MQ=255)
             aagaagGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCgg                                                                                                                                        >  1:1720437/1‑150 (MQ=255)
             aagaagGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCgg                                                                                                                                        <  2:1043165/150‑1 (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                                                                                                                           >  2:1147910/1‑138 (MQ=255)
                                      gAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAAc                                                                                                                           <  1:1147910/138‑1 (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)
                                                       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                                                                                <  2:915852/150‑1 (MQ=255)
                                                                     gATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCgg                                                                                <  1:684055/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:1818227/1‑150 (MQ=255)
                                                                                   ccAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCcgcg                                                                  >  1:1206161/1‑150 (MQ=255)
                                                                                        gcggcgAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGataat                                                             >  2:344894/1‑150 (MQ=255)
                                                                                        gcggcgAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGataat                                                             <  1:1019813/150‑1 (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                              >  1:1128586/1‑135 (MQ=255)
                                                                                                                                      cccGATGCGACACTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTAccc                              <  2:1128586/135‑1 (MQ=255)
                                                                                                                                        cGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCa             >  2:130691/1‑150 (MQ=255)
                                                                                                                                        cGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCa             <  1:1337866/150‑1 (MQ=255)
                                                                                                                                         gaaatgACGCTTTAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCaa            >  2:1950212/6‑150 (MQ=255)
                                                                                                                                         gaaatgACGCTTTAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCAGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCaa            >  1:836251/6‑150 (MQ=255)
                                                                                                                                            gCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACg         <  2:1638848/150‑1 (MQ=255)
                                                                                                                                             cGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGc        <  2:1559314/150‑1 (MQ=255)
                                                                                                                                               aCGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCa      >  2:1690375/1‑150 (MQ=255)
                                                                                                                                                cGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACACCGCTCTGGTGATTCTCAAACGCCAg     >  1:1316044/1‑150 (MQ=255)
                                                                                                                                                  cTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGATGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGAt   >  1:1155710/1‑150 (MQ=255)
                                                                                                                                                   tAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATg  <  1:74861/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATGTGACGCTGAAAGAAGGGATCGAGAAGCAGCTGCTGAACGAATTCCCGGAGCTGAAAGGTGTGCGCGATCTCACCGAACACCAGCGCGGCGAACACTCCTACTACTAAGTTATCCCCTCTTTTGTAGATTGCCCGATGCGACGCTAAAGCGTCATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATG  >  NC_000913/3546089‑3546385

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: