Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,188,406 C→T 13.7% E457E (GAG→GAA phoQ ← sensory histidine kinase in two‑component regulatory system with PhoP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,188,4060CT13.7% 187.2 / 18.3 73E457E (GAG→GAAphoQsensory histidine kinase in two‑component regulatory system with PhoP
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base T (6/4);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GGTATCTGCAACCGATTATAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGAT  >  NC_000913/1188259‑1188554
                                                                                                                                                   |                                                                                                                                                    
ggTATCTGCAACCGATTATAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCt                                                                                                                                                     >  7:322537/1‑149 (MQ=255)
 gTATCTGCAACCGATTATAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTAccc                                                                                                                                                      <  4:179236/147‑1 (MQ=255)
 gTATCTGCAACCGATTATAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTAccc                                                                                                                                                      <  6:264513/147‑1 (MQ=255)
    tCTGCAACCGATTTTAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACAAAAAATCACCTCCATCCGCGCACAGCCCAGCATGCTCTCTCCGGCGACGATTTTAAGCTCata                                                                                                                                                 >  1:272071/1‑149 (MQ=255)
      tGCAACCGATTATAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCTTCATAtt                                                                                                                                               >  7:160931/1‑149 (MQ=255)
      tGCAACCGATTATAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTAATCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCACCATGCTCTCTCAGGCGACGATTTTACCCTCAAAtt                                                                                                                                               >  7:197279/1‑149 (MQ=255)
       gCAACCGATTATAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTg                                                                                                                                              >  6:11034/1‑149 (MQ=255)
         aaCCGATTATAACGGAGGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCCGCATGCGCTCTCCGGCGACGAGTTTACCCTCATATTGCt                                                                                                                                            >  6:101077/1‑149 (MQ=255)
           ccGATTATAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCTTCATATTGCTCg                                                                                                                                          <  4:54784/149‑1 (MQ=255)
           ccGATTATAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCAGGCGACGATTTTACCTTCATATTGCTCg                                                                                                                                          >  4:234261/1‑149 (MQ=255)
            cGATTATAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGTCGATTTTACCCTCATATTGCTCgg                                                                                                                                         >  8:155245/1‑149 (MQ=255)
                       gATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGAATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCgg                                                                                                                              <  4:175432/149‑1 (MQ=255)
                              aaCGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGTCGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTAcc                                                                                                                        <  7:155245/148‑1 (MQ=255)
                               aCGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGc                                                                                                                      <  8:75034/149‑1 (MQ=255)
                                  tAATGGGTGAAGTATGGGCATTTTTATTCATCTTTCGGCGGAGTATGCTGGGGTACAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAg                                                                                                                   <  6:89884/149‑1 (MQ=255)
                                  tAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAg                                                                                                                   <  8:322537/149‑1 (MQ=255)
                                        gTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACTCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTa                                                                                                              >  8:21509/1‑148 (MQ=255)
                                           aaGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCTTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTa                                                                                                              <  5:78581/145‑1 (MQ=255)
                                           aaGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCTTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTa                                                                                                              >  6:78581/1‑145 (MQ=255)
                                           aaGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCTTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACAcc                                                                                                          >  5:438506/1‑149 (MQ=255)
                                           aaGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTacac                                                                                                           <  6:174705/148‑1 (MQ=255)
                                             gTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTcgcg                                                                                                                               >  8:45697/1‑126 (MQ=255)
                                             gTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTcgcg                                                                                                                               <  7:45697/126‑1 (MQ=255)
                                                 gggCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGcc                                                                                                     >  4:253349/1‑148 (MQ=255)
                                                     atatTTATTCATATTTCTGCGCAGAATGCTGTCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTACGGCGACGATTTTACACTCATATTGCTCGGTGATTTCGCTGGCTACAGCCAGCCCTACACCTTGCCCAggg                                                                                                <  2:272071/149‑1 (MQ=255)
                                                          tattCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACATCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTAACCTTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTa                                                                                            >  4:299127/1‑148 (MQ=255)
                                                             tCATCTTTCGGCGCAGAATGCTGGCGACCGAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAg                                                                                         >  8:343514/1‑148 (MQ=255)
                                                               aTCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTAt                                                                                      >  3:399969/1‑149 (MQ=255)
                                                                 cTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGcc                                                                                                     <  5:336279/132‑1 (MQ=255)
                                                                 cTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGcc                                                                                                     >  6:336279/1‑132 (MQ=255)
                                                                       gcgcAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCTTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGt                                                                              <  1:461899/149‑1 (MQ=255)
                                                                           aGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGAc                                                                          <  7:185250/149‑1 (MQ=255)
                                                                                      gACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGaa                                                                >  5:342222/1‑148 (MQ=255)
                                                                                      gACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGaa                                                                >  3:328951/1‑148 (MQ=255)
                                                                                      gACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCAGGCGACGATTTTACCTTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGaaa                                                               <  3:234261/149‑1 (MQ=255)
                                                                                          aaaaaTCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGa                                                           <  4:399969/149‑1 (MQ=255)
                                                                                            aaaTCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGAcc                                                         <  3:293366/149‑1 (MQ=255)
                                                                                             aaTCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACct                                                        <  8:300519/149‑1 (MQ=255)
                                                                                                   cTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgct                                                  >  7:132235/1‑149 (MQ=255)
                                                                                                      cATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgc                                               <  6:342222/149‑1 (MQ=255)
                                                                                                      cATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgc                                               <  8:132235/149‑1 (MQ=255)
                                                                                                       aTCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTc                                                                   >  5:156006/1‑128 (MQ=255)
                                                                                                       aTCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTc                                                                   <  6:156006/128‑1 (MQ=255)
                                                                                                       aTCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgct                                              >  8:72956/1‑149 (MQ=255)
                                                                                                          cgcgCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGAGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTaa                                           >  2:272023/1‑149 (MQ=255)
                                                                                                           gcgcACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAAt                                          >  4:59021/1‑149 (MQ=255)
                                                                                                                  gCCCAGCATGCTCTGTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGTGGGCTACCGCCAGACCTACACCTTACCAAGTGCGTAATGTATCAACCCGTTGACAACGGTCGAATATGACCTCTCGCTTGCTTAATGGAATAc                                   <  8:197279/149‑1 (MQ=255)
                                                                                                                    ccAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgc                                               <  1:104972/135‑1 (MQ=255)
                                                                                                                    ccAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCgcttgc                                               >  2:104972/1‑135 (MQ=255)
                                                                                                                     cAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCCGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATTCCCTCTCGCTTGCTTAATGGAATACCgg                                >  6:5092/1‑149 (MQ=255)
                                                                                                                     cAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACgg                                                                     >  2:183831/1‑112 (MQ=255)
                                                                                                                     cAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACgg                                                                     <  1:183831/112‑1 (MQ=255)
                                                                                                                     cAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCgg                                >  1:356803/1‑149 (MQ=255)
                                                                                                                      atcATGCTCTCTCCGGCGACGATTTTACCATCATATTGCTCGGTGATTTAGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCggg                               <  2:356803/147‑1 (MQ=255)
                                                                                                                      aGCATGCTCTCTCCGGCGACGATTTTACTCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCggg                               <  7:21509/149‑1 (MQ=255)
                                                                                                                          tGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCca                           >  3:370964/1‑149 (MQ=255)
                                                                                                                           gCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCcat                          >  3:268417/1‑149 (MQ=255)
                                                                                                                            ctctctCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCcat                          >  1:175127/1‑148 (MQ=255)
                                                                                                                             tctctcCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGGCCTCTCGCTTGCTTAATGGAATACCGGGGCcatca                        >  2:389017/1‑149 (MQ=255)
                                                                                                                             tctctcCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAAGGACCTCTCGCTTGCTTAATGGAATACCGGGGCcat                          >  8:419466/1‑147 (MQ=255)
                                                                                                                               tctcCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCcatca                        <  3:59021/147‑1 (MQ=255)
                                                                                                                                   cGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTc                   <  4:370964/148‑1 (MQ=255)
                                                                                                                                   cGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGAGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTc                   <  1:272023/148‑1 (MQ=255)
                                                                                                                                    ggCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGaaa                                                               <  1:291316/103‑1 (MQ=255)
                                                                                                                                    ggCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGaaa                                                               >  2:291316/1‑103 (MQ=255)
                                                                                                                                    ggCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGa                 <  5:11034/149‑1 (MQ=255)
                                                                                                                                      cgacgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGAcc               <  5:101077/149‑1 (MQ=255)
                                                                                                                                       gacgaTTTTACCTTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGACca              >  1:222761/1‑149 (MQ=255)
                                                                                                                                       gacgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACCCCTTGCCCAGGGCGTAAAGTATCAACACGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAAAGGCATACCGGGGCCATCATCCTCGACca              >  8:326419/1‑149 (MQ=255)
                                                                                                                                       gacgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGACca              >  2:5822/1‑149 (MQ=255)
                                                                                                                                          gaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTg                                                                            >  5:153046/1‑84 (MQ=255)
                                                                                                                                          gaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTg                                                                            <  6:153046/84‑1 (MQ=255)
                                                                                                                                          gaTGTTACCCTCATATTGCTCGGGGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGTTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGAGGCCATCATCCTCGACCACaa           <  3:99976/149‑1 (MQ=255)
                                                                                                                                                aCCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGACCACAATATag      <  2:175127/148‑1 (MQ=255)
                                                                                                                                                aCCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGACCACAATATag      <  1:5822/148‑1 (MQ=255)
                                                                                                                                                 cccTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGGCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGACCACAATATagag    <  3:253349/149‑1 (MQ=255)
                                                                                                                                                  ccTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAAc                                                                                  >  8:236322/1‑70 (MQ=255)
                                                                                                                                                  ccTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAAc                                                                                  <  7:236322/70‑1 (MQ=255)
                                                                                                                                                  ccTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGATGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCggg                               >  3:144116/1‑121 (MQ=255)
                                                                                                                                                  ccTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGATGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCggg                               <  4:144116/121‑1 (MQ=255)
                                                                                                                                                   cTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGAt  <  4:328951/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGTATCTGCAACCGATTATAACGGATGCTTAACGTAATGCGTGAAGTATGGGCATATTTATTCATCTTTCGGCGCAGAATGCTGGCGACCAAAAATCACCTCCATCCGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGAT  >  NC_000913/1188259‑1188554

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: