Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,275,123 G→T 17.0% Q7K (CAG→AAG)  yjcH ← DUF485 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,275,1230GT17.0% 865.3 / 128.6 295Q7K (CAG→AAG) yjcHDUF485 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (128/117);  new base T (33/17);  total (161/134)
Fisher's exact test for biased strand distribution p-value = 8.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 296 total aligned reads displayed.

ATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATG  >  CP009273/4274974‑4275271
                                                                                                                                                     |                                                                                                                                                    
  tCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGAt                                                                                                                                                    <  2:2413002/150‑1 (MQ=255)
    aGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATaa                                                                                                                                                  >  1:4017762/1‑150 (MQ=255)
    aGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATaa                                                                                                                                                  >  2:2953144/1‑150 (MQ=255)
    aGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATaa                                                                                                                                                  >  1:1740680/1‑150 (MQ=255)
      cggcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCCTTGTCTTCTATCCGCTGATAAAt                                                                                                                                                >  1:2654737/1‑150 (MQ=255)
       ggcgTGCCCAGCCAGCCGGGCGCGATGGCTATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATa                                                                                                                                               <  1:286587/150‑1 (MQ=255)
        gcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAg                                                                                                                                              <  2:446857/150‑1 (MQ=255)
         cgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGACATGCGCATTGTCTTCTATCCGCTGATAAATAGt                                                                                                                                             >  2:2515373/1‑150 (MQ=255)
           tGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAACTGCGCATTGTCTTCTATCCGCTGATAAATAGTGc                                                                                                                                           >  2:4410146/1‑150 (MQ=255)
               cagccagcCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCAtt                                                                                                                                       >  1:1659003/1‑150 (MQ=255)
                  ccagccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATtcat                                                                                                                                    >  1:3733245/1‑150 (MQ=255)
                   cagccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATtcata                                                                                                                                   >  1:1462235/1‑149 (MQ=255)
                    agccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATtcatca                                                                                                                                  <  2:736916/150‑1 (MQ=255)
                    agccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATtcataa                                                                                                                                  <  2:842737/150‑3 (MQ=255)
                    agccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATtcataa                                                                                                                                  >  2:1796961/1‑148 (MQ=255)
                     gccGGGCGCTAAGGCGAACAGTAAAATAAAGCCTATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATcac                                                                                                                                 <  1:25688/150‑1 (MQ=255)
                        gggCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGa                                                                                                                              >  2:3566604/1‑150 (MQ=255)
                             cgAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATtctc                                                                                                                         <  2:3627002/150‑1 (MQ=255)
                               aaGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCa                                                                                                                       >  2:1090031/1‑150 (MQ=255)
                                aGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTAATAACAGATTCTcccg                                                                                                                      >  1:772535/1‑148 (MQ=255)
                                  gCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCagag                                                                                                                    <  1:1628366/150‑1 (MQ=255)
                                  gCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCagag                                                                                                                    <  2:4018548/150‑1 (MQ=255)
                                    gATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGt                                                                                                                  <  1:2853903/150‑1 (MQ=255)
                                     aTCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTa                                                                                                                 <  2:4300931/150‑1 (MQ=255)
                                      tCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGtca                                                                                                                >  2:1221314/1‑148 (MQ=255)
                                       cAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAAt                                                                                                               <  2:2465348/150‑1 (MQ=255)
                                          tAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTa                                                                                                            <  2:488653/150‑1 (MQ=255)
                                           aaaaTAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAg                                                                                                           <  1:2032813/150‑1 (MQ=255)
                                              aTAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGTTTCTCCAGAGGTAATGTAGGGa                                                                                                        >  2:1431711/1‑150 (MQ=255)
                                               tAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGAt                                                                                                       >  2:2192152/1‑150 (MQ=255)
                                                aaaGCCGATATAAACTGCCAGCATACTAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCAGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTCGGGAtt                                                                                                      >  1:3148082/1‑150 (MQ=255)
                                                 aaGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCCCCAGAGGTAATGTAGGGATTg                                                                                                     >  2:3748320/1‑150 (MQ=255)
                                                   gCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGAGTGtt                                                                                                   >  1:2315209/1‑150 (MQ=255)
                                                    ccGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGttt                                                                                                  <  1:2634242/150‑1 (MQ=255)
                                                     cGATATAACCGGCCAGCATAATACTCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATGGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGtttt                                                                                                 <  1:4324990/150‑1 (MQ=255)
                                                     cGATATAACCGGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCAGTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGtttt                                                                                                 <  1:127372/150‑1 (MQ=255)
                                                     cGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGtttt                                                                                                 <  1:1967485/150‑1 (MQ=255)
                                                      gATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTa                                                                                                >  2:377211/1‑150 (MQ=255)
                                                       atatAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGACATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTaa                                                                                               >  2:3929996/1‑150 (MQ=255)
                                                        tataAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGACATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCCGAGGTACTGTAGGGATTGTTTTAAt                                                                                              >  2:1558918/1‑150 (MQ=255)
                                                        tataAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAAt                                                                                              >  2:3162011/1‑150 (MQ=255)
                                                          taAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTGTTTCGACTAACTCCCAGAAATGCGCATTGTCTTCTATCCTCTGATAATTAGTGCAATTCATAAAAGATTCTACAGAGGTAATGTAGGGATTGTTTTAATTc                                                                                            >  2:3086854/1‑150 (MQ=255)
                                                           aaaCTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTcc                                                                                           >  1:1040994/1‑150 (MQ=255)
                                                           aaaCTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTcc                                                                                           >  1:3322217/1‑150 (MQ=255)
                                                           aaaCTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTcc                                                                                           >  2:670795/1‑150 (MQ=255)
                                                            aaCTGCCAGCATAATACTCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTccc                                                                                          <  2:3514265/150‑1 (MQ=255)
                                                            aaCTGACAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTACCTGAAATGCGCATGGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGATGTAATGTAGGGATTGTTTTAATTccc                                                                                          <  1:2165946/150‑1 (MQ=255)
                                                                  cAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTccc                                                                                    >  2:981688/1‑150 (MQ=255)
                                                                  cAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTccc                                                                                    >  2:1855159/1‑150 (MQ=255)
                                                                   aGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCt                                                                                   <  1:1802025/150‑1 (MQ=255)
                                                                    gCATAATACTCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCtt                                                                                  <  1:2723102/150‑1 (MQ=255)
                                                                    gCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCTCAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCtt                                                                                  >  2:3809137/1‑150 (MQ=255)
                                                                    gCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGATGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCtt                                                                                  <  1:1535430/150‑1 (MQ=255)
                                                                     cATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTa                                                                                 >  2:1727932/1‑150 (MQ=255)
                                                                      ataataATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTAt                                                                                >  2:3349038/1‑150 (MQ=255)
                                                                        aataatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGGAATGTAGGGATTGTTTTAATTCACGCTCCCTTATgg                                                                              >  2:2325347/1‑150 (MQ=255)
                                                                         ataatCGACAGGATGGTGGCAAACCGTTGCCTTTTTTCGACTATCTCCCTTAAATGCGAATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATggg                                                                             <  1:1663748/150‑1 (MQ=255)
                                                                           aatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGgag                                                                           <  1:1023677/150‑1 (MQ=255)
                                                                            atCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGACATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCACCCTTATGGgaga                                                                          >  2:3039685/1‑150 (MQ=255)
                                                                              cGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAg                                                                        <  2:2597816/150‑1 (MQ=255)
                                                                               gACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAgg                                                                       >  1:3862067/1‑150 (MQ=255)
                                                                                  aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTa                                                                    <  2:3401051/150‑1 (MQ=255)
                                                                                  aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTa                                                                    >  2:2202269/1‑150 (MQ=255)
                                                                                  aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTa                                                                    >  2:3844787/1‑150 (MQ=255)
                                                                                  aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTa                                                                    <  1:3229374/150‑1 (MQ=255)
                                                                                   ggATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTaa                                                                   >  2:2694796/1‑150 (MQ=255)
                                                                                       ggtggCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCAGGGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCt                                                               <  1:3933001/150‑1 (MQ=255)
                                                                                              aaCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTaa                                                        >  1:269326/1‑150 (MQ=255)
                                                                                                 cGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAAccc                                                     <  2:2384877/150‑1 (MQ=255)
                                                                                                     gCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGc                                                 >  1:4117275/1‑150 (MQ=255)
                                                                                                        gTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGa                                              >  2:1284039/1‑150 (MQ=255)
                                                                                                        gTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGa                                              >  1:4529798/1‑150 (MQ=255)
                                                                                                          ttttCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAAt                                            <  1:559120/150‑1 (MQ=255)
                                                                                                            ttCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGt                                          <  1:3590177/150‑1 (MQ=255)
                                                                                                              cGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAg                                        <  2:1051686/150‑1 (MQ=255)
                                                                                                               gACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAgg                                       >  1:4420012/1‑150 (MQ=255)
                                                                                                               gACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCACGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTCACCCGTGCCGAATGTAgg                                       >  2:4392663/1‑150 (MQ=255)
                                                                                                                   aaCTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCgg                                   >  1:2249702/1‑150 (MQ=255)
                                                                                                                     cTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGAt                                 <  2:2555259/150‑1 (MQ=255)
                                                                                                                     cTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGAt                                 <  1:1324485/150‑1 (MQ=255)
                                                                                                                        ccTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAg                              >  2:3244595/1‑150 (MQ=255)
                                                                                                                        ccTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAg                              >  2:1014137/1‑150 (MQ=255)
                                                                                                                          tGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGc                            >  1:87023/1‑150 (MQ=255)
                                                                                                                           gAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCg                           >  2:2107340/1‑150 (MQ=255)
                                                                                                                            aaaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGt                          <  1:209947/150‑1 (MQ=255)
                                                                                                                             aaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGtt                         <  2:3732243/150‑1 (MQ=255)
                                                                                                                             aaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGtt                         <  1:2988222/150‑1 (MQ=255)
                                                                                                                                 cgcATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcg                     >  2:1447042/1‑150 (MQ=255)
                                                                                                                                  gcATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgc                    >  2:1500205/1‑150 (MQ=255)
                                                                                                                                   cATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCTCAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgcc                   >  2:2356561/1‑150 (MQ=255)
                                                                                                                                        tcttctATCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATc              >  1:3190214/1‑150 (MQ=255)
                                                                                                                                          ttctATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCg            <  1:1888502/150‑1 (MQ=255)
                                                                                                                                          ttctATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCg            >  1:942970/1‑150 (MQ=255)
                                                                                                                                             tATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCACTTATGGGAGAAGGTTAAAGCTCGGGTAACCCTTGCCGAATGTAGGCAGGATAAGTCGTTTACGCCGCATCCGGca         >  1:646247/1‑150 (MQ=255)
                                                                                                                                              aTCAGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaa        >  2:348303/1‑150 (MQ=255)
                                                                                                                                               tCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaat       >  1:2491346/1‑150 (MQ=255)
                                                                                                                                               tCCGCTTATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaat       >  1:3644987/1‑150 (MQ=255)
                                                                                                                                                   cTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatcaat   >  2:3278675/1‑150 (MQ=255)
                                                                                                                                                    tGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATg  <  2:3881757/150‑1 (MQ=255)
                                                                                                                                                    tGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATg  <  1:2393184/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATG  >  CP009273/4274974‑4275271

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