Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,291,284 C→A 5.0% intergenic (+137/+74) narP → / ← ccmH response regulator in two‑component regulatory system with NarQ/heme lyase, CcmH subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,291,2840CA 5.0% 337.3 / 4.6 99intergenic (+137/+74)narP/ccmHresponse regulator in two‑component regulatory system with NarQ/heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base C (47/47);  new base A (2/3);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

CACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCACCGCACGGCGAGTA  >  NC_000913/2291140‑2291430
                                                                                                                                                |                                                                                                                                                  
cacaATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAGGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATGCGGCGTGAAcgc                                                                                                                                                >  3:76179/1‑149 (MQ=37)
   aaTAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATGCGGCGTGAAcgcctt                                                                                                                                             <  5:324614/149‑3 (MQ=255)
   aaTAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                              >  5:225136/1‑148 (MQ=255)
   aaTAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                              >  3:315457/1‑148 (MQ=255)
        aaaTAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTTGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                        <  6:52578/149‑1 (MQ=255)
          aTAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                      <  7:127932/149‑1 (MQ=255)
          aTAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                      <  7:132911/149‑1 (MQ=255)
          aTAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCAGTATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                      <  6:248674/149‑1 (MQ=255)
           tAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                     <  2:6993/149‑1 (MQ=255)
           tAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                     <  7:18313/149‑1 (MQ=255)
             gCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                    >  4:97323/1‑148 (MQ=255)
              cccGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATa                                                                                                                                   >  8:55789/1‑148 (MQ=255)
              cccGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATa                                                                                                                                   >  8:50320/1‑148 (MQ=255)
               ccGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATaaa                                                                                                                                 >  8:148632/1‑149 (MQ=255)
               ccGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATaaa                                                                                                                                 <  7:229856/149‑1 (MQ=255)
                 gATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCATACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAAc                                                                                                                                >  6:298557/1‑148 (MQ=255)
                     gATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTacaaaatt                                                                                                                            <  4:76179/148‑8 (MQ=21)
                     gATTTATAATAGGTCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTTGATCGAGCAATTTATATAATTTACACAATTTTGTAGTCAGGATAAGTCTTTCACGACGAAGACGTAATAAACAAAgc                                                                                                                           >  6:315384/1‑149 (MQ=17)
                          atcatcGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCAc                                                                                                                       <  7:7288/148‑1 (MQ=255)
                           tcatcGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACtt                                                                                                                     >  7:17713/1‑149 (MQ=255)
                             atcGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACttt                                                                                                                    <  1:162962/148‑1 (MQ=255)
                                  gCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAAc                                                                                                              <  3:75191/149‑1 (MQ=255)
                                   cTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACa                                                                                                             <  3:322647/149‑1 (MQ=255)
                                    tGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAg                                                                                                            >  8:218163/1‑149 (MQ=255)
                                    tGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAg                                                                                                            >  3:182626/1‑149 (MQ=255)
                                     gagaTTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGt                                                                                                           >  5:195399/1‑149 (MQ=255)
                                      agaTTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAg                                                                                                            <  5:109058/147‑1 (MQ=255)
                                      agaTTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGt                                                                                                           <  2:269641/148‑1 (MQ=255)
                                       gaTTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTc                                                                                                          >  6:339899/1‑148 (MQ=255)
                                            aTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGccc                                                                                                    >  2:61270/1‑149 (MQ=255)
                                            aTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGcc                                                                                                     <  2:202792/148‑1 (MQ=255)
                                              gACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                     <  6:153682/114‑2 (MQ=18)
                                              gACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                     >  5:153682/1‑113 (MQ=18)
                                                cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATg                                                                                                >  3:284664/1‑149 (MQ=255)
                                                  aaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                    <  8:3552/111‑1 (MQ=255)
                                                  aaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                    >  7:3552/1‑111 (MQ=255)
                                                  aaCGCAAAACTGCCTGATGCGCTACGATTATCAGGCCGAAGTCGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATTAGGCGTTCACGCCGCATCCGGCATTAACTAAGCGCACTTTGTCAACAGTCTAGCCCTATGGc                                                                                              <  6:197476/149‑1 (MQ=18)
                                                     gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATc                                                                                           <  7:148632/149‑1 (MQ=255)
                                                         aaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGAGCACTTTGTCAACACTCTAGCCCGATGGCATCACCa                                                                                       >  7:143886/1‑149 (MQ=255)
                                                          aCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCAt                                                                                      <  3:305784/149‑1 (MQ=255)
                                                             gCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCgg                                                                                   >  8:194994/1‑149 (MQ=255)
                                                              ccTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCgg                                                                                   <  6:263406/148‑1 (MQ=255)
                                                                 gATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCa                                                                                                                <  1:26163/116‑1 (MQ=255)
                                                                 gATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCa                                                                                                                >  2:26163/1‑116 (MQ=255)
                                                                 gATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGAGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCct                                                                               >  5:279758/1‑149 (MQ=255)
                                                                   tGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCctct                                                                             >  2:16555/1‑149 (MQ=255)
                                                                        tACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACa                                                                                                             >  3:16261/1‑112 (MQ=255)
                                                                        tACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACa                                                                                                             <  4:16261/112‑1 (MQ=255)
                                                                          cGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTtttat                                                                       <  4:127138/148‑1 (MQ=255)
                                                                            cTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTtttattta                                                                    <  3:97323/149‑1 (MQ=255)
                                                                            cTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTtttattta                                                                    <  2:173918/149‑1 (MQ=255)
                                                                             ttATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTAc                                                                   >  7:103684/1‑149 (MQ=255)
                                                                               aTCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTAc                                                                   <  5:298557/147‑1 (MQ=255)
                                                                               aTCAGGCCTACGTGGAGCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGAACCTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCAGCACCATCGGGCCTCTTTTTATTTActc                                                                 >  6:32995/1‑149 (MQ=255)
                                                                                  aGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTc                                                                                                                 <  1:274121/98‑1 (MQ=255)
                                                                                  aGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTc                                                                                                                 >  2:274121/1‑98 (MQ=255)
                                                                                  aGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTcc                                                              >  7:200314/1‑149 (MQ=255)
                                                                                      cTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTgcg                                                          >  7:38955/1‑149 (MQ=255)
                                                                                       tACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTgcgg                                                         >  8:260443/1‑149 (MQ=255)
                                                                                       tACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTgcgg                                                         >  3:113732/1‑149 (MQ=255)
                                                                                         cGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCTCTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTgcgg                                                         <  4:213802/147‑1 (MQ=255)
                                                                                           tGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGa                                                      >  3:148700/1‑148 (MQ=255)
                                                                                            ggATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGAc                                                     >  2:266933/1‑148 (MQ=255)
                                                                                              atcgatcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACaaa                                                  >  6:40267/1‑149 (MQ=255)
                                                                                               tcgatcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAAt                                                 <  6:141765/149‑1 (MQ=255)
                                                                                               tcgatcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAAt                                                 <  7:218163/149‑1 (MQ=255)
                                                                                                 gatcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATg                                                <  7:50320/148‑1 (MQ=255)
                                                                                                 gatcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATg                                                <  7:55789/148‑1 (MQ=255)
                                                                                                 gatcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATg                                                <  4:225807/148‑1 (MQ=255)
                                                                                                 gatcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGt                                               <  4:284664/149‑1 (MQ=255)
                                                                                                  atcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGtt                                              >  2:276709/1‑149 (MQ=255)
                                                                                                  atcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGtt                                              <  6:147726/149‑1 (MQ=255)
                                                                                                   tcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTg                                             <  6:77893/149‑1 (MQ=255)
                                                                                                          aTTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGcc                                      >  8:109035/1‑149 (MQ=255)
                                                                                                           ttGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCgg                                                                                   >  5:218951/1‑103 (MQ=255)
                                                                                                           ttGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCgg                                                                                   <  6:218951/103‑1 (MQ=255)
                                                                                                              aaTTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTgcg                                  >  5:245457/1‑149 (MQ=255)
                                                                                                                  tACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGa                               <  1:61270/148‑1 (MQ=255)
                                                                                                                  tACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGa                               >  8:163530/1‑148 (MQ=255)
                                                                                                                   acacAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAta                             >  4:301968/1‑149 (MQ=255)
                                                                                                                    cacaATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAta                             <  4:182626/148‑1 (MQ=255)
                                                                                                                    cacaATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAta                             <  5:337545/148‑1 (MQ=255)
                                                                                                                    cacaATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAta                             <  5:339899/148‑1 (MQ=255)
                                                                                                                     acaATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAtact                           <  2:356521/149‑1 (MQ=255)
                                                                                                                      caATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAtacta                          >  3:258564/1‑149 (MQ=255)
                                                                                                                       aaTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCCCCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGAtactac                         >  1:312842/1‑149 (MQ=255)
                                                                                                                         ttttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGt                       >  2:186115/1‑149 (MQ=255)
                                                                                                                           ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGtt                      <  1:266933/148‑1 (MQ=255)
                                                                                                                              tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCaa                   <  8:143886/148‑1 (MQ=255)
                                                                                                                              tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAAt                  <  3:162030/149‑1 (MQ=255)
                                                                                                                                 gCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATcgc               >  3:93133/1‑147 (MQ=255)
                                                                                                                                 gCCGGATAAGGCGTTCACGACGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCAc               >  8:208284/1‑149 (MQ=255)
                                                                                                                                  ccGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCAcc              <  5:40267/149‑1 (MQ=255)
                                                                                                                                   cGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCACCg             >  4:76363/1‑149 (MQ=255)
                                                                                                                                      aTAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCACCGCAc          >  2:194984/1‑149 (MQ=255)
                                                                                                                                       tAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCACCGCACg         <  6:195399/149‑1 (MQ=255)
                                                                                                                                            ctttCACGCCGCATCCTGCATAAACTAAGCGCACTTTGTCAACTGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCACCGCACGGCGAg    <  5:97498/147‑1 (MQ=255)
                                                                                                                                              ttCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCACCGCACGGCGAGTa  >  1:353722/1‑149 (MQ=255)
                                                                                                                                              ttCACGCCGCATCCGGCATAAACTAAGAGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTGTTTATTTACTCTCCTGCGGAGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCCCCGCACGGCGAgga  >  7:288806/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAATGTTGCATCGCCTGCGCGATACTACGTTCAATCACCGCACGGCGAGTA  >  NC_000913/2291140‑2291430

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: