Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,359,686 C→A 27.7% E182* (GAA→TAA)  yjeH ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,359,6860CA27.7% 397.6 / 117.1 184E182* (GAA→TAA) yjeHputative transporter
Reads supporting (aligned to +/- strand):  ref base C (62/71);  new base A (29/22);  total (91/93)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 186 total aligned reads displayed.

CTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAC  >  CP009273/4359539‑4359833
                                                                                                                                                   |                                                                                                                                                   
cTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTc                                                                                                                                                     <  2:2334598/148‑1 (MQ=255)
 tAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCg                                                                                                                                                    <  1:1136700/148‑1 (MQ=255)
   aTCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACTCAGCAAATAACCCGGTAAGTTCGat                                                                                                                                                  <  1:1027307/148‑1 (MQ=255)
   aTCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGat                                                                                                                                                  <  2:1035508/148‑1 (MQ=255)
   aTCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGat                                                                                                                                                  <  1:879573/148‑1 (MQ=255)
    tCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGata                                                                                                                                                 >  2:913134/1‑148 (MQ=255)
    tCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGata                                                                                                                                                 >  2:816765/1‑148 (MQ=255)
    tCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGata                                                                                                                                                 <  1:1683109/148‑1 (MQ=255)
        gCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTAc                                                                                                                                             <  2:491862/148‑1 (MQ=255)
           agcagcagACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAg                                                                                                                                          <  1:1276103/148‑1 (MQ=255)
           agcagcagACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAg                                                                                                                                          >  2:3223161/1‑148 (MQ=255)
            gcagcagACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAgg                                                                                                                                         >  2:3032293/1‑148 (MQ=255)
                  gACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCgg                                                                                                                                   >  1:1518039/1‑148 (MQ=255)
                  gACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCgg                                                                                                                                   >  1:773101/1‑148 (MQ=255)
                     cAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGaaa                                                                                                                                >  1:2719807/1‑148 (MQ=255)
                      aaTCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAg                                                                                                                               >  2:2600076/1‑148 (MQ=255)
                      aaTCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAg                                                                                                                               <  2:2866488/148‑1 (MQ=255)
                       atcatcAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAgg                                                                                                                              >  1:1183492/1‑148 (MQ=255)
                         catccAAGCACGAGGAAAATCACGCTCGGGAGTTTTAAAGTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAAGAACCCGGGAAGTTCGATATTACCAGGTGCCGGAAAgggg                                                                                                                            <  2:2871191/148‑1 (MQ=255)
                         catcaAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAgggg                                                                                                                            <  2:1092972/148‑1 (MQ=255)
                           tcaAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGat                                                                                                                          <  1:2987258/148‑1 (MQ=255)
                           tcaAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGat                                                                                                                          <  1:353619/148‑1 (MQ=255)
                            caAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGata                                                                                                                         <  2:47991/148‑1 (MQ=255)
                              aaGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATAtt                                                                                                                       <  1:2879774/148‑1 (MQ=255)
                                gCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCg                                                                                                                     >  2:66498/1‑148 (MQ=255)
                                gCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCg                                                                                                                     >  2:1227793/1‑148 (MQ=255)
                                 cACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGACAGGGGATATTCGc                                                                                                                    >  2:3125600/1‑148 (MQ=255)
                                  aCGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCa                                                                                                                   <  1:1840926/148‑1 (MQ=255)
                                    gAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAgg                                                                                                                 <  1:1345462/148‑1 (MQ=255)
                                     aGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGt                                                                                                                >  1:1007471/1‑148 (MQ=255)
                                       gAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGttt                                                                                                              >  2:2533942/1‑148 (MQ=255)
                                        aaaaTCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTg                                                                                                             >  1:981684/1‑148 (MQ=255)
                                          aaTCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGat                                                                                                           <  2:1140900/148‑1 (MQ=255)
                                          aaTCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGat                                                                                                           >  1:2186723/1‑148 (MQ=255)
                                            tCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGatat                                                                                                         >  2:3197498/1‑148 (MQ=255)
                                               cGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATcgc                                                                                                      >  1:377044/1‑148 (MQ=255)
                                                    tGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATcgcccgcc                                                                                                 <  1:2316253/148‑1 (MQ=255)
                                                         ttttAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCccacca                                                                                            >  1:2258308/1‑148 (MQ=255)
                                                         ttttAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCccacca                                                                                            >  1:2258323/1‑148 (MQ=255)
                                                          tttAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAg                                                                                           <  1:3197501/148‑1 (MQ=255)
                                                          tttAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAg                                                                                           >  1:3428798/1‑148 (MQ=255)
                                                           ttAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCTTTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGa                                                                                          <  2:2096655/148‑1 (MQ=255)
                                                           ttAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGa                                                                                          <  1:1396979/148‑1 (MQ=255)
                                                            tAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGAt                                                                                         <  2:1174329/148‑1 (MQ=255)
                                                            tAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGAt                                                                                         <  1:2080862/148‑1 (MQ=255)
                                                              aaTTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAg                                                                                       <  1:1900288/148‑1 (MQ=255)
                                                                 tCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGa                                                                                    >  1:51234/1‑148 (MQ=255)
                                                                 tCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGa                                                                                    >  1:3268971/1‑148 (MQ=255)
                                                                   cGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATACCCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACa                                                                                  <  1:3506290/148‑1 (MQ=255)
                                                                        cGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAg                                                                             >  1:1586765/1‑148 (MQ=255)
                                                                        cGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAg                                                                             >  1:468568/1‑148 (MQ=255)
                                                                            aTGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGcgcg                                                                         >  2:2184837/1‑148 (MQ=255)
                                                                                gCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGcgacga                                                                     <  1:854923/148‑1 (MQ=255)
                                                                                   aaTGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATACCCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATaa                                                                  <  2:1412486/148‑1 (MQ=255)
                                                                                   aaTGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATaa                                                                  <  2:2220019/148‑1 (MQ=255)
                                                                                    aTGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAg                                                                 >  2:3242310/1‑148 (MQ=255)
                                                                                     tGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGt                                                                >  2:2305470/1‑148 (MQ=255)
                                                                                     tGCATCAAGACCGACAAAACACCAGAACATCACTGATAAAGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGTGTATAATAGCAGGTTTGATATCGCAAGCCCACCAGATAGCGACAAAAAGAGCGACGATAAGt                                                                >  2:2904273/1‑148 (MQ=255)
                                                                                      gCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTc                                                               <  2:2250881/148‑1 (MQ=255)
                                                                                       ccTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTcc                                                              >  1:758927/1‑148 (MQ=255)
                                                                                       ccTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTcc                                                              >  1:2819717/1‑148 (MQ=255)
                                                                                         tCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCgg                                                            >  2:685657/1‑148 (MQ=255)
                                                                                         tCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCgg                                                            >  1:3553535/1‑148 (MQ=255)
                                                                                         tCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCgg                                                            <  2:1495944/148‑1 (MQ=255)
                                                                                         tCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCgg                                                            >  2:1637419/1‑148 (MQ=255)
                                                                                            aGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCaa                                                         <  2:530808/148‑1 (MQ=255)
                                                                                              aCCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATa                                                       <  1:2614863/148‑1 (MQ=255)
                                                                                              aCCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATa                                                       <  1:830339/148‑1 (MQ=255)
                                                                                              aCCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATa                                                       <  1:243779/148‑1 (MQ=255)
                                                                                              aCCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATa                                                       <  2:2034454/148‑1 (MQ=255)
                                                                                                cGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAAc                                                     >  2:2081514/1‑148 (MQ=255)
                                                                                                cGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAAc                                                     >  2:283110/1‑148 (MQ=255)
                                                                                                  aCAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACgg                                                   <  2:3452508/148‑1 (MQ=255)
                                                                                                  aCAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACgg                                                   >  2:3475994/1‑148 (MQ=255)
                                                                                                   cAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGt                                                  <  1:3193078/148‑1 (MQ=255)
                                                                                                   cAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGt                                                  <  1:703325/148‑1 (MQ=255)
                                                                                                        cacCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTa                                             >  1:2516948/1‑148 (MQ=255)
                                                                                                          ccAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGa                                           >  2:2089144/1‑148 (MQ=255)
                                                                                                          ccAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGa                                           >  2:1206868/1‑148 (MQ=255)
                                                                                                           cAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGAt                                          <  1:2694766/148‑1 (MQ=255)
                                                                                                            aGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGAtt                                         <  2:1569052/148‑1 (MQ=255)
                                                                                                               aCATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGc                                      >  2:2547537/1‑148 (MQ=255)
                                                                                                                   cACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTg                                  <  1:675519/148‑1 (MQ=255)
                                                                                                                   cACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTg                                  <  1:2908409/148‑1 (MQ=255)
                                                                                                                    aCTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTgg                                 <  2:2698387/148‑1 (MQ=255)
                                                                                                                     cTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGa                                >  1:20395/1‑148 (MQ=255)
                                                                                                                         tAACGCAGCAAATAACCCGGTAAGTTAGATATTACCGGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTg                            >  2:1722843/1‑148 (MQ=255)
                                                                                                                            cGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATAGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCa                         >  1:3121029/1‑148 (MQ=255)
                                                                                                                             gcagcaAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCAc                        >  1:2339582/1‑148 (MQ=255)
                                                                                                                             gcagcaAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCAc                        >  1:2340265/1‑148 (MQ=255)
                                                                                                                              cagcaAATAACCCGGTCAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCAcc                       >  2:1408337/1‑148 (MQ=255)
                                                                                                                              cagcaAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCAcc                       >  2:335649/1‑148 (MQ=255)
                                                                                                                                    aTAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGt                 >  2:79598/1‑148 (MQ=255)
                                                                                                                                      aaCGCGGTAAGTTAGATATTACCAGGTGCCGGATAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTAc               <  2:2376010/148‑1 (MQ=255)
                                                                                                                                      aaCCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTAc               <  2:3270885/148‑1 (MQ=255)
                                                                                                                                        cccGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCg             <  2:377044/148‑1 (MQ=255)
                                                                                                                                             tAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGAtata        <  1:1759836/148‑1 (MQ=255)
                                                                                                                                                gTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACca     <  1:1621648/148‑1 (MQ=255)
                                                                                                                                                  tCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACcaca   >  2:3089420/1‑148 (MQ=255)
                                                                                                                                                  tCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACcaca   >  2:1543794/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAC  >  CP009273/4359539‑4359833

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: