Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 566,354 A→G 8.9% pseudogene (23/279 nt) exoD ← pseudogene, DLP12 prophage; phage‑type exonuclease family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,3540AG 8.9% 252.1 / 13.0 91pseudogene (23/279 nt)exoDpseudogene, DLP12 prophage; phage‑type exonuclease family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (46/36);  new base G (6/2);  total (52/39)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCG  >  NC_000913/566207‑566501
                                                                                                                                                   |                                                                                                                                                   
cGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGCGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGAGTCACCCACATGCTGTACTGCACCTGGGCCATGTAACCCGATTTTATAGCCTCGAAACCACCGCGCCGGAAt                                                                                                                                                   >  2:3775/1‑150 (MQ=255)
 gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGGAtt                                                                                                                                                  >  2:486580/1‑150 (MQ=255)
      aTCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                                                                             <  1:1094609/150‑1 (MQ=255)
      aTCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                                                                             <  1:433450/150‑1 (MQ=255)
            aCATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATcc                                                                                                                                       <  2:411140/150‑1 (MQ=255)
               tAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                                                                    <  2:1009555/150‑1 (MQ=255)
               tAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                                                                                                                    >  2:498460/1‑150 (MQ=255)
                aaTGCATGCCTTAACGCTGCATGCGAGGTTCATTGTTTGCAACGTACCAGGCATCATTTCGCGTCACCCACATGCTGTACTTCACCTTGGCCATGTAAGCCGATTTTATTGCCTCGAAACGACCGAGCCGGAATTTCATGAAATCCCggg                                                                                                                                   <  1:857734/150‑1 (MQ=255)
                aaTGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGTCCATGTAAGCCGATTTATTTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                                                                                                                   <  2:592721/150‑1 (MQ=255)
                  tGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGGATTTCATGAAATTCCggctc                                                                                                                                 <  1:486580/150‑4 (MQ=255)
                  tGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAg                                                                                                                                 <  2:970610/150‑1 (MQ=255)
                   gCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAgg                                                                                                                                >  1:602764/1‑150 (MQ=255)
                   gCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAgg                                                                                                                                <  2:594338/150‑1 (MQ=255)
                    cAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGt                                                                                                                               >  2:270575/1‑150 (MQ=255)
                      ggCCTTCACGCTTCATGCGCGGGTCATAGTTGGCATAGTACCAGGCTTCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGGATTTCATGAAATTCCggctcggtg                                                                                                                             <  2:582930/150‑8 (MQ=255)
                          ttCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGt                                                                                                                         <  2:279113/150‑1 (MQ=255)
                          ttCACGCTTCATGCGCGGGGCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGt                                                                                                                         <  2:259718/150‑1 (MQ=255)
                            cACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCGCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTAt                                                                                                                       >  2:21679/1‑150 (MQ=255)
                                cTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGcc                                                                                                                   >  2:1211099/1‑150 (MQ=255)
                                   cATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCAGGa                                                                                                                >  1:1201355/1‑150 (MQ=255)
                                           ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                                                                        >  1:604159/1‑150 (MQ=255)
                                             tCATAGTTGGGAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCGGTACTGCACCTGGTCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCt                                                                                                      <  1:119203/150‑1 (MQ=255)
                                             tCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACCTTCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGTAATTTCATGAAATCCCGGGAGGTACTAGTATTGCCGGAAGCGTGTCCt                                                                                                      <  2:462964/150‑1 (MQ=255)
                                              cATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTg                                                                                                     <  2:882972/150‑1 (MQ=255)
                                              cATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTg                                                                                                     <  2:869583/150‑1 (MQ=255)
                                                  gTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATc                                                                                                 <  2:1081050/150‑1 (MQ=255)
                                                    tGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCg                                                                                               >  1:909423/1‑150 (MQ=255)
                                                      gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGTGCCATGTAAGCCGAGTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGa                                                                                             >  2:1017342/1‑150 (MQ=255)
                                                      gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGCAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGa                                                                                             >  2:804439/1‑150 (MQ=255)
                                                        aaaGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATg                                                                                           <  1:1284830/150‑1 (MQ=255)
                                                        aaaGTACCAGGCATCTTTTCGCGTCACCCACATGCTGGACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCAACGAGCCGGTATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATg                                                                                           <  2:524022/150‑1 (MQ=255)
                                                           gTACCAGGCAGCTTCTCGCGTCACCCCCATGCTGTACTGCACCTGGGCCATGTAAGCCGAGTTTATTGCCGCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATAGCCGGAAGCGGGGCCTGTATCCGGATGCag                                                                                        <  1:1050742/150‑1 (MQ=255)
                                                              ccAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAACCCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGt                                                                                     >  1:319142/1‑150 (MQ=255)
                                                               cAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTc                                                                                    <  2:596466/150‑1 (MQ=255)
                                                                aGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGtat                                                                                   >  2:642214/1‑148 (MQ=255)
                                                                aGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCt                                                                                   >  1:1027210/1‑150 (MQ=255)
                                                                aGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCt                                                                                   <  1:350123/150‑1 (MQ=255)
                                                                aGGCATATTTTCCCGACACCCACATTCTGTACTGCAACTGGGCCATGAAAGGCCATGTTATTGCCTCGAACCCCCCGCGCCGGAAATTCATCAAATCCCGGGACGTAAGAGAATTGCCGGAACCGTGGCCTGTATCCGGATGCAGAGtgt                                                                                   >  2:981838/1‑148 (MQ=255)
                                                                 ggCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCtt                                                                                  <  2:352593/150‑1 (MQ=255)
                                                                       ttttCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGACTCGAAACCACCGAGCCGCAATTTCATGAATTCCCGGGAGGTACGAGTATTGCCGGAAGCCTGGCCTGTATCCGGATGAAGAGACTTATcagg                                                                            >  1:176897/1‑147 (MQ=255)
                                                                        tttCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTg                                                                           <  1:228256/150‑1 (MQ=255)
                                                                         ttCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTgg                                                                          >  2:400776/1‑150 (MQ=255)
                                                                         ttCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTgg                                                                          <  1:349773/150‑1 (MQ=255)
                                                                         ttCGCGGCACCCCCCTGCTGTACTGCACCTGGGCAATGTAAGCCGATTTTATTGCCTCGCAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTTGCCTGTATCCGGATGCAGAGTCTTATCCGTgg                                                                          <  1:818640/150‑1 (MQ=255)
                                                                                cACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGa                                                                   >  2:7332/1‑150 (MQ=255)
                                                                                          cTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTa                                                         <  1:1172252/150‑1 (MQ=255)
                                                                                          cTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTa                                                         >  1:1264214/1‑150 (MQ=255)
                                                                                           tGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCAGTGGAAATCGAACGCGCATTAc                                                        >  2:398099/1‑150 (MQ=255)
                                                                                            gTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACt                                                       >  1:547088/1‑150 (MQ=255)
                                                                                              gcTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGCGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATCACtgg                                                     >  1:447926/2‑150 (MQ=255)
                                                                                              aCTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtgg                                                     >  1:826681/1‑150 (MQ=255)
                                                                                                tGCCCCTGGGCCATGTAAGCCGAGTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggtt                                                   <  2:107327/150‑1 (MQ=255)
                                                                                                  cACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggttgg                                                 <  1:14880/150‑1 (MQ=255)
                                                                                                    ccTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTggttggtt                                               >  2:1154490/1‑150 (MQ=255)
                                                                                                    ccTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTggttggtt                                               >  1:1199643/1‑150 (MQ=255)
                                                                                                    ccTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCCGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTggttggtt                                               >  1:33268/1‑150 (MQ=255)
                                                                                                     cTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGATCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTa                                              <  1:875749/150‑1 (MQ=255)
                                                                                                     cTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAAACGAACGCGCATTACTGGTTGGTTa                                              >  1:148639/1‑150 (MQ=255)
                                                                                                        ggCCATGTAAGCAGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCa                                           <  1:1215440/150‑1 (MQ=255)
                                                                                                         gCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTTCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCTTGGAAATCGAACGCGCATTACTGGTTGGTTTCCaa                                          >  1:832176/1‑150 (MQ=255)
                                                                                                         gCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCaa                                          >  2:80896/1‑150 (MQ=255)
                                                                                                           cATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACt                                        <  2:947051/150‑1 (MQ=255)
                                                                                                           cATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACt                                        >  2:17169/1‑150 (MQ=255)
                                                                                                               tAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTa                                    >  2:1309766/1‑150 (MQ=255)
                                                                                                                   ccGATTTTATTGCCTCGAAACCACCGCGCCGGACTTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGGTGGTTACCAACTTGTACCAg                                <  2:1281632/150‑1 (MQ=255)
                                                                                                                    cGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGa                               <  2:602764/150‑1 (MQ=255)
                                                                                                                     gATTTTATTGCCTCGATACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGTAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGaa                              <  1:932281/150‑1 (MQ=255)
                                                                                                                     gATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGaa                              >  1:705156/1‑150 (MQ=255)
                                                                                                                         ttATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACAAACTTGTACCAGAACATg                          >  1:666223/1‑150 (MQ=255)
                                                                                                                           atagCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTGATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCTACTTGTCCCAGAACTTGCg                        >  2:1213921/4‑150 (MQ=255)
                                                                                                                           aTTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCg                        >  2:110161/1‑150 (MQ=255)
                                                                                                                           aTTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCg                        <  1:21679/150‑1 (MQ=255)
                                                                                                                              gCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTAGCAGAACATGCGGGc                     >  2:1243346/1‑150 (MQ=255)
                                                                                                                               ccgcGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGcc                    <  1:238655/147‑1 (MQ=255)
                                                                                                                               ccTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCGTGTATCCGGATGCAGAGTCTTTGCCGTGGAATTCGAACGCGGATTACTGGTTGGTTACCATCTGGGCCCCGAACCTGCggggc                    >  2:191909/1‑148 (MQ=255)
                                                                                                                               ccTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGcc                    >  1:667354/1‑150 (MQ=255)
                                                                                                                                       ccaccGAGCCGGAATTTCATGAAATCCCGGGAGGGACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCgct            >  2:856516/1‑150 (MQ=255)
                                                                                                                                        caccGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCgctg           <  1:711990/150‑1 (MQ=255)
                                                                                                                                         accGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGAgctgg          >  1:504261/1‑150 (MQ=255)
                                                                                                                                            gAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCgctggctg       >  1:262163/1‑150 (MQ=255)
                                                                                                                                            ctcccGGGATTTCATGAAATTCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCgctggctg       >  2:1192373/4‑150 (MQ=255)
                                                                                                                                              cccGGGATTTCATGAAATTCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCgctggctgcc     >  2:1126260/2‑148 (MQ=255)
                                                                                                                                              cccGGGATTTCATGAAATTCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGAc     >  2:682293/2‑150 (MQ=255)
                                                                                                                                              cccGGGATTTCATGAAATTCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTAATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTAACAACTTGTACCAGAACATGCGGGCCAATCCGCTGGCTGAc     >  2:92589/2‑150 (MQ=255)
                                                                                                                                               ccGGGATTTCATGAAATTCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGAcg    >  1:320625/1‑150 (MQ=255)
                                                                                                                                               ccGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGAcg    >  1:516837/1‑150 (MQ=255)
                                                                                                                                               ccGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGGGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGACCATGCGGGCCACTGCGCTGACTGAcg    >  2:65205/1‑150 (MQ=255)
                                                                                                                                               ccGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCATGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTCCCACCTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGAcg    >  2:100571/1‑150 (MQ=255)
                                                                                                                                                 ggAATTTCTTGAAATCCCGGGAGGTACGAGTATTGCCGGTAGCGTGGCCTGTATCCGGATGCAGAGTGTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGAcgcg  <  2:1001454/150‑1 (MQ=255)
                                                                                                                                                 ggAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGAcgcg  >  2:360711/1‑150 (MQ=255)
                                                                                                                                                 ggAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGAcgcg  >  2:615803/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCG  >  NC_000913/566207‑566501

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: