Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,474,881 G→A 5.3% intergenic (+25/‑103) 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,8810GA 5.3% 260.2 / 4.3 94intergenic (+25/‑103)yfdQ/yfdRCPS‑53 (KpLE1) prophage; uncharacterized protein/CPS‑53 (KpLE1) prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base G (38/50);  new base A (2/3);  total (40/54)
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.04e-01

GGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACA  >  NC_000913/2474733‑2475029
                                                                                                                                                      |                                                                                                                                                    
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                                                                                                                                    <  1:1139580/150‑1 (MQ=255)
                 aGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACg                                                                                                                                    <  2:1179375/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                                                                                                     <  1:1066172/150‑1 (MQ=255)
                                                cGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAgg                                                                                                     <  2:9176/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)
                                                                                                                                                   c‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAAAAAAATGTCATTTATTAAATCTTTTTCCGGGAAGCATTTTTATTATGACa  >  2:778826/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGC‑‑GGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACA  >  NC_000913/2474733‑2475029

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: