Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,291,234 C→G 5.0% intergenic (+89/+122) 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,2340CG 5.0% 239.6 / 5.8 80intergenic (+89/+122)narP/ccmHresponse regulator in two‑component regulatory system with NarQ/heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base C (36/40);  new base G (2/2);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-02

TCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAA  >  NC_000913/2291090‑2291380
                                                                                                                                                |                                                                                                                                                  
tCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAAGTGgatcgatc                                                                                                                                                >  1:1066339/1‑149 (MQ=255)
  cGCTCACGCGTGGCGGCCACCATTCTGTTCCGGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGaatcgatcca                                                                                                                                              >  2:537905/1‑147 (MQ=255)
       aCGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCCTCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGGTGCGCTACGCCTAGCAGGCCTACGTGGATCGATCAATTTAt                                                                                                                                         >  2:373678/1‑149 (MQ=255)
       aCGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTAt                                                                                                                                         <  1:985623/149‑1 (MQ=255)
        cgcgTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTAtt                                                                                                                                        <  1:195699/149‑1 (MQ=255)
         gcgTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTTTCAGGCCTCCGTGGATCGATCATTTTAttt                                                                                                                                       >  1:196206/1‑148 (MQ=255)
         gcgTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGAGTTATGACAACCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAAGTGGAGCGATCCATTAATTg                                                                                                                                       >  2:116572/1‑149 (MQ=38)
              gcggcCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGAGGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAAttt                                                                                                                                  >  2:618278/1‑149 (MQ=255)
               cggcCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGTATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACCGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTa                                                                                                                                 <  2:105729/149‑1 (MQ=255)
                 gcCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTaca                                                                                                                               >  2:906802/1‑149 (MQ=255)
                  ccaccaTTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCACTTTATTGAATTTacac                                                                                                                              >  1:711592/1‑149 (MQ=255)
                  ccaccaTTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTacac                                                                                                                              >  1:138545/1‑149 (MQ=255)
                        tgcTGTTCCGTCAATAACGCGGGGCACAATAAAAATAACCAGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCGTGTGGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAAtttt                                                                                                                        <  2:500519/147‑1 (MQ=38)
                          cTGTTCCGGCAACAACGCGGGTCACAATAAAAAGAGCCCGATCGATTTATCATCGGGCTGAGTTTTATGACAACCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGt                                                                                                                      <  2:344280/149‑1 (MQ=38)
                           tGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTa                                                                                                                     <  1:369908/149‑1 (MQ=255)
                           tGTTCCGGCAACAACGCGGCGCACAATAAAAATAGCCCTATGTATTTATCATCTGGCTGAGATTTATGACAAACGCAAAACGGCCTGATGCGCTACGCGTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTa                                                                                                                     <  1:523338/149‑1 (MQ=255)
                            gTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCGAGCAATTGATTGAATTTACACAATTTTGTAg                                                                                                                    <  1:280251/149‑1 (MQ=255)
                            gTTCCTGCAACAACGCGGGGCACAAGAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAg                                                                                                                    >  1:688333/1‑149 (MQ=255)
                             ttCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAgg                                                                                                                   >  1:105088/1‑149 (MQ=255)
                               ccTGCAACAACGCGGGGCACGATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGGTTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACACTTTTGTACGcc                                                                                                                 >  2:866446/1‑149 (MQ=255)
                               ccTGCAACAACGCGGGGCACGATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGcc                                                                                                                 >  2:629496/1‑149 (MQ=255)
                                cTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCATGCCTACGTGGATGGATCAATTTATTGAATTTACACAATTTTGTAGGCCg                                                                                                                >  1:737152/1‑149 (MQ=255)
                                cTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCATGCCTACGTGGATGGATCAATTTATTGAATTTACACAATTTTGTAGGCCg                                                                                                                >  1:169690/1‑149 (MQ=255)
                                  gCAACAACGCGGGGCACAATAAAAATAGCCCGATGGAATTAACATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCGTATCAGGCCTACGTGGATCGNTCAATTTAGTGAATTTACACAATTTAGTGGGCCggg                                                                                                              >  2:787398/1‑148 (MQ=37)
                                   cattcacGCGGGGCCTAATAAACATAGTCTGATCGATCTATCAGCGGGCTGAGATTTATGACATAAGCAAAACTTCATGATGCGCTACGCTTATCAGGCCTACGTGGTTCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGac                                                                                                             <  2:140314/144‑2 (MQ=21)
                                      caacGCGGGGCACAATAACAATAGCCCGGTGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAg                                                                                                          <  2:550154/149‑1 (MQ=255)
                                      caacGCGGGGCACAATAAAAATAGCCCGAAGGATTTATCATCGGGCTGAGATTTAAGACAAACGCAAAACTGTCTGATGCGCTACGCTTATCAGGACTACGTGGCGCGATCAATTTATGGAATTTGCACAATATAGGAGGACggcgacg                                                                                                          >  1:724417/1‑144 (MQ=25)
                                        acGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGc                                                                                                        >  1:795960/1‑149 (MQ=255)
                                             gggCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCa                                                                                                   <  1:1126867/149‑1 (MQ=255)
                                                 acaATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                               <  1:541735/149‑1 (MQ=255)
                                                   aaTAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGGTCAcgccgc                                                                                             >  1:1145215/1‑149 (MQ=255)
                                                      aaaaaTAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                          <  2:1095199/149‑1 (MQ=255)
                                                           tAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCATGCCTACGTGGATGGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                     <  2:169690/149‑1 (MQ=255)
                                                           tAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCATGCCTACGTGGATGGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                     <  2:737152/149‑1 (MQ=255)
                                                            aGCCCGATGGATTTATCATCGGGCTGAGATGTATGACAAACGCAAAACTGCCTGAGGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                    >  2:223838/1‑149 (MQ=255)
                                                              cccGATGGATTTATCATCGGGCTGAGATTTATGATTAATGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCCGATAAGGCGTTCACGCCGCATCCGGCATaa                                                                                  <  1:601655/149‑1 (MQ=255)
                                                              cccGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCGATTTATGGAATTTACACCATTTTGTAGGCGGGATCAGGCGTTGACGCGGCATCCGCGATaa                                                                                  >  2:215939/1‑149 (MQ=21)
                                                               ccGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACCGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATaaa                                                                                 <  2:835976/149‑1 (MQ=255)
                                                                  aTGGTTTTAACACCGGGCTGAGATTTATTACAAACACAAAACTGCCGGTGGGGCTACGCTTATCCGGCCTACGTGGATCGATCAATTTATTGAAGTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTa                                                                              <  2:581699/149‑1 (MQ=21)
                                                                   gggATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTaa                                                                             <  2:1062220/148‑1 (MQ=255)
                                                                    ggATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAg                                                                            >  2:856184/1‑149 (MQ=255)
                                                                     gATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAgc                                                                           <  1:755157/149‑1 (MQ=255)
                                                                     gATTTATCATCGGGATGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAgc                                                                           <  1:560982/149‑1 (MQ=255)
                                                                       tttATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAgcgc                                                                         >  1:945194/1‑149 (MQ=255)
                                                                        ttctcatcGGGCGGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCa                                                                        <  2:838771/146‑1 (MQ=255)
                                                                        ttATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCa                                                                        <  2:832874/149‑1 (MQ=255)
                                                                           tcatcGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACtt                                                                     <  1:393568/149‑1 (MQ=255)
                                                                           tcatcGGGCTGAGAATTATGACAAACGCAAGACTGCCTGGTGCGCTACGGTTGTCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTTTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACtt                                                                     <  2:120956/149‑1 (MQ=21)
                                                                            catcGGCCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCTGCATAAACTAAGCGCACttt                                                                    <  2:13745/149‑1 (MQ=255)
                                                                                    tGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAg                                                            <  2:1126587/149‑1 (MQ=255)
                                                                                      agaTTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTc                                                          <  1:76467/149‑1 (MQ=255)
                                                                                         tttATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCCACAGTCTAg                                                       >  2:1098182/1‑149 (MQ=255)
                                                                                          ttATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGc                                                      >  1:1041930/1‑149 (MQ=255)
                                                                                          ttATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGc                                                      <  1:1138068/149‑1 (MQ=255)
                                                                                                cagaCGCAAAACTTCCTGTTGCGCTACGCTTATCAGGCCTACGTCGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCCTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATg                                                <  2:149411/146‑1 (MQ=21)
                                                                                                cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATg                                                <  1:615030/149‑1 (MQ=255)
                                                                                                 aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACTAAGCGCACTTTGTCCACAGTCTAGCCCGAtcg                                               >  1:186479/1‑147 (MQ=255)
                                                                                                 aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTACGCGCACTTTGTCAACAGTCTAGCCCGATgg                                               >  1:567953/1‑149 (MQ=255)
                                                                                                 aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATgg                                               >  1:125879/1‑149 (MQ=255)
                                                                                                 aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATgg                                               >  2:442238/1‑149 (MQ=255)
                                                                                                     gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATc                                           <  1:1045465/149‑1 (MQ=255)
                                                                                                              ccTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCggg                                  <  2:688333/149‑1 (MQ=255)
                                                                                                               cTGATGGGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGc                                 <  1:138601/149‑1 (MQ=255)
                                                                                                               cTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGc                                 <  1:1107104/149‑1 (MQ=255)
                                                                                                                tGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTCATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGcc                                >  1:827523/1‑149 (MQ=255)
                                                                                                                   tGCGCTACGCTTATCAGGCCTACGTGGATCGAGCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAAGCGTTCACGCCGCATCCGGCATAAACTAAGAGCACTTTGTAAACAGTCTTACCCGATGTCATCACCATCGGGCctct                             >  1:620382/1‑149 (MQ=21)
                                                                                                                   tGCGCCACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCctct                             <  1:1120340/149‑1 (MQ=255)
                                                                                                                      gcTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCtttt                          >  1:1094948/1‑149 (MQ=255)
                                                                                                                       cTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCttttt                         >  2:826290/1‑149 (MQ=255)
                                                                                                                        tACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTttta                        >  1:355042/1‑149 (MQ=255)
                                                                                                                        tACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCCCTTttta                        >  2:218168/1‑149 (MQ=255)
                                                                                                                        tACGCTTATCAGGACTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAACGCACATTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTttta                        >  1:784672/1‑149 (MQ=255)
                                                                                                                         aCGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTtttat                       <  2:568806/149‑1 (MQ=255)
                                                                                                                           gCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTtttattt                     >  1:534203/1‑149 (MQ=255)
                                                                                                                               aTCAGGCCTACGTGGATCGATCAATTTAGGGAATTTACACAATTTTGTAGGCCGGCTAAGGCGTTCACGCCGCCTCCGGCATAAACTAAGCGCACTTTGTCAACAGGCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTActc                 >  1:617523/1‑149 (MQ=255)
                                                                                                                                 cAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACtctc               >  2:202989/1‑149 (MQ=255)
                                                                                                                                    gCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCGTCTTTTTATTTACTctgctg            >  1:332722/1‑149 (MQ=255)
                                                                                                                                     ccTACGTTGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTgc           <  1:866446/149‑1 (MQ=255)
                                                                                                                                     ccTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTgc           <  1:629496/149‑1 (MQ=255)
                                                                                                                                      cTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTgcg          <  2:557647/149‑1 (MQ=255)
                                                                                                                                          gTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCCCCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGa      <  2:945194/149‑1 (MQ=255)
                                                                                                                                          gTGGATCGATCAAGTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGa      <  2:713807/149‑1 (MQ=255)
                                                                                                                                              atcgatcAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACaaa  >  2:331091/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGGGGCACAATAAAAATAGCCCGATGGATTTATCATCGGGCTGAGATTTATGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCGATCAATTTATTGAATTTACACAATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACTAAGCGCACTTTGTCAACAGTCTAGCCCGATGGCATCACCATCGGGCCTCTTTTTATTTACTCTCCTGCGGCGACAAA  >  NC_000913/2291090‑2291380

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: