Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,235,687 A→T 9.0% 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.0% 223.2 / 8.2 77intergenic (+30/+22)fadR/ycgBfatty acid metabolism regulon transcriptional regulator/SpoVR family stationary phase protein
Reads supporting (aligned to +/- strand):  ref base A (40/30);  new base T (3/4);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

TGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCC  >  NC_000913/1235540‑1235835
                                                                                                                                                   |                                                                                                                                                     
tGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTAAGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCaaa                                                                                                                                                      >  8:179254/1‑149 (MQ=255)
  aaGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAAc                                                                                                                                                     >  3:231181/1‑148 (MQ=255)
   aGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGTTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGa                                                                                                                                                   <  2:183115/149‑1 (MQ=255)
   aGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCGCTTCCGTTTAAAGAGCAAACCCCTCAAAc                                                                                                                                                     >  1:179488/1‑147 (MQ=255)
   aGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCTAACCCCTCAAACGa                                                                                                                                                   <  8:116766/149‑1 (MQ=255)
   aGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGa                                                                                                                                                   <  2:172302/149‑1 (MQ=255)
         gcACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGG‑tt                                                                                                                                            <  3:191729/149‑1 (MQ=255)
           aCGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑tttt                                                                                                                                          >  5:308386/1‑149 (MQ=255)
            cGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑ttttt                                                                                                                                         <  7:138013/149‑1 (MQ=255)
              aTCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCGTTTGAGGGG‑tttgctc                                                                                                                                       <  3:368773/149‑5 (MQ=255)
               tCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTttgt                                                                                                                                      >  5:89358/1‑149 (MQ=255)
                       aagAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTa                                                                                                                              >  2:9618/3‑149 (MQ=255)
                        cGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTa                                                                                                                              >  5:149381/1‑148 (MQ=255)
                        cGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTa                                                                                                                              <  7:183490/148‑1 (MQ=255)
                         gAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTa                                                                                                                              >  6:27120/1‑147 (MQ=255)
                           aaCAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑tttt                                                                                                                                          >  2:4572/1‑133 (MQ=255)
                           aaCAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑tttt                                                                                                                                          <  1:4572/133‑1 (MQ=255)
                                tGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTcc                                                                                                                     >  8:274761/1‑149 (MQ=255)
                                     cGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTc                                                                                                                >  5:290382/1‑149 (MQ=255)
                                       cTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCtt                                                                                                              >  4:333253/1‑149 (MQ=255)
                                          tGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCtt                                                                                                              >  5:266536/1‑146 (MQ=255)
                                          tGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCtt                                                                                                              <  6:266536/146‑1 (MQ=255)
                                          tGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGc                                                                                                           >  3:188229/1‑149 (MQ=255)
                                            ggCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCgg                                                                                                         <  7:274761/149‑1 (MQ=255)
                                             gCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGc                                                                                                           <  6:399332/146‑1 (MQ=255)
                                             gCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGc                                                                                                           >  5:399332/1‑146 (MQ=255)
                                             gCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCggg                                                                                                        <  1:9618/149‑1 (MQ=255)
                                             gCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCggg                                                                                                        >  2:352868/1‑149 (MQ=255)
                                                      tGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTcc                                                                                               <  6:308386/149‑1 (MQ=255)
                                                        gCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGGGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAg                                                                                             <  2:24327/149‑1 (MQ=255)
                                                         cGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAg                                                                                             >  7:356730/1‑148 (MQ=255)
                                                             aTTTGGCAGCGGATGCAGAAATATCTGCCGGTTGATTTAGCGATTCAGGGGCGATAATCCCTTCAGTTTAAAGAGCTAACCCCTCAAACGAGTGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAAc                                                                                         <  8:92913/148‑1 (MQ=255)
                                                              tttGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTc                                                                                       >  1:301046/1‑149 (MQ=255)
                                                                 ggCACCGTATGCAGAAAAATCTGCCGGGTTATTTAGCCATTCAGGGGCTATAATCCCTTCCGTTTAAACCGCAAACCCCACAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGa                                                                                     <  2:179488/148‑1 (MQ=255)
                                                                  gCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGAt                                                                                    <  7:179254/148‑1 (MQ=255)
                                                                    aCCGGATGCAGAAAAAACAGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TGTTTTGTTGATTATACAGCGTTCCCATTCTTGGCGGGCAACGTTCCAGCACCTCGAtgcg                                                                                 >  6:304061/1‑148 (MQ=255)
                                                                             aGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTTCCGTCTTc                                                                        >  7:83410/1‑149 (MQ=255)
                                                                             aGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTc                                                                        <  5:210338/149‑1 (MQ=255)
                                                                             aGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTc                                                                        >  7:238055/1‑149 (MQ=255)
                                                                              gAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCg                                                                       >  2:19473/1‑149 (MQ=255)
                                                                                aaaaTCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGtt                                                                     <  6:149381/149‑1 (MQ=255)
                                                                                    tCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGc                                                                 >  3:330892/1‑149 (MQ=255)
                                                                                      tGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTg                                                               >  3:255306/1‑149 (MQ=255)
                                                                                        ccGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTc                                                                                                <  5:63018/114‑1 (MQ=255)
                                                                                        ccGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTc                                                                                                >  6:63018/1‑114 (MQ=255)
                                                                                           ggTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGa                                                          <  4:255306/149‑1 (MQ=255)
                                                                                               aTTTAGCCATTCAGGGGCGATAATCACTTCCGTTTAAAGAGCAAACCCCTCAAACGAGCGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcat                                                      >  6:375496/1‑149 (MQ=255)
                                                                                                 ttAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTg                                                                  >  5:46550/1‑135 (MQ=255)
                                                                                                 ttAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTg                                                                  <  6:46550/135‑1 (MQ=255)
                                                                                                       aTTCATGGTCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGcc                                                                              <  8:294301/117‑1 (MQ=255)
                                                                                                       aTTCAGGGGCGATAATCCCTTCCGTTTAAAGATCAAACCCCTCAAACGATGGG‑TTTTTTGGTGTTTTTACAGATTTCCCAGTCTTGGCGGGCAACGTTCCAGCAACTAGATGCTGcc                                                                              >  7:294301/1‑117 (MQ=255)
                                                                                                          cAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATc                                           >  6:1246/1‑149 (MQ=255)
                                                                                                           aGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTAACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATcc                                          >  6:300891/1‑149 (MQ=255)
                                                                                                            ggggCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATccc                                         <  8:356730/149‑1 (MQ=255)
                                                                                                               gCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCcac                                      <  4:231181/149‑1 (MQ=255)
                                                                                                                  aTAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAgg                                   >  6:380282/1‑149 (MQ=255)
                                                                                                                      tCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacat                                                 <  8:65541/131‑1 (MQ=255)
                                                                                                                      tCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacat                                                 >  7:65541/1‑131 (MQ=255)
                                                                                                                       cccTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGa                                >  1:205757/1‑147 (MQ=255)
                                                                                                                       cccTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGa                                <  2:205757/147‑1 (MQ=255)
                                                                                                                        ccTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATg                              <  5:380282/148‑1 (MQ=255)
                                                                                                                            ccGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGcac                           <  4:188229/147‑1 (MQ=255)
                                                                                                                                   aaGAGCAATCCCCTCAAACGAGGGG‑TTTTTTGTTGTTGTTACAGAATTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATACCCACAGGCGATGCACATGCTTCAg                  <  6:290382/149‑1 (MQ=255)
                                                                                                                                     gagCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCACCTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGa                >  5:151656/1‑149 (MQ=255)
                                                                                                                                     gagCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGa                >  5:157687/1‑149 (MQ=255)
                                                                                                                                     gagCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTATTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGa                >  7:182489/1‑149 (MQ=255)
                                                                                                                                      aaaaaaCCCCTCGTTTGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatc                                                <  6:128837/113‑1 (MQ=255)
                                                                                                                                      aaaaaaCCCCTCGTTTGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGcatcacatc                                                >  5:128837/4‑116 (MQ=255)
                                                                                                                                        aaaaCCCCTCGTTTGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGc                                                                 >  6:399745/2‑97 (MQ=255)
                                                                                                                                        aaaaCCCCTCGTTTGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGc                                                                 <  5:399745/96‑1 (MQ=255)
                                                                                                                                        aaaaCCCCTCGTTTGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGAct              >  2:42975/2‑148 (MQ=255)
                                                                                                                                         aaaCCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCACCTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGActtc            <  1:83666/149‑1 (MQ=255)
                                                                                                                                          aaCCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGActtct           >  8:227048/1‑149 (MQ=255)
                                                                                                                                           aCCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGActtctt          >  3:108443/1‑149 (MQ=255)
                                                                                                                                            ccccTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCttt         <  8:377397/149‑1 (MQ=255)
                                                                                                                                             cccTCAAACGAGGGGTTTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCttt         >  4:383601/1‑149 (MQ=255)
                                                                                                                                                tCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTgcgc     >  6:156186/1‑149 (MQ=255)
                                                                                                                                                 cAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGcc    >  5:389178/1‑149 (MQ=255)
                                                                                                                                                   aaCGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGcccc  <  4:330892/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGG‑TTTTTTGTTGTTTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCC  >  NC_000913/1235540‑1235835

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