Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,763,057 C→A 5.8% intergenic (+212/‑5) yfjQ → / → yfjR CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,763,0570CA 5.8% 1142.1 / 29.2 298intergenic (+212/‑5)yfjQ/yfjRCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (169/111);  new base A (15/3);  total (184/114)
Fisher's exact test for biased strand distribution p-value = 7.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 307 total aligned reads displayed.

TGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCACGCTGCGGCGTGATTTTC  >  CP009273/2762909‑2763206
                                                                                                                                                    |                                                                                                                                                     
   cTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAAtctctc                                                                                                                                                   <  1:1857762/150‑1 (MQ=255)
   cTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCGCTTTCAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAAtctctc                                                                                                                                                   >  1:2337200/1‑150 (MQ=255)
     cAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCAt                                                                                                                                                 <  1:3414769/150‑1 (MQ=255)
      aaTCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTGTCCTGGAGAATCTCTCATg                                                                                                                                                >  1:191435/1‑150 (MQ=255)
         cccccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGaca                                                                                                                                             <  1:3401755/150‑1 (MQ=255)
          ccccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGacac                                                                                                                                            >  2:3058215/1‑150 (MQ=255)
             cGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAgg                                                                                                                                         <  1:2205007/150‑1 (MQ=255)
              ggTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGc                                                                                                                                        <  1:2470341/150‑1 (MQ=255)
                ttATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAg                                                                                                                                      <  1:568549/150‑1 (MQ=255)
                 tataGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGa                                                                                                                                     >  2:1240891/1‑150 (MQ=255)
                       ttttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAAcgccg                                                                                                                               <  2:1455662/150‑1 (MQ=255)
                         ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACgccgcc                                                                                                                             <  2:4149981/150‑1 (MQ=255)
                          tAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCa                                                                                                                            >  2:3071151/1‑150 (MQ=255)
                           aaCCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCAt                                                                                                                           >  1:4151574/1‑150 (MQ=255)
                           aaCCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCAt                                                                                                                           <  1:4628190/150‑1 (MQ=255)
                            aCCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATg                                                                                                                          <  2:1024576/150‑1 (MQ=255)
                             cccccGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGa                                                                                                                         >  2:2254450/1‑150 (MQ=255)
                              ccccGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGAAACAGGCAGAACGCCGCCATGAc                                                                                                                        >  1:1145921/1‑150 (MQ=255)
                                ccGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCg                                                                                                                      <  2:4298772/150‑1 (MQ=255)
                                ccGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCg                                                                                                                      <  1:3740160/150‑1 (MQ=255)
                                    tACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCggctg                                                                                                                  >  2:2863253/1‑150 (MQ=255)
                                       aTCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTCATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCggctggct                                                                                                               >  1:1300650/1‑150 (MQ=255)
                                       aTCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCggctggct                                                                                                               >  1:216712/1‑150 (MQ=255)
                                        tCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGgctggctg                                                                                                              >  2:2130335/1‑150 (MQ=255)
                                          tGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtc                                                                                                            <  2:1342262/150‑1 (MQ=255)
                                            gCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcag                                                                                                          >  2:432740/1‑150 (MQ=255)
                                            gCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTgactgtcag                                                                                                          >  2:1105681/1‑150 (MQ=255)
                                             cTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcagg                                                                                                         >  2:1215687/1‑150 (MQ=255)
                                               tttGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcaggct                                                                                                       <  1:3041939/150‑1 (MQ=255)
                                               tttGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcaggctg                                                                                                      <  2:3825807/150‑1 (MQ=255)
                                                 tGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcaggctgt                                                                                                     >  1:809273/1‑150 (MQ=255)
                                                 tGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcaggctgt                                                                                                     <  1:1629881/150‑1 (MQ=255)
                                                    aGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCAc                                                                                                  >  1:1085257/1‑150 (MQ=255)
                                                    aGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCAc                                                                                                  >  1:3835908/1‑150 (MQ=255)
                                                     gAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACt                                                                                                 >  1:3173678/1‑150 (MQ=255)
                                                       aaTAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGa                                                                                               >  1:1415545/1‑150 (MQ=255)
                                                         tAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATa                                                                                             >  1:2881577/1‑150 (MQ=255)
                                                           aaaaaTAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAAt                                                                                           >  2:599633/1‑150 (MQ=255)
                                                             aaaTAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCa                                                                                         <  1:2674145/150‑1 (MQ=255)
                                                             aaaTAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCa                                                                                         >  2:344821/1‑150 (MQ=255)
                                                              aaTAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAg                                                                                        >  1:2041943/1‑150 (MQ=255)
                                                                tAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGcc                                                                                      >  1:3637994/1‑150 (MQ=255)
                                                                   tttCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATATCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTc                                                                                   >  1:272722/1‑150 (MQ=255)
                                                                   tttCTGCGTTGTCCATACCCTGTCCGCCCCACTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGATAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTc                                                                                   >  1:1785020/1‑150 (MQ=255)
                                                                     tCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTg                                                                                 <  1:1032055/150‑1 (MQ=255)
                                                                      cTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTgg                                                                                >  2:3168913/1‑150 (MQ=255)
                                                                      cTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTgg                                                                                >  1:4356214/1‑150 (MQ=255)
                                                                      cTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTgg                                                                                <  1:3624388/150‑1 (MQ=255)
                                                                         cGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTg                                                                             >  1:2816555/1‑150 (MQ=255)
                                                                           ttGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCg                                                                           <  1:374507/150‑1 (MQ=255)
                                                                            tGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCgg                                                                          >  2:1182706/1‑150 (MQ=255)
                                                                            tGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCgg                                                                          <  2:2418016/150‑1 (MQ=255)
                                                                             gTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAAACACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATATCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCggg                                                                         >  2:2188834/1‑150 (MQ=255)
                                                                               ccATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTg                                                                       >  2:2484413/1‑150 (MQ=255)
                                                                               ccATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTg                                                                       >  2:273121/1‑150 (MQ=255)
                                                                                cATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGa                                                                      <  2:3112711/150‑1 (MQ=255)
                                                                                cATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGa                                                                      >  1:294938/1‑150 (MQ=255)
                                                                                  tACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATATCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGaaa                                                                    >  2:3028011/1‑150 (MQ=255)
                                                                                   aCCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAAc                                                                   >  2:104174/1‑150 (MQ=255)
                                                                                     ccTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGt                                                                 >  1:2233752/1‑150 (MQ=255)
                                                                                        gTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGACTCTCTCATGACACAGGCAGAACGCCGCCATGACCGTCTGTCTGTCAGTCTGTCACTGATAATAAGCCGACTGGTTGCGGGTGAAACGTTGa                                                              >  2:250921/1‑150 (MQ=17)
                                                                                         tCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAg                                                             >  1:2717376/1‑150 (MQ=255)
                                                                                          ccGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGc                                                            <  1:4256578/150‑1 (MQ=255)
                                                                                           cGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCg                                                           >  2:3070299/1‑150 (MQ=255)
                                                                                           cGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCg                                                           <  1:1511101/150‑1 (MQ=255)
                                                                                           cGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCg                                                           <  1:3905116/150‑1 (MQ=255)
                                                                                            gCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGt                                                          <  1:2611421/150‑1 (MQ=255)
                                                                                                      ttAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTg                                                <  1:3721902/150‑1 (MQ=255)
                                                                                                       tAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTgg                                               >  2:1896195/1‑150 (MQ=255)
                                                                                                          aGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTggcgg                                            <  1:4453364/150‑1 (MQ=255)
                                                                                                          aGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTggcgg                                            >  1:1470575/1‑150 (MQ=255)
                                                                                                            tAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCt                                          >  2:3003589/1‑150 (MQ=255)
                                                                                                            tAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCt                                          >  2:2395964/1‑150 (MQ=255)
                                                                                                              aTCACATCATTTTCAGTCAGTTAACTGTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGa                                        <  2:3710281/150‑1 (MQ=255)
                                                                                                               tCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAg                                       <  1:2768171/150‑1 (MQ=255)
                                                                                                                cacaTCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGt                                      >  2:198951/1‑150 (MQ=255)
                                                                                                                  catcatTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGttt                                    >  1:3857755/1‑150 (MQ=255)
                                                                                                                   atcatTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTg                                   <  1:4051488/150‑1 (MQ=255)
                                                                                                                   atcatGTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACCGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTg                                   <  1:4058385/150‑1 (MQ=255)
                                                                                                                    tcatTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTgg                                  <  1:1824525/150‑1 (MQ=255)
                                                                                                                       ttttCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCgt                               <  2:1509926/150‑1 (MQ=255)
                                                                                                                         ttCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCgtgt                             >  1:2367283/1‑150 (MQ=255)
                                                                                                                          tccgtcagtTAACGTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTc                            <  1:358279/147‑1 (MQ=255)
                                                                                                                          tcagtcagTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTc                            >  1:150837/1‑150 (MQ=255)
                                                                                                                          tcagtcagTTAACTGTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTc                            <  1:1706779/150‑1 (MQ=255)
                                                                                                                            agtcagtTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAg                          >  2:349924/1‑150 (MQ=255)
                                                                                                                              tcagtTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTg                        <  1:4396195/150‑1 (MQ=255)
                                                                                                                                agtTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTgcg                      >  1:1381800/1‑150 (MQ=255)
                                                                                                                                 gtTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTgcgc                     >  1:58455/1‑150 (MQ=255)
                                                                                                                                  ttAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCa                    >  2:184897/1‑150 (MQ=255)
                                                                                                                                  ttAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCa                    >  2:4404427/1‑150 (MQ=255)
                                                                                                                                    aaCTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCACg                  >  2:859823/1‑150 (MQ=255)
                                                                                                                                      cTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCACGCt                <  2:1152278/150‑1 (MQ=255)
                                                                                                                                        ttCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCACGCTgc              >  1:702663/1‑150 (MQ=255)
                                                                                                                                            tGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCACGCTgcggcg          >  1:716773/1‑150 (MQ=255)
                                                                                                                                             ggAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCACGCTGCGGCGt         >  1:3669300/1‑150 (MQ=255)
                                                                                                                                             ggAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCACGCTGCGGCGt         >  1:3096989/1‑150 (MQ=255)
                                                                                                                                             ggAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCACGCTGCGGCGt         >  2:512781/1‑150 (MQ=255)
                                                                                                                                              gagaATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCACGCTGCGGCGTg        <  2:4004002/150‑1 (MQ=255)
                                                                                                                                               agaATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCACGCTGCGGCGTGa       >  1:311707/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGTGAAACGTTGAGCGTGCGCAAACTGGCGGCTGAGTTTGGCGTGTCAGTGCGCACGCTGCGGCGTGATTTTC  >  CP009273/2762909‑2763206

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: