Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,481 T→G 5.4% intergenic (+19/‑334) gtrS → / → tfaS serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage/pseudogene, CPS‑53 (KpLE1) prophage; tail fiber assembly protein fragment;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,470,4810TG 5.4% 512.8 / 14.4 178intergenic (+19/‑334)gtrS/tfaSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage/pseudogene, CPS‑53 (KpLE1) prophage; tail fiber assembly protein fragment;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (77/91);  new base G (4/6);  total (81/97)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 180 total aligned reads displayed.

ATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATG  >  NC_000913/2470333‑2470630
                                                                                                                                                     |                                                                                                                                                     
atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAAAATCCATGTAAACCACTAAGTTTAAATTGCCCGCAAAta                                                                                                                                                      >  2:1437927/1‑149 (MQ=255)
    agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg                                                                                                                                                  >  2:395951/1‑150 (MQ=255)
     gagaTTTAAGCAAGCATTTATTTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg                                                                                                                                                  <  2:929551/150‑1 (MQ=255)
     gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc                                                                                                                                                 >  2:883556/1‑150 (MQ=255)
            aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaat                                                                                                                                          <  2:1163990/150‑3 (MQ=255)
               caatcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc                                                                                                                                       <  1:629499/150‑1 (MQ=255)
               caatcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc                                                                                                                                       <  2:1082166/150‑1 (MQ=255)
               caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAGTGCGCGCCATGATTGGCGCGcaattta                                                                                                                                       <  2:1601159/150‑6 (MQ=255)
                aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt                                                                                                                                      >  2:1000570/1‑150 (MQ=255)
                aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAAGAAATCAAGACAAAATCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAATCATGGCGCGCACaatca                                                                                                                                      >  2:1331924/1‑147 (MQ=255)
                 agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaatttaaa                                                                                                                                     <  1:371473/150‑8 (MQ=255)
                 agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCAAATCATGGCGCGCACAAGCta                                                                                                                                     >  1:303085/1‑150 (MQ=255)
                    atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAATAGCGATGTAAAGCAATAAGTTTAAATTGCGCGCCAATCATGGCGCGCAAAAGCTataa                                                                                                                                  >  1:903136/1‑150 (MQ=255)
                       tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAc                                                                                                                               <  1:1281884/150‑1 (MQ=255)
                        atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAcc                                                                                                                              >  2:356564/1‑150 (MQ=255)
                          tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCa                                                                                                                             >  2:1161936/1‑150 (MQ=255)
                            tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAc                                                                                                                           >  1:393007/1‑150 (MQ=255)
                             tgAATGTATTTTTCGTTCGAACTTTGCTTAAACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAcc                                                                                                                          <  1:863854/149‑1 (MQ=255)
                             ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAcc                                                                                                                          >  2:1678366/1‑150 (MQ=255)
                               aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTa                                                                                                                        >  2:686472/1‑150 (MQ=255)
                                  gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAAAATAGCCAAGTAAAGCACTAAGTTTAAATTTCTCGCCAATCATGGCTCGCACAATATATAATAACAACCTAAtt                                                                                                                     >  1:872149/1‑150 (MQ=255)
                                   tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAttt                                                                                                                    >  2:41211/1‑150 (MQ=255)
                                       tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                                                                <  1:139336/150‑1 (MQ=255)
                                       tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                                                                <  1:1002922/150‑1 (MQ=255)
                                       tttCGTTCGACCTTTGCTTATAAATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAACCACTAAGTTAAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                                                                <  1:1287620/150‑1 (MQ=255)
                                          cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcctcc                                                                                                             <  2:1651923/150‑1 (MQ=255)
                                             tcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                                                                          <  1:950091/150‑1 (MQ=255)
                                                  cTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTagag                                                                                                     <  2:1710900/150‑1 (MQ=255)
                                                  atttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTagag                                                                                                     <  1:297145/149‑1 (MQ=255)
                                                      gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCTCAAATCATGGCGCGCACAAGCTATAATACCAAACTAATTTCTCCTCCTCTTAGAGTGAc                                                                                                 >  1:248451/1‑150 (MQ=255)
                                                      gCTTATACATAACAGTCGTTTTTTTATTTATAAAATAATTCAAAGAAATCAAGACAATTGCGATGTAAAGCAGTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCGATAATACCAACATAATTTCTCCGCCTCTTAGAGTGACt                                                                                                >  2:151280/1‑150 (MQ=255)
                                                       cTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTAAAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACt                                                                                                <  2:863286/150‑1 (MQ=255)
                                                         tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAATTTTAAATTGCCCGCCAATCAAGGCGCGCACAAGCTAAAATACCAACATAATTTCTACTCCTCTTAGAGTGACtat                                                                                              >  1:1097699/1‑150 (MQ=255)
                                                         tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCAAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCACTTAGAGTGACtat                                                                                              >  1:834104/1‑150 (MQ=255)
                                                             aaTAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtct                                                                                          <  1:52382/149‑1 (MQ=255)
                                                                aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                                                                       <  2:578091/150‑1 (MQ=255)
                                                                    ttCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATTTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                                                                                   <  1:1275843/149‑1 (MQ=255)
                                                                     tCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAg                                                                                  <  2:1378738/150‑1 (MQ=255)
                                                                        tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                                                               <  1:1732212/150‑1 (MQ=255)
                                                                        tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                                                                               <  2:1259820/150‑1 (MQ=255)
                                                                            tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                                                                           <  2:236268/150‑1 (MQ=255)
                                                                              tAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                                                                         <  1:1127823/150‑1 (MQ=255)
                                                                              aaaTTTAAAAAATAATTCAAATAAATAAACACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGt                                                                         <  1:1523043/149‑1 (MQ=255)
                                                                               aaTTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTa                                                                        <  2:112324/150‑1 (MQ=255)
                                                                                 tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGAACTAAGTTTAAATTGCGCGACATTCATGGCGCGCAAAAGCTATTATACCAACATTATTTCTCCTAGTCTTAGAGTGACTATATATCCTGATAGAATTGAGGTAtt                                                                      >  1:1147796/1‑150 (MQ=255)
                                                                                  ttATAAAATAATTCAAAGAAATCAAGACAATAGCTATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTg                                                                     <  1:1444661/150‑1 (MQ=255)
                                                                                  ttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCAAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTg                                                                     <  2:303085/150‑1 (MQ=255)
                                                                                   tataAAATAATTCAAATAAATCAAGACAATAGCGATTTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                                                                    <  1:1331924/150‑1 (MQ=255)
                                                                                   tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                                                                    >  2:978505/1‑150 (MQ=255)
                                                                                   tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                                                                    >  1:156502/1‑150 (MQ=255)
                                                                                   tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                                                                    >  1:732364/1‑150 (MQ=255)
                                                                                     taAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACACGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACt                                                                  <  2:1681111/150‑1 (MQ=255)
                                                                                       aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt                                                                <  2:240306/150‑1 (MQ=255)
                                                                                       aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt                                                                <  1:1095127/150‑1 (MQ=255)
                                                                                              tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc                                                         >  2:1351151/1‑150 (MQ=255)
                                                                                              tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc                                                         <  2:1282671/150‑1 (MQ=255)
                                                                                               cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc                                                        >  1:1243677/1‑150 (MQ=255)
                                                                                               cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc                                                        >  2:677032/1‑150 (MQ=255)
                                                                                                aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCACCCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATAGCGGTATTGACTATCAAATGcca                                                       >  2:1699200/1‑149 (MQ=255)
                                                                                                aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGccc                                                       >  1:1633197/1‑150 (MQ=255)
                                                                                                aaaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcca                                                       >  1:1537954/1‑149 (MQ=255)
                                                                                                 aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt                                                      >  2:525286/1‑150 (MQ=255)
                                                                                                 aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt                                                      <  2:612556/150‑1 (MQ=255)
                                                                                                 aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt                                                      >  1:1061129/1‑150 (MQ=255)
                                                                                                 aaGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCt                                                      <  1:1231789/150‑1 (MQ=255)
                                                                                                     aaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAtt                                                  <  2:441835/150‑1 (MQ=255)
                                                                                                      aTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTc                                                 <  1:455366/150‑1 (MQ=255)
                                                                                                        cAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATAtct                                                                                          >  1:401371/1‑107 (MQ=255)
                                                                                                         atgACAATAGCGATGTAAAGCACTAAGTTTTTATTGCGCGCCAATCATGGCGCGCACTATCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATTCCCTGATTCgtt                                              <  1:129634/148‑1 (MQ=255)
                                                                                                         aaGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgtt                                              <  1:1138241/150‑1 (MQ=255)
                                                                                                          aGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttg                                             <  1:1000410/150‑1 (MQ=255)
                                                                                                            aCAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttctt                                           >  2:1645769/1‑150 (MQ=255)
                                                                                                            aCAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCACTGATTCgttgtt                                           >  2:726156/1‑150 (MQ=255)
                                                                                                                 aGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTg                                      <  2:623957/150‑1 (MQ=255)
                                                                                                                  tcGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGt                                     <  1:477463/149‑1 (MQ=255)
                                                                                                                     aTGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAat                                  <  2:104278/150‑1 (MQ=255)
                                                                                                                     aTGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCCCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAat                                  <  2:1468247/150‑1 (MQ=255)
                                                                                                                        tAAATCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATAtc                               >  1:1624964/1‑150 (MQ=255)
                                                                                                                        tAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTTTATCTCCTGATAGAATTGCGGTATTGACTATCAAAGGCCCTGATTCGTTGTTTATTGTAatatg                               >  2:823490/1‑149 (MQ=255)
                                                                                                                           aGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcc                            <  2:70353/150‑1 (MQ=255)
                                                                                                                           aGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcc                            <  1:1026880/150‑1 (MQ=255)
                                                                                                                            gCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcct                           <  1:848383/150‑1 (MQ=255)
                                                                                                                             cACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcctc                          <  2:779811/150‑1 (MQ=255)
                                                                                                                             cACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcctc                          >  1:1435958/1‑150 (MQ=255)
                                                                                                                              aCTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCctct                         <  1:717321/150‑1 (MQ=255)
                                                                                                                               cTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcctcaa                        >  2:489664/1‑148 (MQ=255)
                                                                                                                                 aaGTTTAAATTGCGCTCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTTACTATATCTCCTGTTAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATc                      <  1:130572/150‑1 (MQ=255)
                                                                                                                                   gTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTg                    >  1:585585/1‑150 (MQ=255)
                                                                                                                                     ttAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCa                  <  2:710560/150‑1 (MQ=255)
                                                                                                                                      gcttgtGCGCGCCATGATTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAg                 >  1:166118/6‑150 (MQ=255)
                                                                                                                                        aaTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGAc               <  2:1720366/150‑1 (MQ=255)
                                                                                                                                         aTTGCGCGCCAATCATGGCGAGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACg              <  2:872149/150‑1 (MQ=255)
                                                                                                                                          gtGCGCGCCATGATTGGCGCGCACAATCTATAATACAAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGa             >  1:347698/2‑150 (MQ=255)
                                                                                                                                           tGCGCGCCATGATTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTg                                                                                                   >  2:873161/1‑63 (MQ=37)
                                                                                                                                           tGCGCGCCATGATTGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGAt            >  1:1368452/1‑150 (MQ=255)
                                                                                                                                              gcgcCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAAc         <  2:329207/150‑1 (MQ=255)
                                                                                                                                               cgcCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGCCTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACt        >  1:613716/1‑150 (MQ=255)
                                                                                                                                                   aaTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTaaa    <  1:883556/150‑1 (MQ=255)
                                                                                                                                                   aaTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTaaa    <  1:1043089/150‑1 (MQ=255)
                                                                                                                                                   aaTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTaaa    <  2:1561363/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATG  >  NC_000913/2470333‑2470630

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: