Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,575,216 G→A 5.2% intergenic (‑14/+31) pqqL ← / ← yddB putative periplasmic M16 family zinc metalloendopeptidase/putative TonB‑dependent outer membrane receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,575,2160GA 5.2% 206.9 / 2.4 78intergenic (‑14/+31)pqqL/yddBputative periplasmic M16 family zinc metalloendopeptidase/putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base G (31/42);  new base A (2/2);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-01

CGCAAGCCATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCGG  >  NC_000913/1575069‑1575359
                                                                                                                                                   |                                                                                                                                               
cGCAAGCCATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGa                                                                                                                                                >  1:963161/1‑149 (MQ=255)
  cAAGCCATTGTCCAGTTGACCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTTACGTCAGTAAGAAACAGAGGTTGCTCATATTTATCGACATGCGAAAACCGGGCGAAt                                                                                                                                              >  2:751053/1‑149 (MQ=255)
     gCCATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAAGCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGCGCGAATTTa                                                                                                                                           <  1:931902/149‑1 (MQ=255)
      ccATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGAGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTAc                                                                                                                                          <  1:836837/149‑1 (MQ=255)
       cATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATTATTATCTCCATGCGAAAACCGGGCGAATTTAcc                                                                                                                                         >  1:824150/1‑149 (MQ=255)
       cATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTAcc                                                                                                                                         >  1:385979/1‑149 (MQ=255)
        aTTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAAGCAGCCGCCCCGGGAGCAACAGAGTTGCCACTAACGTCCGTAAGAAACAGCGGATTCTCATAAGTATCACCATGCGAAAACCGGGCGAAGTTAccc                                                                                                                                        >  2:997411/1‑149 (MQ=255)
             ccAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGTAAATTCGCCCGGttt                                                                                                                                   >  1:734646/1‑148 (MQ=255)
              cAGTTGCCCGTTAATTAACTTTTCATACTGCGGTAAGGCGGCGGCAATAAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGTAAATTCGCCCGGtttt                                                                                                                                  <  2:371406/149‑3 (MQ=255)
              cAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTaa                                                                                                                                  <  2:830732/149‑1 (MQ=255)
               aGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCACTCAGCCGCCCGGGGCGCAACAGAGTTGCCACTAACGTCAGTAAGAAAAAGAGGTTTCTCATAATTATCTCCTTGCGAAAACCTGGCGAATTTACCCGGTTAAg                                                                                                                                 >  2:226182/1‑149 (MQ=255)
                  tGCCCGGTAATTAACTTTTCATCCTGCGGTACGGCGGCGGCAATCAGCCGCACGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGTAAATTCGCCCGGttttcgca                                                                                                                              <  1:52622/149‑7 (MQ=255)
                   gCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGACACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTaaa                                                                                                                             >  1:644123/1‑149 (MQ=255)
                   gCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTATCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTaaa                                                                                                                             >  1:644126/1‑149 (MQ=255)
                   gCCCGCTAATTAACTTTTCATCCAGCGGTAAGGCGGCGGCAATCAGCCTCCCGTGTAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTaaa                                                                                                                             <  1:338366/149‑1 (MQ=255)
                      cGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATc                                                                                                                          <  2:329484/149‑1 (MQ=255)
                       ggTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATcc                                                                                                                         <  1:73532/149‑1 (MQ=255)
                          aaTTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGaa                                                                                                                      <  1:403148/149‑1 (MQ=255)
                           tttAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGACGCCCGTGGAGCAACAGAGTTGCCACTAACGTCAGTTAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAAc                                                                                                                     <  1:578276/148‑1 (MQ=255)
                             tAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTa                                                                                                                   >  1:350705/1‑149 (MQ=255)
                                cTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTa                                                                                                                >  1:460772/1‑149 (MQ=255)
                                     cATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAAtt                                                                                                           <  2:433871/149‑1 (MQ=255)
                                       tCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTc                                                                                                         <  1:176760/149‑1 (MQ=255)
                                        ccTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCTGGCGAATTTACCCGGTTAAGTAAAATCAGAACTATTAAAATTTCa                                                                                                        >  2:708146/1‑149 (MQ=255)
                                        ccTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCa                                                                                                        >  1:370598/1‑149 (MQ=255)
                                        ccTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCa                                                                                                        >  1:462330/1‑149 (MQ=255)
                                           gCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGc                                                                                                     <  1:119020/149‑1 (MQ=255)
                                              ttAAGGCGGCGGCAATCAGCCGCCAGGGGAGCAACAGATTTGCCACTAACGTGAGTAAGAAACAGAGTTTTCTCATAATTATCTCCATGCGAAAACCGGGCGATTTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGa                                                                                                  <  1:38176/148‑1 (MQ=255)
                                               tAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGAc                                                                                                 >  2:809096/1‑149 (MQ=255)
                                                     ggcggcAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAg                                                                                           <  1:897588/149‑1 (MQ=255)
                                                      gcggcAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGc                                                                                          <  2:999548/149‑1 (MQ=255)
                                                           aaTCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGaa                                                                                     >  1:424850/1‑149 (MQ=255)
                                                             tCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAg                                                                                   <  2:629342/149‑1 (MQ=255)
                                                                gccgccCGGGTAGCAACAGAGTTGCCACTAACGTCTGTAAGAAACAGAGGTTTCTCATATTTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTAc                                                                                >  2:324019/1‑149 (MQ=255)
                                                                 ccgccCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTATAATTTCATGCTGACATCAAGCCAGAAAGTACg                                                                               >  1:690596/1‑149 (MQ=255)
                                                                   gccCTGGGAGCAACATAGTTGCCTCTAACGTCAGTAATAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCAGAAATATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACTGc                                                                             >  2:1004919/1‑149 (MQ=255)
                                                                   gccCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGc                                                                             <  2:248787/149‑1 (MQ=255)
                                                                   gccCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGc                                                                             <  2:126790/149‑1 (MQ=255)
                                                                              aaCAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAg                                                                  <  2:118465/149‑1 (MQ=255)
                                                                              aaCAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTAGCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAATATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAg                                                                  >  1:692095/1‑149 (MQ=255)
                                                                                cAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGATCTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGt                                                                >  2:130748/1‑149 (MQ=255)
                                                                                cAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTAACCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGt                                                                >  1:722721/1‑149 (MQ=255)
                                                                                    gTTGCCACTAACGTCCGTAAGAAACAGACGTTTCTCATAATTATCTCCATGCGTAAACCGTGCGAATATACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGcc                                                            <  1:712728/149‑1 (MQ=255)
                                                                                     ttGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGccc                                                           >  2:713243/1‑148 (MQ=255)
                                                                                          aCTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGAACGGCCACTGGCGTAGGTCGCCACACCg                                                      >  1:94520/1‑149 (MQ=255)
                                                                                           cTAACGTCACAAAGAAACAGAGGTTTCTCATAATAATCACCAAGCAAAAACCGCGCTAATTTACCCGGTTAAGTAAAATCAAAACTAATAAAATATCATTATAAAATCAAACCAGAAAGTACGCACACTCGAGTAGATCGCcaccactt                                                     >  2:1002455/1‑144 (MQ=255)
                                                                                            tAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGgt                                                    <  2:46803/149‑1 (MQ=255)
                                                                                            tAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGgt                                                    <  1:809096/149‑1 (MQ=255)
                                                                                            aaaCGCCAGTAAGAAACAGAGGTTTCTCATAATTATCACCCCGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGgt                                                    <  1:1028264/148‑1 (MQ=255)
                                                                                              aCGTCAGTAAGAAACAGAGGTTTCTCAAAATTATCTCCATGCGAAAACCGGGCGAATATACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGgtgt                                                  >  2:845619/1‑149 (MQ=255)
                                                                                                gTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAACTACGGCCACTGTCGTAGGTCGCAACACCGgtgtaa                                                >  2:94379/1‑147 (MQ=255)
                                                                                                 gaaGTAAGAAACAGAGGTCTCTCATAATAATCTCCATGCTAAAACCGTGCGAATTAACCCGGATAAGTAAAATCCGAACTATTAAAATTTCCTGCTGACATCAAGCCCGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCg                                               <  1:836566/147‑1 (MQ=255)
                                                                                                    gTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTg                                            <  2:855953/149‑1 (MQ=255)
                                                                                                    gTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTAGTTg                                            >  1:850082/1‑149 (MQ=255)
                                                                                                     tAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAATATCCGAACTATTAAAATTTCATGCTGACATCATGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGt                                           <  1:461960/149‑1 (MQ=255)
                                                                                                      aaGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTa                                          >  1:88888/1‑149 (MQ=255)
                                                                                                        gAAACAGAGGTTTCTCATAATTATCTCCATGCGAATACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATc                                        <  1:708146/149‑1 (MQ=255)
                                                                                                        gAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATc                                        <  2:765702/149‑1 (MQ=255)
                                                                                                          aaCAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCaa                                      <  2:482441/149‑1 (MQ=255)
                                                                                                             agagGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCACCAg                                   >  1:821004/1‑149 (MQ=255)
                                                                                                               tgttttCTCATAATAATCTCCATGCGAAAACCGGCCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGTCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCg                                 <  2:187893/146‑1 (MQ=255)
                                                                                                                       cATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCt                         <  2:644123/149‑1 (MQ=255)
                                                                                                                       cATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCt                         <  2:644126/149‑1 (MQ=255)
                                                                                                                          aaTTATCTCCATTAGAAAACCGGGCGAATTTAACAGGTTAAGTCAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACTGCAACTGGAGTATGTCGACACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATc                      <  2:488299/149‑1 (MQ=255)
                                                                                                                          aaTTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATc                      <  2:462330/149‑1 (MQ=255)
                                                                                                                          aaTTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATc                      <  1:228946/149‑1 (MQ=255)
                                                                                                                              aTCTCCATGCGAAAACCGGGCGAATTTACCTGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACTGACACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGt                  <  1:863372/149‑1 (MQ=255)
                                                                                                                               tctcCATGCGAAAACCGTGCGAATTTACCCGGGTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTa                 <  2:633388/149‑1 (MQ=255)
                                                                                                                                    tttacttaaCCGGGTAAATTCGCCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTGTCCAGTACATTg            >  1:948210/8‑149 (MQ=255)
                                                                                                                                    aTGCGAAATCCGGGCGAATTTACCCGGTTAAGTAATATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTg            <  1:997411/149‑1 (MQ=255)
                                                                                                                                      gCGAATACGGGGCGAATTTTACTGTTTAAGTAAAATCCGACATATTAAAATTTCATGCTGACATCATGCCAGAAAGTACGGGCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTTTGCCGTACATTgag          <  1:496965/149‑1 (MQ=255)
                                                                                                                                      gCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGTCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTgag          <  2:139919/149‑1 (MQ=255)
                                                                                                                                        gAAAACCGGGCGAATTTACCCGGTTAAGTACACTCCGAACTATTAAATTTTCATGCTGACATCAATCCAGAAAGTACGGCAACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACATCGGTTTTGCTATCCAGATCATTGAGaa        >  1:739857/1‑149 (MQ=255)
                                                                                                                                        gAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGACCTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTAGCAACAGAGGTTTTGCTATCCAGTACATTGAGaa        >  1:567656/1‑149 (MQ=255)
                                                                                                                                        gAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGaa        >  1:743283/1‑149 (MQ=255)
                                                                                                                                           aaCCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAatat     <  2:690596/149‑1 (MQ=255)
                                                                                                                                             ccGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCg   <  1:182963/149‑1 (MQ=255)
                                                                                                                                              cGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCgg  <  2:75937/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CGCAAGCCATTGTCCAGTTGCCCGGTAATTAACTTTTCATCCTGCGGTAAGGCGGCGGCAATCAGCCGCCCGGGGAGCAACAGAGTTGCCACTAACGTCAGTAAGAAACAGAGGTTTCTCATAATTATCTCCATGCGAAAACCGGGCGAATTTACCCGGTTAAGTAAAATCCGAACTATTAAAATTTCATGCTGACATCAAGCCAGAAAGTACGGCCACTGGCGTAGGTCGCCACACCGGTGTTCGTTGTATCAACAGCGGTTTTGCTATCCAGTACATTGAGAATATCGG  >  NC_000913/1575069‑1575359

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: