Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,471,252 T→C 15.8% N89S (AAC→AGC)  yfdK ← CPS‑53 (KpLE1) prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,471,2520TC15.8% 2852.7 / inf 1115N89S (AAC→AGC) yfdKCPS‑53 (KpLE1) prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base T (425/507);  new base C (96/79);  total (527/588)
Fisher's exact test for biased strand distribution p-value = 2.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-02

Only 100 of 1130 total aligned reads displayed.

TAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGAAAACTCAATA  >  NC_000913/2471103‑2471401
                                                                                                                                                     |                                                                                                                                                     
 aGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGtt                                                                                                                                                      <  1:9455206/150‑1 (MQ=255)
  ggCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTg                                                                                                                                                     >  2:9386267/1‑150 (MQ=255)
  ggCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTg                                                                                                                                                     >  1:8696324/1‑150 (MQ=255)
  ggCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTg                                                                                                                                                     <  2:692045/150‑1 (MQ=255)
  ggCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGCTg                                                                                                                                                     >  1:4399571/1‑150 (MQ=255)
   gCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTgg                                                                                                                                                    <  1:8008337/150‑1 (MQ=255)
     cATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGcc                                                                                                                                                  <  2:4098540/150‑1 (MQ=255)
     cATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGcc                                                                                                                                                  >  1:10746289/1‑150 (MQ=255)
        tCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGa                                                                                                                                               >  2:11347685/1‑150 (MQ=255)
        tCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGa                                                                                                                                               >  2:6573090/1‑150 (MQ=255)
        tCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGa                                                                                                                                               <  1:5915140/150‑1 (MQ=255)
         ccagatCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGAt                                                                                                                                              <  2:7408842/146‑1 (MQ=255)
         cAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGAt                                                                                                                                              <  2:11536780/150‑1 (MQ=255)
          aaTATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGAtt                                                                                                                                             <  2:3481470/150‑1 (MQ=255)
            tatCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTaa                                                                                                                                           <  2:1919742/150‑1 (MQ=255)
             atCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAAt                                                                                                                                          >  2:5691784/1‑150 (MQ=255)
              tCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATc                                                                                                                                         >  1:4657497/1‑150 (MQ=255)
               cTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCa                                                                                                                                        >  2:1538409/1‑150 (MQ=255)
                  gcgcACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAAtt                                                                                                                                     <  1:11947209/150‑1 (MQ=255)
                   cgcACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATtg                                                                                                                                    <  1:6793673/150‑1 (MQ=255)
                    gcACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATtgc                                                                                                                                   >  1:1801597/1‑150 (MQ=255)
                    gcACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATtgc                                                                                                                                   <  2:7167875/150‑1 (MQ=255)
                      aCTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATtgctg                                                                                                                                 >  1:9091708/1‑150 (MQ=255)
                       cTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATtgctgc                                                                                                                                <  1:1546722/150‑1 (MQ=255)
                        tGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATTgctgct                                                                                                                               >  1:1165922/1‑150 (MQ=255)
                              tATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCa                                                                                                                         >  2:10286745/1‑150 (MQ=255)
                              tATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAAGAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCa                                                                                                                         <  2:3864817/150‑1 (MQ=255)
                                tCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTa                                                                                                                                            <  1:4264556/129‑1 (MQ=255)
                                  gACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTc                                                                                                                     <  2:2902114/150‑1 (MQ=255)
                                    ccAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCgg                                                                                                                   >  1:3887927/1‑150 (MQ=255)
                                    ccAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCgg                                                                                                                   <  2:8731899/150‑1 (MQ=255)
                                    ccAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCgg                                                                                                                   <  1:9728787/150‑1 (MQ=255)
                                     cAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGc                                                                                                                  >  1:5373122/1‑150 (MQ=255)
                                      aGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCt                                                                                                                 >  2:4250791/1‑150 (MQ=255)
                                       gCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTg                                                                                                                >  1:5079520/1‑150 (MQ=255)
                                        cTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTgc                                                                                                               <  1:7769469/150‑1 (MQ=255)
                                         tCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTgca                                                                                                              >  1:4853992/1‑150 (MQ=255)
                                         tCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTgca                                                                                                              <  2:5077054/150‑1 (MQ=255)
                                                  tCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATAGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGt                                                                                                     >  1:9959272/1‑150 (MQ=255)
                                                   ccAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGtt                                                                                                    >  1:2007594/1‑150 (MQ=255)
                                                    cAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTc                                                                                                   <  2:5889225/150‑1 (MQ=255)
                                                    cAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTc                                                                                                   <  1:9873666/150‑1 (MQ=255)
                                                     aGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTcc                                                                                                  >  1:731012/1‑150 (MQ=255)
                                                      gATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCt                                                                                                 >  2:11231558/1‑150 (MQ=255)
                                                          aTCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCAtg                                                                                             <  1:9564601/150‑1 (MQ=255)
                                                                  acaAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtg                                                                                     <  2:11890282/150‑1 (MQ=255)
                                                                  acaAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtg                                                                                     >  2:475306/1‑150 (MQ=255)
                                                                    aaaTTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtggt                                                                                   >  2:2713998/1‑150 (MQ=255)
                                                                    aaaTTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtggt                                                                                   <  2:3721583/150‑1 (MQ=255)
                                                                    aaaTTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtggt                                                                                   <  2:8406395/150‑1 (MQ=255)
                                                                     aaTTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtggtg                                                                                  >  2:3117202/1‑150 (MQ=255)
                                                                      aTTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTtggtggtgg                                                                                 <  2:12032221/150‑1 (MQ=255)
                                                                            ttGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTc                                                                           <  1:10176520/150‑1 (MQ=255)
                                                                               cgcCTGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCTTATTATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGa                                                                        <  2:12216724/150‑1 (MQ=255)
                                                                               cgcCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGa                                                                        >  2:7724505/1‑150 (MQ=255)
                                                                               cgcCAGGTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGa                                                                        <  1:2466061/150‑1 (MQ=255)
                                                                                   aGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCAt                                                                    >  2:10706898/1‑150 (MQ=255)
                                                                                     ttCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGc                                                                  <  2:8207131/150‑1 (MQ=255)
                                                                                      tCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCTTATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCa                                                                 <  1:11053278/150‑1 (MQ=255)
                                                                                        cTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAgg                                                               >  2:616591/1‑150 (MQ=255)
                                                                                           aCCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGaaa                                                            <  2:9747207/150‑1 (MQ=255)
                                                                                             cTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAAc                                                          <  1:6551403/150‑1 (MQ=255)
                                                                                                   gacgacCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCAtttt                                                    <  1:3299195/150‑1 (MQ=255)
                                                                                                    acgacCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTtc                                                   <  2:10508935/150‑1 (MQ=255)
                                                                                                    acgacCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTtc                                                   <  1:3459605/150‑1 (MQ=255)
                                                                                                      gacCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTtctc                                                 <  1:5626114/150‑1 (MQ=255)
                                                                                                          aaTAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCAGGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGc                                             >  2:2387921/1‑150 (MQ=255)
                                                                                                           aTAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCg                                            >  2:9550698/1‑150 (MQ=255)
                                                                                                               ccGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATAc                                        >  2:1301180/1‑150 (MQ=255)
                                                                                                               ccGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATAc                                        >  2:9132463/1‑150 (MQ=255)
                                                                                                                   tttACCAGCCCATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGAt                                    <  1:2199923/150‑1 (MQ=255)
                                                                                                                      aCCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGAtttt                                 <  1:10019474/150‑1 (MQ=255)
                                                                                                                      aCCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGAtttt                                 <  1:6932356/150‑1 (MQ=255)
                                                                                                                        cAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTc                               >  1:7147840/1‑150 (MQ=255)
                                                                                                                        cAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTc                               <  1:307630/150‑1 (MQ=255)
                                                                                                                        cAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTc                               >  1:11161275/1‑150 (MQ=255)
                                                                                                                        cAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTc                               >  1:277175/1‑150 (MQ=255)
                                                                                                                        cAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTc                               <  1:12085918/150‑1 (MQ=255)
                                                                                                                        cAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTc                               <  2:6369441/150‑1 (MQ=255)
                                                                                                                         aGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTcc                              <  1:3719445/150‑1 (MQ=255)
                                                                                                                          gCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCt                             >  1:9215757/1‑150 (MQ=255)
                                                                                                                          gCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCt                             <  2:8846671/150‑1 (MQ=255)
                                                                                                                          gCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTACt                             >  2:7961663/1‑150 (MQ=255)
                                                                                                                            ccATTGTTTACTGTTCATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCttt                           >  2:4419223/1‑150 (MQ=255)
                                                                                                                                tGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcggc                       <  1:8287045/150‑1 (MQ=255)
                                                                                                                                tGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcggc                       <  1:7606842/150‑1 (MQ=255)
                                                                                                                                tGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcggc                       <  1:3821555/150‑1 (MQ=255)
                                                                                                                                 gTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcggcg                      <  2:5288539/150‑1 (MQ=255)
                                                                                                                                 gTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcggcg                      <  1:1052291/150‑1 (MQ=255)
                                                                                                                                  tttACTGTTCATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTcggcgg                     >  1:8454302/1‑150 (MQ=255)
                                                                                                                                       tGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATc                <  2:5529921/150‑1 (MQ=255)
                                                                                                                                        gTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATcc               <  1:6294626/150‑1 (MQ=255)
                                                                                                                                          tCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCgg             >  1:8904545/1‑150 (MQ=255)
                                                                                                                                           cATATAATCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGa            <  1:11320533/150‑1 (MQ=255)
                                                                                                                                            atatAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGaa           >  2:2215616/1‑150 (MQ=255)
                                                                                                                                            atatAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGaa           <  1:6126027/150‑1 (MQ=255)
                                                                                                                                              ataATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGaaaa         <  2:3311111/150‑1 (MQ=255)
                                                                                                                                                aaTCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGAAAACt       <  1:12279026/150‑1 (MQ=255)
                                                                                                                                                aaTCGCTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTACTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGAAAACt       >  2:10930790/1‑150 (MQ=255)
                                                                                                                                                   cGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGAAAACTCaa    <  1:11460138/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTCATATAATCGTTGGCCTGATTAATCAATTGCTGCTTTTTCAGTTCGGCTGCAGCAATCTGTTCCTCATGTGTTGGTGGTGGAATTTCAGACCATGCAGGAAAACCATTTTCTCCAGCGATACGGATTTTTCCTTTCGGCGGTAATCCGGAAAACTCAATA  >  NC_000913/2471103‑2471401

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