Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,761,733 T→A 31.0% V222E (GTA→GAA)  yfjP → CP4‑57 prophage; 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,761,7330TA31.0% 397.2 / 263.4 264V222E (GTA→GAA) yfjPCP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (90/91);  new base A (29/51);  total (122/142)
Fisher's exact test for biased strand distribution p-value = 5.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 283 total aligned reads displayed.

CAAAATATCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGC  >  CP009273/2761585‑2761882
                                                                                                                                                    |                                                                                                                                                     
 aaaaTATCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGTTTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGgaag                                                                                                                                                     >  2:4934300/1‑150 (MQ=255)
 aaaaTATCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGgtag                                                                                                                                                     >  1:4296135/1‑150 (MQ=255)
 aaaaTATCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGgtag                                                                                                                                                     >  2:3875706/1‑150 (MQ=255)
 aaaaTATCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGgtag                                                                                                                                                     >  1:792601/1‑150 (MQ=255)
 aaaaTATCAGCCGCAAAATCTGCCTGCTGCATCAGCAAATCCAGCCCGTGAATACGGTGAGTGCCGTGTCGGTTCGCCTGCAGTGGGGACAAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGgtag                                                                                                                                                     >  1:1735276/1‑150 (MQ=21)
  aaaTATCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGgtagt                                                                                                                                                    <  2:3365124/150‑1 (MQ=255)
  aaaTATCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGgaagt                                                                                                                                                    <  2:3119537/150‑1 (MQ=255)
   aaTATCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGgtagta                                                                                                                                                   <  1:4683158/150‑1 (MQ=255)
   aaTATCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGgtagta                                                                                                                                                   >  2:2326001/1‑150 (MQ=255)
    atatCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGtagtag                                                                                                                                                  <  2:1267110/150‑1 (MQ=255)
     tatCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAgc                                                                                                                                                 <  2:2770539/150‑1 (MQ=255)
     gatCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAgc                                                                                                                                                 <  1:4657446/149‑1 (MQ=255)
      atCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAgcg                                                                                                                                                >  1:2424215/1‑150 (MQ=255)
        cAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCt                                                                                                                                              >  2:1302260/1‑150 (MQ=255)
        cAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCt                                                                                                                                              >  1:1390809/1‑150 (MQ=255)
         aGCCGCAAAATCTGCCTGCTGCAGGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTa                                                                                                                                             <  2:1649207/150‑1 (MQ=255)
              cAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCa                                                                                                                                        >  2:2844216/1‑150 (MQ=255)
              cAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCa                                                                                                                                        >  2:907255/1‑150 (MQ=255)
               aaaaTCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCaa                                                                                                                                       <  2:1145599/150‑1 (MQ=255)
                 aaTCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACa                                                                                                                                     >  2:670408/1‑150 (MQ=255)
                   tCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATc                                                                                                                                   >  2:4645572/1‑150 (MQ=255)
                     tgcctgcTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATccc                                                                                                                                 >  2:3058339/1‑150 (MQ=255)
                     tgcctgcTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATccc                                                                                                                                 >  2:4429728/1‑150 (MQ=255)
                       cctgcTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCtt                                                                                                                               <  2:2725292/150‑1 (MQ=255)
                       cctgcTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCtt                                                                                                                               <  2:3450191/150‑1 (MQ=255)
                         tgctgcATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTc                                                                                                                             <  2:420062/150‑1 (MQ=255)
                             gcATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTTGCGCTACTGCAACATCCCTTTCGaac                                                                                                                         <  1:3450714/150‑1 (MQ=255)
                               aTGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGaacaa                                                                                                                       >  2:2689677/1‑150 (MQ=255)
                                tGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGaacaac                                                                                                                      <  1:1572672/150‑1 (MQ=255)
                                tGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTCCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGaacaac                                                                                                                      >  1:681049/1‑150 (MQ=255)
                                 gAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACg                                                                                                                     >  1:5019520/1‑150 (MQ=255)
                                  aGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACgn                                                                                                                    <  1:590132/150‑2 (MQ=255)
                                  aGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACgg                                                                                                                    >  1:4628005/1‑150 (MQ=255)
                                    cTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTa                                                                                                                  <  1:3122341/150‑1 (MQ=255)
                                    cTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTa                                                                                                                  >  2:2447181/1‑150 (MQ=255)
                                       ttCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGcc                                                                                                               >  1:4414734/1‑150 (MQ=255)
                                       ttCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGcc                                                                                                               >  1:213338/1‑150 (MQ=255)
                                        tCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGTCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACTACGGTAGccc                                                                                                              <  1:493933/150‑1 (MQ=255)
                                        tCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGccc                                                                                                              <  2:1817481/150‑1 (MQ=255)
                                        tCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGccc                                                                                                              <  2:1228737/150‑1 (MQ=255)
                                         ccAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCg                                                                                                             <  2:2129116/150‑1 (MQ=255)
                                           tgCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGAGTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCggg                                                                                                           >  2:2536887/2‑150 (MQ=255)
                                           aGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCggg                                                                                                           >  1:1178141/1‑150 (MQ=255)
                                           aGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCggg                                                                                                           >  1:2608282/1‑150 (MQ=255)
                                           aGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCggg                                                                                                           >  1:2821996/1‑150 (MQ=255)
                                           aGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCggg                                                                                                           >  1:2697226/1‑150 (MQ=255)
                                            gCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGa                                                                                                          >  2:3191133/1‑150 (MQ=255)
                                            gCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGa                                                                                                          <  2:3984608/150‑1 (MQ=255)
                                               cGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCa                                                                                                       >  1:1048328/1‑150 (MQ=255)
                                                       cGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTg                                                                                               <  2:3474933/150‑1 (MQ=255)
                                                       cGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTg                                                                                               <  2:642007/150‑1 (MQ=255)
                                                       cGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTg                                                                                               <  2:1616729/150‑1 (MQ=255)
                                                         gtgtgtGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGAc                                                                                             <  1:2870421/150‑1 (MQ=255)
                                                         gtgtgtGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGAc                                                                                             <  1:3334625/150‑1 (MQ=255)
                                                            tgtgCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAt                                                                                          >  1:102480/1‑150 (MQ=255)
                                                             gtgCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAtt                                                                                         >  2:4841592/1‑150 (MQ=255)
                                                              tgCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAttt                                                                                        >  2:214593/1‑150 (MQ=255)
                                                                  gtgtCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGt                                                                                    <  2:2558823/150‑1 (MQ=255)
                                                                   tgtCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTg                                                                                   <  2:1502026/150‑1 (MQ=255)
                                                                    gtCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGa                                                                                  >  1:1290504/1‑150 (MQ=255)
                                                                    gtCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGa                                                                                  >  1:2564391/1‑150 (MQ=255)
                                                                    gtCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGa                                                                                  >  2:4767305/1‑150 (MQ=255)
                                                                      cGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGaaa                                                                                <  2:2308982/150‑1 (MQ=255)
                                                                        gTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAcc                                                                              <  1:3875706/150‑1 (MQ=255)
                                                                               tGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTg                                                                       <  2:370804/150‑1 (MQ=255)
                                                                                 cAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGcc                                                                     <  1:4175846/150‑1 (MQ=255)
                                                                                   gTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCAt                                                                   <  1:2684790/150‑1 (MQ=255)
                                                                                   gTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCAt                                                                   <  1:1871350/150‑1 (MQ=255)
                                                                                   gTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCAt                                                                   <  1:917157/150‑1 (MQ=255)
                                                                                     ggggACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                                                                 >  1:4865333/1‑150 (MQ=255)
                                                                                       ggACTCAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                                                               <  1:3439849/150‑1 (MQ=255)
                                                                                       ggACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                                                               <  1:890115/150‑1 (MQ=255)
                                                                                          cTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGAt                                                            <  2:4854043/150‑1 (MQ=255)
                                                                                            aaGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATAc                                                          <  2:777149/150‑1 (MQ=255)
                                                                                              gggTGATGGCAGAGCGGATGATTAAGTTTCTGCCCCGTGAGTCCAGCAGCCCCGAAGTAGCGCTACTGCAACATCCCTTACGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGAAGGTGCCATACTCGATACgg                                                        >  1:33338/1‑150 (MQ=255)
                                                                                              gggTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACgg                                                        <  1:1609353/150‑1 (MQ=255)
                                                                                              gggTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACgg                                                        >  1:33749/1‑150 (MQ=255)
                                                                                                 tGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTaa                                                     >  1:1510319/1‑150 (MQ=255)
                                                                                                  gATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAg                                                    >  2:3001686/1‑150 (MQ=255)
                                                                                                   aTGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGc                                                   <  2:3395396/150‑1 (MQ=255)
                                                                                                   aTGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGc                                                   <  2:138852/150‑1 (MQ=255)
                                                                                                    tGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCa                                                  >  1:3194926/1‑150 (MQ=255)
                                                                                                         gagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCttt                                             <  1:4114535/150‑1 (MQ=255)
                                                                                                             ggATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcccc                                         <  2:1985471/150‑1 (MQ=255)
                                                                                                               atgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCtt                                       <  1:2011894/150‑1 (MQ=255)
                                                                                                                  atTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAtt                                    >  2:3114693/1‑150 (MQ=255)
                                                                                                                  atTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAtt                                    <  2:1723059/150‑1 (MQ=255)
                                                                                                                  atTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAtt                                    >  2:1705265/1‑150 (MQ=255)
                                                                                                                      aGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCg                                <  2:4296135/150‑1 (MQ=255)
                                                                                                                               cgcgTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGc                       <  1:2313500/150‑1 (MQ=255)
                                                                                                                                  gTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGa                    <  2:3394761/150‑1 (MQ=255)
                                                                                                                                   tGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGAc                   <  2:1094427/150‑1 (MQ=255)
                                                                                                                                       gCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGAGACGATTTCTGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTATCCCCTTATTCCCGCCCCGGTGCGGACCata               >  1:990413/1‑149 (MQ=255)
                                                                                                                                            cagcCCGGTCGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCa          >  1:2976355/1‑150 (MQ=255)
                                                                                                                                            cagcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCa          >  2:2891242/1‑150 (MQ=255)
                                                                                                                                            cagcCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCa          <  1:2676302/150‑1 (MQ=255)
                                                                                                                                             agcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAg         >  1:3085985/1‑150 (MQ=255)
                                                                                                                                                 cGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCg     <  1:290365/150‑1 (MQ=255)
                                                                                                                                                    tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGc  >  1:3537280/1‑150 (MQ=255)
                                                                                                                                                    tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACATGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGc  >  2:405246/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAAAATATCAGCCGCAAAATCTGCCTGCTGCATGAGCTATTCCAGCCCGTGCATCCGGTGTGTGCCGTGTCGGTTCGCCTGCAGTGGGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGC  >  CP009273/2761585‑2761882

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: