Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,235,687 A→T 9.4% intergenic (+30/+22) 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,6870AT 9.4% 210.8 / 8.7 83intergenic (+30/+22)fadR/ycgBfatty acid metabolism regulon transcriptional regulator/SpoVR family stationary phase protein
Reads supporting (aligned to +/- strand):  ref base A (38/37);  new base T (3/5);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCC  >  NC_000913/1235538‑1235834
                                                                                                                                                     |                                                                                                                                                   
aGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCaa                                                                                                                                                     <  1:488159/150‑1 (MQ=255)
 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                                                                                           <  1:1287922/150‑1 (MQ=255)
                                                          gCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGc                                                                                           >  2:974034/1‑150 (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                                       >  1:527993/1‑150 (MQ=255)
                                                                                                              ggggCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATcccc                                       >  2:390287/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:300383/150‑1 (MQ=255)
                                                                                                                       aTCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGAt                              >  1:1041774/1‑150 (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                                              >  2:411309/1‑119 (MQ=255)
                                                                                                                                      agagCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatca                                              <  1:411309/119‑1 (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                                       <  2:1094122/121‑1 (MQ=255)
                                                                                                                                          aaaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATcccc                                       >  1:1094122/2‑122 (MQ=255)
                                                                                                                                          aaaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGActtc           >  2:632756/2‑150 (MQ=255)
                                                                                                                                              ccccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTg       >  2:956960/1‑150 (MQ=255)
                                                                                                                                              ccccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTg       <  1:404290/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                   
AGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCC  >  NC_000913/1235538‑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: