Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,235,688 A→T 9.4% intergenic (+31/+21) fadR → / ← ycgB fatty acid metabolism regulon transcriptional regulator/SpoVR family stationary phase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,235,6880AT 9.4% 201.1 / 8.7 83intergenic (+31/+21)fadR/ycgBfatty acid metabolism regulon transcriptional regulator/SpoVR family stationary phase protein
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base T (3/5);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

GTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCC  >  NC_000913/1235539‑1235834
                                                                                                                                                     |                                                                                                                                                  
gTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCaaa                                                                                                                                                    >  2:238740/1‑150 (MQ=255)
  gAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACg                                                                                                                                                  <  1:870798/150‑1 (MQ=255)
      gcgcgcACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCTTGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGACATTCAGGGGCTATTATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAggg                                                                                                                                              <  2:504447/150‑1 (MQ=255)
         cgcACGATCAGGTGTACGAAACAGTGGGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCTGAAAATTCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAATAGCAATCCCCGCAAACGAGGGGtt                                                                                                                                           >  2:299467/1‑150 (MQ=255)
           cACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGGtttg                                                                                                                                         <  1:1310758/150‑2 (MQ=255)
            aCGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATACAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGttttt                                                                                                                                        <  2:23149/150‑1 (MQ=255)
             cGATCAGGTGTACGAATCAGCGCGTCGCTATGGGGATTAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGGtttgct                                                                                                                                       <  1:1357616/150‑4 (MQ=255)
             cGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGGtttgct                                                                                                                                       <  2:1150235/150‑4 (MQ=255)
             cGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGtttttt                                                                                                                                       <  2:960114/150‑1 (MQ=255)
             cGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGtttttt                                                                                                                                       >  1:238185/1‑150 (MQ=255)
              gATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttg                                                                                                                                      <  2:238185/150‑1 (MQ=255)
              gATCAGGGGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGGtttgctc                                                                                                                                      <  2:720367/150‑5 (MQ=255)
               aTCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttgt                                                                                                                                     >  2:1545221/1‑150 (MQ=255)
                     tgtACGAAACAGTGCGTCGCTATGGGCATGAGAGAGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGtttt                                                                                                                               >  1:1389921/1‑150 (MQ=255)
                      gtACGACACAGTGCGTCGCTATGGGCATGAGAGCCGCGAGAGTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGGTTATAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGttttt                                                                                                                              <  1:275064/150‑1 (MQ=255)
                      gtACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGttttt                                                                                                                              >  2:300383/1‑150 (MQ=255)
                        aCGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTAc                                                                                                                            <  1:41187/150‑1 (MQ=255)
                         cGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACTAGGGGTTATTTGTTGTTTTTACa                                                                                                                           >  2:147834/1‑150 (MQ=255)
                         cGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACa                                                                                                                           <  2:1197373/150‑1 (MQ=255)
                         cGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCACTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTTCa                                                                                                                           >  1:677349/1‑150 (MQ=255)
                           aaaCAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGa                                                                                                                         <  1:1255776/150‑1 (MQ=255)
                                gTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTcc                                                                                                                    <  2:1293530/150‑1 (MQ=255)
                                  gCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCa                                                                                                                  <  2:1328933/150‑1 (MQ=255)
                                     tCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTc                                                                                                               >  2:1136866/1‑150 (MQ=255)
                                       gCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCtt                                                                                                             <  2:889285/150‑1 (MQ=255)
                                           tGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCATAAATATCTGCCGGGTGATTTAGCCTTCCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCg                                                                                                         <  2:141368/150‑1 (MQ=255)
                                             ggCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCggg                                                                                                       <  1:1545221/150‑1 (MQ=255)
                                             ggCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCggg                                                                                                       >  2:1634696/1‑150 (MQ=255)
                                               cATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCa                                                                                                     <  1:1275325/150‑1 (MQ=255)
                                                     agTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTc                                                                                               >  2:921816/1‑150 (MQ=255)
                                                      gTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTcc                                                                                              <  1:1693576/150‑1 (MQ=255)
                                                       tGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCa                                                                                             <  1:1304277/150‑1 (MQ=255)
                                                         gCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGc                                                                                           >  2:974034/1‑150 (MQ=255)
                                                         gCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGc                                                                                           <  1:1287922/150‑1 (MQ=255)
                                                         gCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGc                                                                                           <  1:1294538/150‑1 (MQ=255)
                                                          cGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCa                                                                                          >  2:1219926/1‑150 (MQ=255)
                                                                    cACCGGCTGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGAtgct                                                                                >  1:69383/1‑150 (MQ=255)
                                                                      ccGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGAtgatgc                                                                              >  1:841255/1‑150 (MQ=255)
                                                                            gCAGAAAAATCTGCCGGGTGATTTATCCATTCAGGGGAGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCtt                                                                        >  1:285129/1‑150 (MQ=255)
                                                                             cAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTc                                                                       >  1:794306/1‑150 (MQ=255)
                                                                                           gggTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGa                                                         >  1:1551338/1‑150 (MQ=255)
                                                                                                 tttAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatca                                                   >  1:1544298/1‑150 (MQ=255)
                                                                                                  ttAGCCATTCAGGGGCGATAATCCCTTCCGATTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcac                                                  <  1:53691/150‑1 (MQ=255)
                                                                                                      ccATTCAGTGGCGATTATCCCTTCCGGTTAAAGAGCAAACCACTCAAACGAGGGGTTTTGTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatca                                              <  1:1050845/150‑1 (MQ=255)
                                                                                                      ccATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatca                                              <  2:794306/150‑1 (MQ=255)
                                                                                                      ccATTCACTGGCGAAACTCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTATGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatca                                              <  1:1056430/150‑1 (MQ=255)
                                                                                                        aTTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCaaa                                            <  1:1219926/150‑1 (MQ=255)
                                                                                                           cAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATcc                                         >  1:881567/1‑150 (MQ=255)
                                                                                                            cggggCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAACCGAGGGGTTTTTTGTTGTTTTTACAGATTACCCATTCTTGGCGGGCAACGATCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATccc                                        <  1:1335802/149‑1 (MQ=255)
                                                                                                             ggggCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATcccc                                       >  2:390287/1‑150 (MQ=255)
                                                                                                             ggggCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATcccc                                       >  1:527993/1‑150 (MQ=255)
                                                                                                              gggCGATAATCCCTTCTGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCAATCTTGGCGGGCAACGTTCCAGCTACTCGCTGCTGCCGTCTTAGTTTTGCTGTTAGAGCATCACATCAAATccccc                                      >  2:1299495/1‑149 (MQ=255)
                                                                                                                gCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCcacc                                    >  2:1122712/1‑149 (MQ=255)
                                                                                                                  gATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGATGTTCGAGCATCACATCAAATCCCCACAgg                                  >  2:1393040/1‑150 (MQ=255)
                                                                                                                   aTAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGc                                 >  1:1322659/1‑150 (MQ=255)
                                                                                                                     aaTCCGTTCAGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTCCAGATTTCCCATTCTTGTCGGGCAACGTTCGAGCAACTCGCTGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGa                               <  1:603122/150‑1 (MQ=255)
                                                                                                                      aTCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGAt                              >  1:1047278/1‑150 (MQ=255)
                                                                                                                      aTCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGAt                              >  1:1041774/1‑150 (MQ=255)
                                                                                                                      aTCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGAt                              <  1:300383/150‑1 (MQ=255)
                                                                                                                         ccTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGca                           <  1:1634696/150‑1 (MQ=255)
                                                                                                                         ccTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGca                           >  2:1085414/1‑150 (MQ=255)
                                                                                                                            tCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACAt                        <  1:655030/150‑1 (MQ=255)
                                                                                                                             ccGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATg                       <  2:1544298/150‑1 (MQ=255)
                                                                                                                              cGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTg                                                                 <  1:528751/107‑1 (MQ=255)
                                                                                                                              cGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGc                      >  2:25488/1‑150 (MQ=255)
                                                                                                                              cGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGACGTCTTCGTTTTg                                                                 >  2:528751/1‑107 (MQ=255)
                                                                                                                                   caacaaaaaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAg                 >  2:1588728/8‑150 (MQ=255)
                                                                                                                                    aaGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAgg                <  2:1540482/150‑1 (MQ=255)
                                                                                                                                     agagCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatca                                              <  1:411309/119‑1 (MQ=255)
                                                                                                                                     agagCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatca                                              >  2:411309/1‑119 (MQ=255)
                                                                                                                                     agagCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGa               >  2:769060/1‑150 (MQ=255)
                                                                                                                                       agCAAACCCCTCAAATGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACTGGCGATGCACATGCTTCAGGAct             <  2:285129/150‑1 (MQ=255)
                                                                                                                                       agCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGAct             <  2:1566034/150‑1 (MQ=255)
                                                                                                                                        gCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCCCAGGCGATGCACATGCTTCAGGActt            >  2:1426214/1‑150 (MQ=255)
                                                                                                                                         cAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGCTGCACATGCGTCAGGActtc           >  1:1452421/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATcccc                                       >  1:1094122/2‑122 (MQ=255)
                                                                                                                                         aaaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATcccc                                       <  2:1094122/121‑1 (MQ=255)
                                                                                                                                         aaaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGActtc           >  2:632756/2‑150 (MQ=255)
                                                                                                                                             ccccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTg       <  1:404290/150‑1 (MQ=255)
                                                                                                                                             ccccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTg       >  2:956960/1‑150 (MQ=255)
                                                                                                                                              cccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTATgc      >  1:429660/1‑150 (MQ=255)
                                                                                                                                               ccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCAGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATACACATGCTTCAGGACTTCTTTgcg     >  1:1180921/1‑150 (MQ=255)
                                                                                                                                                 tCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGcc   <  2:1322659/150‑1 (MQ=255)
                                                                                                                                                  cAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGccc  >  2:789916/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCC  >  NC_000913/1235539‑1235834

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