Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,474,878 T→C 5.2% intergenic (+22/‑106) yfdQ → / → yfdR CPS‑53 (KpLE1) prophage; uncharacterized protein/CPS‑53 (KpLE1) prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,474,8780TC 5.2% 295.2 / 5.3 97intergenic (+22/‑106)yfdQ/yfdRCPS‑53 (KpLE1) prophage; uncharacterized protein/CPS‑53 (KpLE1) prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base T (37/54);  new base C (2/3);  total (40/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.75e-01

TACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATG  >  NC_000913/2474729‑2475026
                                                                                                                                                     |                                                                                                                                                      
tACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCt                                                                                                                                                        <  2:1041257/150‑1 (MQ=255)
 aCTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTg                                                                                                                                                       <  2:496788/150‑1 (MQ=255)
 aCTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAAAGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTg                                                                                                                                                       <  1:663356/150‑1 (MQ=255)
    ggTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑gg                                                                                                                                                  >  2:349934/1‑150 (MQ=255)
       tCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑ggggg                                                                                                                                               <  1:753546/150‑1 (MQ=255)
          gcgcATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCACCTGC‑‑GGGGGCAt                                                                                                                                            >  2:1282461/1‑150 (MQ=255)
           cgcATTATTCAGCTGGAAGCGGTGCAGGAAGCTATGGCTAACGAATTTCGTGATCTGCTTGTTGCGAAATTCAAAGACAGCAAAGTAGCAACCTGTATTGGTAGTTTAACCTCCTGATTTCCTTACTGCAAATGCCCCTTC‑‑GGGGGCAtt                                                                                                                                           >  1:794796/1‑150 (MQ=255)
           cgcATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCAtt                                                                                                                                           <  2:722041/150‑1 (MQ=255)
            gcATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGC‑‑AGGGGCAttt                                                                                                                                          <  1:1070864/150‑1 (MQ=255)
            gcATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGC‑‑AGGGGCAttt                                                                                                                                          <  1:215681/150‑1 (MQ=255)
             cATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTa                                                                                                                                         >  2:175293/1‑150 (MQ=255)
              attattCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTAt                                                                                                                                        <  1:496773/150‑1 (MQ=255)
              attaatCCGCTGGCAGCGGTGCAGGAAGATAAGGCTAACGAATTTCCTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAAAGCCCCTGC‑‑GGGGGCATTTAt                                                                                                                                        <  1:621296/150‑1 (MQ=255)
               ttattCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCAGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATg                                                                                                                                       <  1:385915/150‑1 (MQ=255)
               ttaatCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAAATTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTCGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATg                                                                                                                                       <  1:835502/150‑1 (MQ=255)
                 attCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGACTTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGa                                                                                                                                     <  1:333954/150‑1 (MQ=255)
                  ttCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGC‑‑AGGGGCAtttgcagta                                                                                                                                    <  1:1265925/150‑7 (MQ=255)
                     aGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACg                                                                                                                                 <  2:1179375/150‑1 (MQ=255)
                     aGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACg                                                                                                                                 <  1:1139580/150‑1 (MQ=255)
                      gCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGt                                                                                                                                >  1:77157/1‑150 (MQ=255)
                      gCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGAAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGt                                                                                                                                >  2:438096/1‑150 (MQ=255)
                       cTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTTAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGTCATTTATGGAAACGTa                                                                                                                               <  1:832200/150‑1 (MQ=255)
                       cTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAACGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTa                                                                                                                               >  1:806682/1‑150 (MQ=255)
                       cTGGAAGCGGTGCAGGAAGATATGGCTAACGAATATCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTa                                                                                                                               <  1:1170458/150‑1 (MQ=255)
                          gAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAAtt                                                                                                                            <  1:478809/150‑1 (MQ=255)
                             gCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAAc                                                                                                                         <  1:748313/150‑1 (MQ=255)
                              cgctGCAGGAAGATATGGATAACGATTTTCGTTCTCTGGTTTTTGAGATATTCAAAGACAGCAAAGAAGAAACCTTTCTTGGTACTTTCTCCGCCTGATTTCATTACTTCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACt                                                                                                                        <  1:209448/147‑1 (MQ=255)
                                gTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCa                                                                                                                      >  1:83557/1‑150 (MQ=255)
                                  gCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGATACCTTTTTTGGTACTTTCTCCGCCTGATTTCATTACTGCAAATCCCCCTGC‑‑GGGGGCATTTATGGAAACGGAATTAACTCaat                                                                                                                    >  1:246761/1‑150 (MQ=255)
                                   cAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCaata                                                                                                                   <  1:38516/150‑1 (MQ=255)
                                        aGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                                                                              <  1:964453/150‑1 (MQ=255)
                                         tatatGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCAcc                                                                                                             <  2:1240847/149‑1 (MQ=255)
                                           aatGGCTAACTAATTTAGTTATCTTCTTGTTGAGAAATTCAAAGACATCAAAGTATAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCACCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCgg                                                                                                           <  2:209152/149‑1 (MQ=255)
                                              ggCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATg                                                                                                        >  1:768531/1‑150 (MQ=255)
                                               gCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATgg                                                                                                       <  2:360973/150‑1 (MQ=255)
                                                 tAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTg                                                                                                     >  1:1095535/1‑150 (MQ=255)
                                                    cGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAgg                                                                                                  <  2:9176/150‑1 (MQ=255)
                                                    cGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAgg                                                                                                  <  1:1066172/150‑1 (MQ=255)
                                                     gAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATTGTGAg                                                                                                   <  2:307489/150‑1 (MQ=255)
                                                        tttCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTAACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCtt                                                                                              >  1:51027/1‑150 (MQ=255)
                                                          tCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTcc                                                                                            >  1:697750/1‑150 (MQ=255)
                                                           cGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCt                                                                                           <  2:1058936/150‑1 (MQ=255)
                                                             tGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCAACGGATGGTGAGGGCTTCCttt                                                                                         >  2:152086/1‑150 (MQ=255)
                                                              gATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCtttt                                                                                        >  2:1216717/1‑150 (MQ=255)
                                                                 cTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTAcc                                                                                     <  2:1198676/150‑1 (MQ=255)
                                                                      tgttgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATAAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAAc                                                                                <  1:245844/150‑1 (MQ=255)
                                                                       gttgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACt                                                                               >  1:1094484/1‑150 (MQ=255)
                                                                        ttgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTc                                                                              >  1:1259010/1‑150 (MQ=255)
                                                                         tgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCa                                                                             >  1:316933/1‑150 (MQ=255)
                                                                           agaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAgc                                                                           <  1:771219/150‑1 (MQ=255)
                                                                            gaAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAgcg                                                                          <  1:39933/150‑1 (MQ=255)
                                                                               aTTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCTAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCgg                                                                       <  1:727857/150‑1 (MQ=255)
                                                                                ttCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATGTATGGAAACGTAATTACCTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGt                                                                      >  1:968046/1‑150 (MQ=255)
                                                                                  cAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGc                                                                    <  1:313408/150‑1 (MQ=255)
                                                                                    aaGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAg                                                                  >  2:112427/1‑150 (MQ=255)
                                                                                          gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATTATCACCGGATGGTGAGGGCTTCCTTTTACCAAAACTCAGCGCGGTGCAGCGCata                                                            >  1:859302/1‑150 (MQ=255)
                                                                                          gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCata                                                            >  1:747831/1‑150 (MQ=255)
                                                                                              aGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACg                                                        <  2:354001/150‑1 (MQ=255)
                                                                                                tAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTg                                                      <  2:77157/150‑1 (MQ=255)
                                                                                                  gAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGga                                                    >  2:989702/1‑150 (MQ=255)
                                                                                                   aaaCCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGgag                                                   <  2:768531/150‑1 (MQ=255)
                                                                                                   aaaCCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGgag                                                   <  2:707376/150‑1 (MQ=255)
                                                                                                   aaaCCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCAAAACTCAGCGCGGTGCAGCGCATATACGTGgag                                                   >  2:275134/1‑150 (MQ=255)
                                                                                                      ccTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGTAGCGCATATACGTGGAGAAc                                                >  2:696074/1‑150 (MQ=255)
                                                                                                      acTTTATTGGTCCATACACCGCCTGAGTTCATTACTGCAACTGCCCCTGC‑‑GGGGGCATTTATGGAAACGTTATTAACTCAATAATCACCGGTTGGTGAGGGCATCCTTTTACCCAAACTCAGCGCGGTGCAGAGCATATACGTGGAGAAc                                                <  2:499853/149‑1 (MQ=255)
                                                                                                        tttATTGGTACTATCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACaa                                              >  1:1057747/1‑150 (MQ=255)
                                                                                                           ataggTACTTTCACCGCCTGAATTCAATACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGATGGCTTCCTTTTACCCAAACTCAGCGCGCTGCAGCGCATATACGTGGATAACAAAAt                                           <  1:503263/147‑1 (MQ=255)
                                                                                                            ttGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATg                                          <  1:1273796/150‑1 (MQ=255)
                                                                                                              ggTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTc                                        >  2:452411/1‑150 (MQ=255)
                                                                                                                  cTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCAttt                                    >  1:316704/1‑150 (MQ=255)
                                                                                                                    ttCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATttat                                  <  2:1317947/150‑1 (MQ=255)
                                                                                                                      cACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATttatta                                >  2:1221052/1‑150 (MQ=255)
                                                                                                                         cGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTaaaa                             >  2:669313/1‑150 (MQ=255)
                                                                                                                         cGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTaaaa                             <  2:303140/150‑1 (MQ=255)
                                                                                                                          gCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAAc                            <  1:1083618/150‑1 (MQ=255)
                                                                                                                               aTTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACttttt                       <  1:1257186/150‑1 (MQ=255)
                                                                                                                               aTTTCATTACTGCAAATCCCCCGGC‑‑GGGGGCATTTATGGAAACGTAAATAACTCAAATATCACCGGATGGTGAGGGGGTCATTTTACACAAACTCAGCGCGGTGCAGCGAATATTCGTGGAGAACAAAATCTCTGTTATTTACACtattt                       >  2:140781/1‑150 (MQ=255)
                                                                                                                                tttCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTc                      <  2:83557/150‑1 (MQ=255)
                                                                                                                                 ttCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTcc                     <  2:1094484/150‑1 (MQ=255)
                                                                                                                                 ttCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTcc                     <  1:346080/150‑1 (MQ=255)
                                                                                                                                  tCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCg                    <  1:382598/150‑1 (MQ=255)
                                                                                                                                   cATAACTCCAAATCCCCCTGC‑‑GGGGGCATTTATGGAAACGTACAAAACTCAATAATCACCGGACGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCgg                   <  1:1294355/150‑1 (MQ=255)
                                                                                                                                    aTTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCggg                  >  2:67907/1‑150 (MQ=255)
                                                                                                                                     ttACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGa                 <  1:586550/150‑1 (MQ=255)
                                                                                                                                       aCTGCAAATGCCCCTGC‑‑TTGTGCATTTATGGAAACGTAATTTACTCAATAATCACCGTATGGTGATGGCTTCCTTTTACCCAAACTCAGCGCGGTGCTGCGCATATACGTGGAGATCAAAATGTCATTTATTAAAACTTTTTCCGGGAAg               <  1:438096/150‑1 (MQ=255)
                                                                                                                                        ccataaaTGACCCCGC‑‑AGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGc              >  2:745914/5‑150 (MQ=255)
                                                                                                                                         tGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCa             >  1:1074518/1‑150 (MQ=255)
                                                                                                                                         cataaaTGCCCCCGC‑‑AGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCa             >  2:1057895/4‑150 (MQ=255)
                                                                                                                                          gCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCAt            >  2:368345/1‑150 (MQ=255)
                                                                                                                                          gCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAAGAATCACCGGATGGTGAGGGCGTCCTTTTCCCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCAGTTATTAAAACTTTTTCCGGGAAGCAt            >  1:320007/1‑150 (MQ=255)
                                                                                                                                           cAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGGTGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCAtt           >  2:691311/1‑150 (MQ=255)
                                                                                                                                              aTGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCAttttt        <  1:452411/150‑1 (MQ=255)
                                                                                                                                                 ccccTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACGAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTttatt     >  2:485207/1‑150 (MQ=255)
                                                                                                                                                 ccccTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTttatt     >  2:1063587/1‑150 (MQ=255)
                                                                                                                                                  cccTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTttatta    <  2:316933/150‑1 (MQ=255)
                                                                                                                                                   ccTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTtattat   <  1:446417/150‑1 (MQ=255)
                                                                                                                                                    cTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATg  <  2:1259010/150‑1 (MQ=255)
                                                                                                                                                    cTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGGCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATg  <  2:1156735/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATG  >  NC_000913/2474729‑2475026

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