Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,857 A→T 6.5% intergenic (+24/+22) qseC → / ← ygiZ quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,171,8570AT 6.5% 149.1 / 2.9 62intergenic (+24/+22)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base T (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-01

AAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAAT  >  NC_000913/3171712‑3172004
                                                                                                                                                 |                                                                                                                                                   
aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                  >  2:423705/1‑149 (MQ=255)
 aaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                                                                 <  2:843069/149‑1 (MQ=255)
 aaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                                                                 >  1:559132/1‑149 (MQ=255)
  aCCGCTACCGGCAGCGGGCTTGGGCTATAGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGTAATGCGGGACAAGCTGGATTTGAGGAGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGTGCAGTTAAAACg                                                                                                                                                >  2:771258/1‑149 (MQ=255)
    cGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACgcc                                                                                                                                              >  2:926444/1‑147 (MQ=255)
      cTACCTGCAGCGGGCTTGGGATATCTATTGGCCAGCTAATCGCAAAATTGCATGGCATGAATGTTGAATTTGGGAATCCGGAACAAGGTGGATTTGCGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                                                            <  2:873228/149‑1 (MQ=255)
      cTACCGGCTGCTGGCTTGGGCTTTCGATTGTCCCGCTAATCGCCAAATTGCATGGCATGAATGTTCAATTTGGGAATGCGGAACAAGGTGTATTTGTGGCGAAGGTAACCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                                                            <  1:1020570/149‑1 (MQ=255)
           ggCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                                                                       <  1:21309/149‑1 (MQ=255)
             cAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGATGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                                                                                     <  1:1150291/149‑1 (MQ=255)
                  ggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                                                                                <  2:57797/149‑1 (MQ=255)
                  ggCTTGGGCTAGCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCGGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                                                                                <  1:62141/149‑1 (MQ=255)
                   gCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGTCATGAATGTTGAAATTGGGAATGCGGAACAAGGTTGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                                                               <  1:541408/149‑1 (MQ=255)
                     ttGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                                                                             >  1:604341/1‑149 (MQ=255)
                         gCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTTGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                                                                         >  1:204460/1‑149 (MQ=255)
                         gCTATCCTTTGTCCAGCGAATACCCAAATTGCATGGCATGAATTTTGAATTTGGGAATGCGCAACAAGGGCGAATTGACGCGAAGGTAAGCTGGAAACTTCTCTGTCTATGCGCGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                                                                         <  1:324208/149‑1 (MQ=255)
                               gATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                                                                                   >  1:647938/1‑149 (MQ=255)
                               gATTGTCCAGCGAATCGCCAAATTGCATGGCAGGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTATCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                                                                                   >  1:259132/1‑149 (MQ=255)
                                aTTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTATAAGCAAAata                                                                                                                  >  1:985728/1‑149 (MQ=255)
                                   gTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAACGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                                                                               <  2:654230/149‑1 (MQ=255)
                                        gCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGt                                                                                                          <  2:1010376/149‑1 (MQ=255)
                                           aaTCGCCAAATTGCATGTCTGGAAGGTTGAATATGGTAATGCGGAACAATGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGAAAAATATAAACCGTccc                                                                                                       <  2:1070671/149‑1 (MQ=255)
                                               gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTTACTTCTCTGTCTATGAGAGCCGTGAAAACGACTCTCATAGATTTTACGAATAGCAAAATATAAACCGTACCCaaa                                                                                                   >  2:713087/1‑149 (MQ=255)
                                                ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaa                                                                                                  <  1:948062/149‑1 (MQ=255)
                                                  aaaTTGCATGGCATGAATGTGGAATTTGGGAATGCGCAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                <  2:264766/149‑1 (MQ=255)
                                                     ttGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcc                                                                                             >  1:624358/1‑149 (MQ=255)
                                                         atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcctcca                                                                                         >  1:765181/1‑149 (MQ=255)
                                                              atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                                                                    >  2:900790/1‑149 (MQ=255)
                                                              atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                                                                    <  1:423705/149‑1 (MQ=255)
                                                              atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATACATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                                                                    >  1:924841/1‑149 (MQ=255)
                                                                     ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                                                                             >  1:1100054/1‑149 (MQ=255)
                                                                     ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATGTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAACCACCAACCACAAACCAg                                                                             >  2:514193/1‑149 (MQ=255)
                                                                      tGAATTTGGGAATGCGGAACAAGGTTGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTATCATAGATTTTACTAATTGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                                                                            <  1:1076641/149‑1 (MQ=255)
                                                                       gcaTTTGGGAATGCGGAACAAGGTGGATTTGCGGCGAAGGTAAGCTGGTAACTTATCTGTCTATGAGAGCCGTTAAAACGCCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAACGCCACCAACCACAAACCAGAc                                                                           >  1:887786/3‑149 (MQ=255)
                                                                          tttGGTCATGCGGAACAAGGTGTTTTTGAGTCGAAGGAAAGCTGGTAACTTCTCTGTCTAGGAGAGCCGTTAAAACGACACTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaa                                                                        <  2:765181/149‑1 (MQ=255)
                                                                                  tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                >  1:617232/1‑149 (MQ=255)
                                                                                   tcGGAAGAAGGTGGATTTGATGCGATGGTAAGCTGTTAACTTCTCTGTATATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                                                               <  1:1145526/148‑1 (MQ=255)
                                                                                   gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                                                               >  2:222336/1‑149 (MQ=255)
                                                                                   gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                                                               >  2:279130/1‑149 (MQ=255)
                                                                                    ctgAACATGGTGGATTTGAGGCTAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                                                              <  2:106334/147‑1 (MQ=255)
                                                                                    cGGACCAAGGTGGATTTGATGCGAAGGTAAGCTGGTAACTGCTCTGTCTCTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTTATAGCAAAATATAAACCGTCCACAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                                                              <  2:523090/149‑1 (MQ=255)
                                                                                       aaCATGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAATAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa                                                           <  2:553146/149‑1 (MQ=255)
                                                                                          aagggtgATTTTCGGCGAAGGTAAGCTGTTAACTTATCTGTCTATGAGAGCCGTTAAAACTACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc                                                        <  2:1075019/143‑1 (MQ=255)
                                                                                           aGGTGGATTTGAGGCGAAGGTAAGCTGGTATCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAATCCGTCCCCAAAAACGCCACCAACCACAAACCAGACAAATACCCAGAGAGAATACa                                                       >  1:1031435/1‑149 (MQ=255)
                                                                                             gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata                                                     >  2:85781/1‑149 (MQ=255)
                                                                                                  tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAccc                                                >  1:292923/1‑149 (MQ=255)
                                                                                                       ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAAc                                           <  1:411688/149‑1 (MQ=255)
                                                                                                        gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACt                                          >  1:480721/1‑149 (MQ=255)
                                                                                                          gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCg                                        >  1:336029/1‑149 (MQ=255)
                                                                                                           aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGa                                       >  2:23492/1‑149 (MQ=255)
                                                                                                           aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCAAAAACCAGACAAATACCGAGCGAGAATAAATATACCCAAAACTCGa                                       >  1:851139/1‑149 (MQ=255)
                                                                                                              gTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACAGAGCGAGAATACATATACCCAAAACTCGaaaa                                    >  2:1106201/1‑147 (MQ=255)
                                                                                                               tAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAAAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACAt                                   <  2:43176/149‑1 (MQ=255)
                                                                                                                 aGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACAttt                                 <  2:608985/149‑1 (MQ=255)
                                                                                                                             ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACGAATAGCAAAATATAAACCGTCCCCAAAAAATCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAATCTCGAACATTTCCCGCATAAaga                     >  1:904561/1‑149 (MQ=255)
                                                                                                                              tctGTCTATGAGAGCCGATAAAACGACTCTCATAGATATTACTAATAGCAAAATATAAACCGACCCCAAAAAATCCACCAACCACAAACCAGACAAATAACGAGCGAGAATACAAATACACAAAACACGAACATTACCCTCATAAagat                    >  2:28359/1‑148 (MQ=255)
                                                                                                                               aaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGt                   >  2:1097633/4‑149 (MQ=39)
                                                                                                                                tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAACCCACCAACCACAAACCAGACAAATACCCAGCGAGAATACATATAACCAAAACTCGAACATTTCCCGCAAAAAGAGtt                  >  2:41392/1‑149 (MQ=255)
                                                                                                                                    tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACACACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCt              >  1:27118/1‑149 (MQ=255)
                                                                                                                                     aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt             <  2:76735/149‑1 (MQ=255)
                                                                                                                                      tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTta            <  2:878778/149‑1 (MQ=255)
                                                                                                                                       tagagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAACGTCCCCAAAAAATCAACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAATCCTCGAACATTTCCCGCATAAAGAGTTTTCTtaa           <  1:713087/148‑1 (MQ=255)
                                                                                                                                       gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAAGACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa           >  2:571819/1‑149 (MQ=255)
                                                                                                                                                tAAATCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACATATACCGACCGAGAATACATATACCCAATACTCGAACAAATCCCGCATAAAGAGTTTTCTTAATATAAGaaa  <  2:69646/149‑2 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAAT  >  NC_000913/3171712‑3172004

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