Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,889,992 T→A 5.9% N181Y (AAT→TAT)  yeaY ← Slp family lipoprotein YeaY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,889,9920TA 5.9% 294.3 / 5.2 94N181Y (AAT→TAT) yeaYSlp family lipoprotein YeaY
Reads supporting (aligned to +/- strand):  ref base T (39/49);  new base A (2/4);  total (41/53)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCATCACCATAAATTTATAGGGCG  >  NC_000913/1889843‑1890138
                                                                                                                                                     |                                                                                                                                                  
tACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGAtt                                                                                                                                                    >  2:849049/1‑150 (MQ=255)
 aCCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCTTAAGAATGCAGAAACAGCGGACCAGCCTCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATta                                                                                                                                                   <  1:75295/150‑1 (MQ=255)
  cccAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATtat                                                                                                                                                  <  2:121299/150‑1 (MQ=255)
  cccAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTAAAACTGTTTGTACTCTCGCGGGGCCGGGATtat                                                                                                                                                  >  2:653630/1‑150 (MQ=255)
    cAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAc                                                                                                                                                >  2:294491/1‑150 (MQ=255)
    cAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAc                                                                                                                                                >  2:1418612/1‑150 (MQ=255)
     aaCGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAcc                                                                                                                                               <  1:184348/150‑1 (MQ=255)
     aaCGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAcc                                                                                                                                               >  2:1050041/1‑150 (MQ=255)
        gAAAAGGTTGCGAAGCGCGTCACTATTTATATTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGc                                                                                                                                            >  2:425738/1‑150 (MQ=255)
          aaaGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGccc                                                                                                                                          >  1:1068268/1‑150 (MQ=255)
            aGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCTGGTATAATCCCGGccccg                                                                                                                                        <  2:5595/150‑2 (MQ=255)
             ggTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCAt                                                                                                                                       >  1:552289/1‑150 (MQ=255)
              gTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCTGGTATAATCCCGGccccgcg                                                                                                                                      <  2:1494578/150‑4 (MQ=255)
                tGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCTGGTATAATCCCGGccccgcgag                                                                                                                                    <  1:957337/150‑6 (MQ=255)
                tGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATccg                                                                                                                                    >  1:1380046/1‑150 (MQ=255)
                  cGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCTGGTATAATCCCGGccccgcgagag                                                                                                                                  <  1:99407/150‑8 (MQ=255)
                  cGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATccgcc                                                                                                                                  <  2:99990/150‑1 (MQ=255)
                  cGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGACTATACCAGCCCCATccgcc                                                                                                                                  >  1:951860/1‑150 (MQ=255)
                  cGAAGCGCGTCACTATTTATTTTTATCTTTACAGTAAGAATGCAGAAACAGCGGACCAGAAGCTGTTTCTTTTTCTTTAAAAACAAGCAAATATTCAGTTACAACAGTTTGTACTCTAGCGCGGCCGGGATAATACCAGCCCCATcagcc                                                                                                                                  >  2:626312/1‑150 (MQ=255)
                   gAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATccgccg                                                                                                                                 >  2:55283/1‑150 (MQ=255)
                    aaGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGt                                                                                                                                <  2:825965/150‑1 (MQ=255)
                      gcgcGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAg                                                                                                                              <  1:1418612/150‑1 (MQ=255)
                      gcgcGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAg                                                                                                                              <  1:674126/150‑1 (MQ=255)
                               tatttattTTTATCTTTACCGTAAGAATGCAGAAACATCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgc                                                                                                                     <  2:241026/150‑1 (MQ=255)
                                  ttatttTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgccag                                                                                                                  <  2:1174274/150‑1 (MQ=255)
                                    atttTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCTCCGTAGCCATAGGgccagcc                                                                                                                <  2:552289/150‑1 (MQ=255)
                                     tttttATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgccagccc                                                                                                               >  1:353062/1‑149 (MQ=255)
                                         tATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGAc                                                                                                           <  2:1040671/150‑1 (MQ=255)
                                            cTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGc                                                                                                        <  2:1477096/150‑1 (MQ=255)
                                             tttACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGcc                                                                                                       <  1:1353903/150‑1 (MQ=255)
                                             tttACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGcc                                                                                                       >  1:1233796/1‑150 (MQ=255)
                                             tttACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGcc                                                                                                       <  1:1534398/150‑1 (MQ=255)
                                              ttACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCa                                                                                                      <  2:1318044/150‑1 (MQ=255)
                                              ttACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCa                                                                                                      <  2:1388497/150‑1 (MQ=255)
                                                aCCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATa                                                                                                    <  1:445904/150‑1 (MQ=255)
                                                   gTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATaaaa                                                                                                 >  2:1145075/1‑150 (MQ=255)
                                                    tAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAAc                                                                                                <  2:224514/150‑1 (MQ=255)
                                                    tAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAAc                                                                                                <  1:928505/150‑1 (MQ=255)
                                                         aTGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATgg                                                                                           <  1:618303/150‑1 (MQ=255)
                                                         aTGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATgg                                                                                           <  2:595008/150‑1 (MQ=255)
                                                          tGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGa                                                                                          >  1:1140806/1‑150 (MQ=255)
                                                           gCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATg                                                                                            >  2:756801/1‑147 (MQ=255)
                                                           gCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATg                                                                                            <  1:756801/147‑1 (MQ=255)
                                                           aaaGAAACAGCGGCTGGTCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGAt                                                                                         >  1:1179561/3‑150 (MQ=255)
                                                            cAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATc                                                                                        <  1:192191/150‑1 (MQ=255)
                                                            cAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATc                                                                                        >  1:637763/1‑150 (MQ=255)
                                                            cAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATc                                                                                        <  2:433790/150‑1 (MQ=255)
                                                              caaaCAGCGGACAAGCCGATGTTTCTTTTTCTTTAAAAACAAGAAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCaa                                                                                      >  2:1529033/2‑150 (MQ=255)
                                                                 aCAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCg                                                                                   <  2:960909/150‑1 (MQ=255)
                                                                  cAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTATAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCgg                                                                                  >  1:264105/1‑150 (MQ=255)
                                                                    gCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCt                                                                                <  1:520221/150‑1 (MQ=255)
                                                                    gCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCt                                                                                <  1:425738/150‑1 (MQ=255)
                                                                       gACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAg                                                                             >  2:1464986/1‑150 (MQ=255)
                                                                       gACCAGCCGCTGTTTCTTCTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAg                                                                             <  2:1431567/150‑1 (MQ=255)
                                                                        aCCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAgg                                                                            >  2:481116/1‑150 (MQ=255)
                                                                           aGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAggcgg                                                                         >  1:1369169/1‑150 (MQ=255)
                                                                             ccGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCa                                                                       <  1:502727/150‑1 (MQ=255)
                                                                              cGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCAt                                                                      <  1:481116/150‑1 (MQ=255)
                                                                               gCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATa                                                                     <  1:849049/150‑1 (MQ=255)
                                                                                   tttctttttcttTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCa                                                                 >  2:1522464/1‑150 (MQ=255)
                                                                                     tctttttctttAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCa                                                                       >  1:139442/1‑142 (MQ=255)
                                                                                     tctttttctttAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCa                                                                       <  2:139442/142‑1 (MQ=255)
                                                                                     tctttttctttAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCAcc                                                               <  1:128414/150‑1 (MQ=255)
                                                                                     tctttttctttAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCAcc                                                               <  1:1393325/150‑1 (MQ=255)
                                                                                     tcttttactttAAAAACAAGCAAATATTAAGTTAAAACTGTTTTTACTATCGCGGGGCCGGGATTATACAATCAACATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCAcc                                                               >  2:801857/1‑150 (MQ=255)
                                                                                       tttttctttAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAAt                                                                   <  2:1198758/144‑1 (MQ=255)
                                                                                       tttttctttAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAAt                                                                   >  1:1198758/1‑144 (MQ=255)
                                                                                               tAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTa                                                     >  1:105095/1‑150 (MQ=255)
                                                                                               tAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCAGCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTa                                                     <  2:1489512/150‑1 (MQ=255)
                                                                                                   aaCAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATg                                                 <  1:55283/150‑1 (MQ=255)
                                                                                                   aaCAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGCCGGCATAATCACCTGCTGGGTTAAATg                                                 >  2:1403621/1‑150 (MQ=255)
                                                                                                           attattCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGtt                                         <  2:637763/150‑1 (MQ=255)
                                                                                                              attCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGt                                      >  2:15841/1‑150 (MQ=255)
                                                                                                                 cAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAAc                                   >  1:1128722/1‑150 (MQ=255)
                                                                                                                  aGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAAcc                                  <  2:1233796/150‑1 (MQ=255)
                                                                                                                  aGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAAcc                                  >  2:1173894/1‑150 (MQ=255)
                                                                                                                   gTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAAccc                                 >  2:796660/1‑150 (MQ=255)
                                                                                                                   gTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAAccc                                 >  2:725074/1‑150 (MQ=255)
                                                                                                                   gTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAAccc                                 <  1:1173894/150‑1 (MQ=255)
                                                                                                                    ttACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCg                                <  1:1529033/150‑1 (MQ=255)
                                                                                                                     tACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAAcc                                  <  1:1015890/147‑1 (MQ=255)
                                                                                                                     tACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGt                               >  1:865240/1‑150 (MQ=255)
                                                                                                                     tACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGAATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAAcc                                  >  2:1015890/1‑147 (MQ=255)
                                                                                                                       cAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAACCCATGGATCAGTCGGCTGAGGCGGCATAATCACCTGCGGGGTTAAATGCCACGGTTTGTGGCCCGtgg                             >  2:1260635/1‑148 (MQ=255)
                                                                                                                       cAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTa                             >  2:944424/1‑150 (MQ=255)
                                                                                                                         aCTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACt                           >  1:984095/1‑150 (MQ=255)
                                                                                                                           tGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTg                         <  1:1453243/150‑1 (MQ=255)
                                                                                                                              ttGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCAt                      >  1:1445018/1‑150 (MQ=255)
                                                                                                                              ttGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCAt                      <  1:1042203/150‑1 (MQ=255)
                                                                                                                               tGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCATc                     <  2:1380046/150‑1 (MQ=255)
                                                                                                                               tGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCATc                     <  2:1209786/150‑1 (MQ=255)
                                                                                                                                gTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCATCa                    <  2:984095/150‑1 (MQ=255)
                                                                                                                                       gatggggCTGGTATAATCCCGGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCCCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCATCACCATaaa             >  2:728226/4‑150 (MQ=255)
                                                                                                                                       cgcgGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCATCACCATaaa             >  2:473785/1‑150 (MQ=255)
                                                                                                                                             gCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCATCACCATAAATTtata       <  1:15841/150‑1 (MQ=255)
                                                                                                                                              ccGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCATCACCATAAATTTATAg      <  1:1070213/150‑1 (MQ=255)
                                                                                                                                               cGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCATCACCATAAATTTATAgg     <  1:944424/150‑1 (MQ=255)
                                                                                                                                                  gATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCATCACCATAAATTTATAGGGCg  <  1:725074/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGCATAATCACCTGCTGGGTTAAATGCCAACGTTTGTAACCCGTTACTTGCATCACCATAAATTTATAGGGCG  >  NC_000913/1889843‑1890138

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