Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 216,074 T→A 5.0% intergenic (+95/+105) nlpE → / ← yaeF lipoprotein involved with copper homeostasis and adhesion/putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913216,0740TA 5.0% 279.7 / 4.3 100intergenic (+95/+105)nlpE/yaeFlipoprotein involved with copper homeostasis and adhesion/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (53/42);  new base A (3/2);  total (56/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.37e-01

Only 100 of 102 total aligned reads displayed.

GATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCC  >  NC_000913/215926‑216222
                                                                                                                                                    |                                                                                                                                                    
gATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGttt                                                                                                                                                      <  1:1269927/149‑1 (MQ=255)
    cGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                                  <  1:261292/149‑1 (MQ=255)
     ggCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                 >  2:122302/1‑149 (MQ=255)
        aGGGAAATTTTACCCCCCCCAGGAGTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAACCAGGCCAGAAGGATATATTTCAACCTTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTAcgccgaa                                                                                                                                              >  1:90705/1‑147 (MQ=39)
          ggAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGCAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                            >  2:1228673/1‑149 (MQ=255)
          ggAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCACCATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACTCCGCATc                                                                                                                                            >  1:1021620/1‑149 (MQ=255)
          ggAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                            >  2:825945/1‑149 (MQ=255)
          ggAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                            >  1:322379/1‑149 (MQ=255)
              aTTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                        <  1:793323/149‑1 (MQ=255)
               ttttACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                       <  1:570508/149‑1 (MQ=255)
                 ttACCCCAACCATGATTGCAGTACTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAATGATATATATCAACATTTTCAATTTGCACGTTTTTTGTAGGCCAGATGCGGCGTAAACGCCTTATCCGGccta                                                                                                                                     <  1:161732/149‑4 (MQ=255)
                    ccccAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGt                                                                                                                                  >  2:98524/1‑149 (MQ=255)
                      ccATCCTGGATTGCAGTAGTTTGGGGCTGTAACCCTTCTTGAGTCAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGc                                                                                                                                <  2:150202/149‑1 (MQ=38)
                           cAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAAc                                                                                                                           <  1:443088/149‑1 (MQ=255)
                              gtttGCAGTTGTTTGGGGCAGTAACCCGTCTTGAGACAGTAACAAACGTAAAACAGGCCAGATGGATATATTTCATCATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGcc                                                                                                                        <  2:711801/147‑1 (MQ=38)
                              gATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTGGTGCTCAACGcc                                                                                                                        >  2:709721/1‑149 (MQ=255)
                                 tGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCCAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGa                                                                                                                     <  2:1204423/149‑1 (MQ=255)
                                  gCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGACACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGAt                                                                                                                    <  1:937490/149‑1 (MQ=255)
                                  gCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCTACATTTTGAATTTGCACTTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAGGTGCTCAACGACTGAt                                                                                                                    >  1:1383091/1‑149 (MQ=255)
                                  gCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGAt                                                                                                                    <  1:776670/149‑1 (MQ=255)
                                     gtagtTTGGGGCAGTAACCCGTCTTGAGACAGCAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGTTTCTCAACGCCTGATGCg                                                                                                                 >  1:396436/1‑149 (MQ=255)
                                        gtTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACg                                                                                                              >  2:548122/1‑149 (MQ=255)
                                         tttGGGGCAGTACCCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGc                                                                                                             >  1:1499144/1‑149 (MQ=255)
                                             gggCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGgc                                                                                                         >  2:1152537/1‑149 (MQ=255)
                                                    aaCCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCtt                                                                                                  <  1:122302/149‑1 (MQ=255)
                                                      cccGTCTTGAGACAGAAACAAACTCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTTTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACTCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                <  1:120198/149‑1 (MQ=255)
                                                      cccGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                <  1:357115/149‑1 (MQ=255)
                                                       ccGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTAAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGAGTTTACGCCGCATCAGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGAGTCTTAta                                                                                               >  1:47708/1‑148 (MQ=255)
                                                        cGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                              <  1:441396/149‑1 (MQ=255)
                                                         gTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAt                                                                                             >  2:20830/1‑149 (MQ=255)
                                                           cTTGAGACAGAAACAAACGCAAAACAGGCCAGAATGAAATATTTCAACATTTTGAATTTGCACGTTTTTTTTAGGCCGGATAAGTCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGc                                                                                           <  2:567601/149‑1 (MQ=255)
                                                           cTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGc                                                                                           <  1:1226522/149‑1 (MQ=255)
                                                           cTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGc                                                                                           <  1:847858/149‑1 (MQ=255)
                                                            ttGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGcc                                                                                          <  1:1156386/149‑1 (MQ=255)
                                                            ttGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGcc                                                                                          >  1:1506080/1‑149 (MQ=255)
                                                            ttGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGcc                                                                                          <  1:548122/149‑1 (MQ=255)
                                                              gagaAAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGAGACGCTGGCGCGTCTTATCATGCCTa                                                                                        >  2:719865/1‑149 (MQ=255)
                                                               agaCAGAAACACACGCAAAACACGCCAGAAGCATATATTTCACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATACCGCGTTTACGCCGCATCCGGCAATGGTGCGCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACa                                                                                      <  1:490192/149‑1 (MQ=25)
                                                                  cAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAc                                                                                    >  1:283459/1‑149 (MQ=255)
                                                                   aGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTATGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGCTGCTCAACGCCTGATCCGACGCTGGAGCGTCTTATCATGCCTACAAcc                                                                                   >  1:481234/1‑149 (MQ=255)
                                                                     aaacaaacGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAcccc                                                                                 >  1:987646/1‑149 (MQ=255)
                                                                      aacaaacCCAAAACAGGCCAGAACGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAccccc                                                                                <  2:355643/149‑1 (MQ=255)
                                                                        caaacGCAAAACAGGCCAGAATGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCt                                                                              <  2:416854/149‑1 (MQ=255)
                                                                        aaaacGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTGTTTTGGAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCt                                                                              <  2:1535429/148‑1 (MQ=255)
                                                                         aaacGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACTCCGCATCCGGCAATCGTCCTCAACGCCTGATGCGACGCTGGCGCGTCATATCATGCATACAACCCCCCTc                                                                             >  2:1086835/1‑149 (MQ=38)
                                                                           ccGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCAt                                                                           >  2:1100770/2‑149 (MQ=255)
                                                                           acGCAAAACAGGCCAGAAGGATATATATCAACATTTTGAATTTGCACGTTATTAGTAGGCCGGATAAGCCGTATACGCCGCATCCGGCAATGGTGATCAAAGCCTGATGCGACTCTGGCGCGTCTTATCATGCCTACAACCCCcatcat                                                                           >  1:1147500/1‑149 (MQ=25)
                                                                             gCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTCCAACCCCCCTCATAc                                                                         >  2:747257/1‑149 (MQ=255)
                                                                             gCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATAc                                                                         >  2:741382/1‑149 (MQ=255)
                                                                             gCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATAc                                                                         >  1:179458/1‑149 (MQ=255)
                                                                             gCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCATTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATAc                                                                         <  1:220552/149‑1 (MQ=255)
                                                                               aaaaCAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCAGACCt                                                                       >  2:1487906/1‑149 (MQ=255)
                                                                                aaaCAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATCCCtt                                                                      >  2:150536/1‑149 (MQ=255)
                                                                                  aCAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTaa                                                                    <  1:542695/149‑1 (MQ=255)
                                                                                   ctggCCAGAAGGTTTTATTTCATCATTTTGAATTTGCACGTTTTTTGTAGGCCTGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAg                                                                   <  2:936556/147‑1 (MQ=255)
                                                                                         aGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGAt                                                             >  1:415216/1‑149 (MQ=255)
                                                                                         aGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGAt                                                             >  1:900308/1‑149 (MQ=255)
                                                                                         aGAAGGATATATTTCAACATTNTGAATTTGCNNGTTTTTTGTAGNCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGAt                                                             >  1:421100/1‑149 (MQ=255)
                                                                                          gAAGGATATATTTCTACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCTATGGTTCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAGCCCCCTCATACCTTAAGGCGGATa                                                            <  1:962383/149‑1 (MQ=255)
                                                                                           aacgATATAATTCAACATTTTGAATTTGCACGTTTTATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATaa                                                           <  2:355251/146‑1 (MQ=255)
                                                                                             ggATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAgg                                                         <  1:566033/149‑1 (MQ=255)
                                                                                              gATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGc                                                        <  1:1345456/149‑1 (MQ=255)
                                                                                              gATATATTTCAACATTTAGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGc                                                        <  2:921916/149‑1 (MQ=255)
                                                                                               atatatTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCa                                                       >  2:467649/1‑149 (MQ=255)
                                                                                                 atatTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCAtt                                                     <  1:1192549/149‑1 (MQ=255)
                                                                                                   atTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACTCAGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGTCGCGTCTTATCATGCCTACACCCCACCTCATACCTTAAGTCGGATAAGGCATTTa                                                   >  1:328282/1‑149 (MQ=37)
                                                                                                   atTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTTCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACATTAAGGCGGATAAGGCATTTa                                                   >  1:1079365/1‑149 (MQ=255)
                                                                                                    tttCAACATTTTGAATTTGCACGTTTTTTTTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTAc                                                  <  2:1142794/149‑1 (MQ=255)
                                                                                                    tttCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTAc                                                  >  1:664772/1‑149 (MQ=255)
                                                                                                     ttCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATTAGGCATTTACt                                                 >  1:1164580/1‑149 (MQ=255)
                                                                                                     ttCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACt                                                 >  2:786284/1‑149 (MQ=255)
                                                                                                     ttCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACt                                                 >  2:555495/1‑149 (MQ=255)
                                                                                                      tCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACTACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACtt                                                >  1:982542/1‑149 (MQ=255)
                                                                                                      ccAACATTTGAATTTGCACGTTTTAGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGGGAGCTCAACGCCTGATCCGACGCGGGCGCGGCGTATCATGCCTACACCCCCCCTCATCCCTTAAGGCGGATAAGGCATTTACtgt                                               <  2:1046057/148‑3 (MQ=25)
                                                                                                       cAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACttt                                               <  1:20830/149‑1 (MQ=255)
                                                                                                         aCATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTAt                                             >  1:601790/1‑149 (MQ=255)
                                                                                                          cATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGAATTTACTTAAta                                            >  1:244179/1‑148 (MQ=255)
                                                                                                               tGAATTTGCACGTTTTTTGTAGGCCGGATAAGTCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCgg                                       <  2:1021620/149‑1 (MQ=255)
                                                                                                               tGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCgg                                       <  1:709721/149‑1 (MQ=255)
                                                                                                                tactttGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGTAATGGGGCTCAACGCCTGATGCGTCGCTGGCGCGTCTTATCATGCCTTCAACCGCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCggg                                      <  2:1447548/146‑1 (MQ=37)
                                                                                                                        aCGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACca                              >  1:52246/1‑149 (MQ=255)
                                                                                                                        aCGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACca                              >  2:329255/1‑149 (MQ=255)
                                                                                                                          gTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATCCCTTACGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACcaac                            <  2:1429680/149‑1 (MQ=255)
                                                                                                                           ttttttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCCTACCTTAAGGCGGATAAGGCATTAACTTTATCACCGGGTCAAACcaacc                           >  2:500671/1‑148 (MQ=255)
                                                                                                                           ttttttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAACCcaaca                           >  1:1012068/1‑149 (MQ=255)
                                                                                                                           ttttttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACcaaca                           >  2:1276296/1‑149 (MQ=255)
                                                                                                                            tttttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACcaacaa                          >  1:172900/1‑149 (MQ=255)
                                                                                                                            tttttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGCGCTCAACGCCTGAAGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACcaacaa                          <  2:662610/149‑1 (MQ=255)
                                                                                                                             ttttGTAGGCAGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCATGATCCGACGCTGGAGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACaaa                         >  1:386514/1‑149 (MQ=255)
                                                                                                                                  tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGAc                    >  1:651673/1‑149 (MQ=255)
                                                                                                                                   ttgCCGGATGCGGCGTAAACGCCTTATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAACCCACCCAAGCGAcc                   >  2:1142239/3‑149 (MQ=37)
                                                                                                                                   ttgCCGGATGCGGCGTAAACGCCTTATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTAAAACCAACAAAGCGccc                   >  2:509441/3‑149 (MQ=38)
                                                                                                                                       cGGATGCGGCGTAAACGCCTTATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGcc               >  1:1522100/1‑149 (MQ=38)
                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGcc               >  2:369485/1‑149 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCtt             >  2:491266/1‑149 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCtt             <  2:361870/149‑1 (MQ=255)
                                                                                                                                           tAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGa           <  1:188045/149‑1 (MQ=255)
                                                                                                                                           tAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGAGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGa           <  2:599554/149‑1 (MQ=255)
                                                                                                                                             aGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGata         >  1:522760/1‑149 (MQ=255)
                                                                                                                                                    tACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGCCCCGCCTTGATATAAATcc  >  1:293862/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCC  >  NC_000913/215926‑216222

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: