Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,235,687 A→T 9.9% intergenic (+30/+22) fadR → / ← ycgB DNA‑binding transcriptional dual regulator FadR/PF04293 family protein YcgB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,235,6870AT 9.9% 247.3 / 10.2 87intergenic (+30/+22)fadR/ycgBDNA‑binding transcriptional dual regulator FadR/PF04293 family protein YcgB
Reads supporting (aligned to +/- strand):  ref base A (37/41);  new base T (4/5);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCG  >  NC_000913/1235541‑1235831
                                                                                                                                                  |                                                                                                                                                
gAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACg                                                                                                                                               >  2:28258/1‑150 (MQ=255)
       cgcACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGtt                                                                                                                                        >  1:612010/1‑150 (MQ=255)
          aCGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCCTGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAACCCCCTCGTTTGAGGGGtttgc                                                                                                                                     <  1:764086/150‑3 (MQ=255)
          aCGATCAGGTGTACGAAACAGTGCGTCGCTATGGACATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGttttt                                                                                                                                     <  1:1449257/150‑1 (MQ=255)
           cGATCAGGTGTACGCAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGACATTCAGGGGCGATAATCCCGTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGgttttt                                                                                                                                    >  2:1538742/1‑150 (MQ=255)
           cGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGGtttgct                                                                                                                                    <  1:1327309/150‑4 (MQ=255)
           cGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGtttttt                                                                                                                                    >  1:1529274/1‑150 (MQ=255)
            gATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGGtttgctc                                                                                                                                   <  2:300932/150‑5 (MQ=255)
            gATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttg                                                                                                                                   <  1:691824/150‑1 (MQ=255)
             aTCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTttgt                                                                                                                                  >  1:911204/1‑150 (MQ=255)
             aTCAGGTGTACGAAACAGTGACTCGCTATGGGCATGAGAGTGCCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGAATTAGACATTCACGGGCGATAATCCCTTCCGTTTAAAGACCAAACCCCTCAAAAGACGCGTTTTttgt                                                                                                                                  >  2:506181/1‑150 (MQ=255)
             aTCAGGGGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGGtttgctct                                                                                                                                  <  1:142855/150‑6 (MQ=255)
                 ggTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACGGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTgttgtt                                                                                                                              >  1:686140/1‑150 (MQ=255)
                 ggTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTgttgtt                                                                                                                              >  1:508238/1‑150 (MQ=255)
                     tACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGACAGCTTTTTTCTTGtttttc                                                                                                                          >  2:582347/1‑149 (MQ=255)
                      aCGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTAc                                                                                                                         >  1:1398968/1‑150 (MQ=255)
                       cGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACa                                                                                                                        >  2:831796/1‑150 (MQ=255)
                             aGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTc                                                                                                                  <  2:2724/150‑1 (MQ=255)
                              gTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTcc                                                                                                                 >  2:93608/1‑150 (MQ=255)
                              gTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTcc                                                                                                                 <  2:1328350/150‑1 (MQ=255)
                                gCGTCTCTATGGGCCTGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTCCAGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCa                                                                                                               <  2:517624/150‑1 (MQ=255)
                                  gTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCAtt                                                                                                             <  2:706375/150‑1 (MQ=255)
                                  gTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCAtt                                                                                                             >  1:953039/1‑150 (MQ=255)
                                     gCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGTCGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTATTGTTGTTTTTACAGATTTCCCATTCtt                                                                                                          >  2:338966/1‑150 (MQ=255)
                                     gCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCtt                                                                                                          <  2:906594/150‑1 (MQ=255)
                                       tATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTgg                                                                                                        <  2:741768/150‑1 (MQ=255)
                                             cATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCa                                                                                                  <  2:950960/150‑1 (MQ=255)
                                               tGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAAc                                                                                                >  1:1208667/1‑150 (MQ=255)
                                               tGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAAc                                                                                                <  2:1398968/150‑1 (MQ=255)
                                                  gagTGGCGAGATTTGTCACCGGATGCAGAAAAATCTGCCGGGTGATTTTGCCATTCAGGGGAGATAATCCCTTCCGTTTAAAGAGCAATCGGCTCAAACGAGGTGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGtt                                                                                             <  2:368743/150‑1 (MQ=255)
                                                  gagTGGCGAGATTTGGCACGGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGtt                                                                                             <  2:686140/150‑1 (MQ=255)
                                                      ggCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAg                                                                                         >  1:276524/1‑150 (MQ=255)
                                                       gCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGc                                                                                        <  1:89594/150‑1 (MQ=255)
                                                        cGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCa                                                                                       <  1:831796/150‑1 (MQ=255)
                                                        cGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCa                                                                                       >  1:1506350/1‑150 (MQ=255)
                                                             tttGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACTTTCCAGCAACTCg                                                                                  >  2:810590/1‑150 (MQ=255)
                                                              ttGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGa                                                                                 >  2:1445930/1‑150 (MQ=255)
                                                                 gCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGAtgc                                                                              <  2:1107824/150‑1 (MQ=255)
                                                                   aCCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGAtgctg                                                                            >  2:300280/1‑150 (MQ=255)
                                                                    ccGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGAtgctgc                                                                           <  1:360990/150‑1 (MQ=255)
                                                                      ggATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCg                                                                         <  2:671467/150‑1 (MQ=255)
                                                                      ggATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCAGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAAAGTTCCAGCAACTCGATGCTGCCg                                                                         >  2:246840/1‑150 (MQ=255)
                                                                          gCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCtt                                                                     >  1:218033/1‑150 (MQ=255)
                                                                           cAGAAAAATCTGCCGGGTGATTTAGCCATTTAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTc                                                                    <  1:208711/150‑1 (MQ=255)
                                                                           cAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTc                                                                    >  1:196370/1‑150 (MQ=255)
                                                                            aGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCg                                                                   <  2:1046562/150‑1 (MQ=255)
                                                                               aaaaTCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGttt                                                                <  2:72925/150‑1 (MQ=255)
                                                                                aaaTCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGtttt                                                               <  2:630899/150‑1 (MQ=255)
                                                                                   tCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCt                                                            >  1:1427797/1‑150 (MQ=255)
                                                                                   tCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCt                                                            <  2:94554/150‑1 (MQ=255)
                                                                                    cTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTg                                                           <  1:216259/150‑1 (MQ=255)
                                                                                     tGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTATTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGt                                                          >  2:167183/1‑150 (MQ=255)
                                                                                         gggTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGa                                                      <  2:508238/150‑1 (MQ=255)
                                                                                            tGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGca                                                   <  1:1383436/150‑1 (MQ=255)
                                                                                             gATTTAGCCATTCAGGGGCGATAATCCCTTCCGTATAAAGAGCAAACCCCTCAAACGAGCGGATTTTTGATGTTTTTACAGATTTCCCATACTTGGCGGGCAACGTTCCAGCAACACGATGCTGCCGTATTCGTTTTGCTGTTCGAGcat                                                  >  2:1350447/1‑150 (MQ=255)
                                                                                               tttAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTATTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatca                                                <  1:167183/150‑1 (MQ=255)
                                                                                                ttAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcac                                               <  1:314951/150‑1 (MQ=255)
                                                                                                   gCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatc                                            <  1:15796/150‑1 (MQ=255)
                                                                                                   gCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTATTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatc                                            >  1:713347/1‑150 (MQ=255)
                                                                                                    ccATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatca                                           <  1:1074632/150‑1 (MQ=255)
                                                                                                       ttCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAAt                                        >  2:518159/1‑150 (MQ=255)
                                                                                                           gcggCGATAATCAATTCCGTTTAAATAGCAAACCCCTCAAACCAGGGGTTTATTTATCTTTTTAAAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATcccc                                    >  2:1521844/3‑150 (MQ=255)
                                                                                                                gATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAgg                               <  1:28258/150‑1 (MQ=255)
                                                                                                                gATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAgg                               <  1:1445930/150‑1 (MQ=255)
                                                                                                                 aTAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGc                              <  2:218033/150‑1 (MQ=255)
                                                                                                                   aaTCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGa                            <  1:1213535/150‑1 (MQ=255)
                                                                                                                      cccTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGc                         <  1:999182/150‑1 (MQ=255)
                                                                                                                        cTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGcac                       <  1:907755/150‑1 (MQ=255)
                                                                                                                        cTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGcac                       >  1:622097/1‑150 (MQ=255)
                                                                                                                             gTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCt                  <  1:1141356/150‑1 (MQ=255)
                                                                                                                               ttAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTc                <  1:1188331/150‑1 (MQ=255)
                                                                                                                                tAAAGAGCAAACCCCTGAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCATCTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCAGATCAAATCCCCACAGGCGATGCACATGCTTCa               <  1:498286/150‑1 (MQ=255)
                                                                                                                                 caacaaaaaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAg              >  1:912764/8‑150 (MQ=255)
                                                                                                                                  aacaaaaaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAgg             >  2:977428/7‑150 (MQ=255)
                                                                                                                                  aacaaaaaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAgg             >  2:1263098/7‑150 (MQ=255)
                                                                                                                                  aaGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAgg             >  1:1435733/1‑150 (MQ=255)
                                                                                                                                  aaGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAgg             >  2:955028/1‑150 (MQ=255)
                                                                                                                                  aaGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAgg             <  1:1406974/150‑1 (MQ=255)
                                                                                                                                    gagCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGAc           >  2:121063/1‑150 (MQ=255)
                                                                                                                                      gCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGActt         >  2:836697/1‑150 (MQ=255)
                                                                                                                                      aaaaaCCCCTCGTTTGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGActt         >  1:1266441/3‑150 (MQ=255)
                                                                                                                                        aaaCCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGActtct       >  2:123416/1‑150 (MQ=255)
                                                                                                                                            cccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTgc   <  1:93608/150‑1 (MQ=255)
                                                                                                                                            cccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGATGCCGTCTTCGTTTTGCTGTTCGAGAATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTgc   >  2:902906/1‑150 (MQ=255)
                                                                                                                                             ccTCATACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCACCGGTCCCGCAACTCGATGCTGCCGGCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTgcg  <  1:1538742/150‑1 (MQ=255)
                                                                                                                                             ccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTgcg  >  1:127619/1‑150 (MQ=255)
                                                                                                                                             ccTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTgcg  <  2:1094173/150‑1 (MQ=255)
                                                                                                                                                  aaCGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGcac                       <  1:857148/124‑1 (MQ=255)
                                                                                                                                                  aaCGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGcac                       >  2:857148/1‑124 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCG  >  NC_000913/1235541‑1235831

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