Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,819,892 C→A 47.5% intergenic (‑10/+320) spy ← / ← astE periplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,819,8920CA47.5% 149.1 / 307.1 240intergenic (‑10/+320)spy/astEperiplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (60/66);  new base A (61/53);  total (121/119)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 254 total aligned reads displayed.

CGGACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAA  >  CP009273/1819745‑1820037
                                                                                                                                                   |                                                                                                                                                 
  gACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtact                                                                                                                                                 >  1:1291210/1‑148 (MQ=255)
   aCCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtcctt                                                                                                                                                >  1:278529/1‑148 (MQ=255)
   aCCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtcctt                                                                                                                                                <  1:2423647/148‑1 (MQ=255)
     cGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtccttcc                                                                                                                                              <  1:2959472/148‑1 (MQ=255)
     cGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtccttcc                                                                                                                                              <  2:1687988/148‑1 (MQ=255)
      gAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATAtacttcct                                                                                                                                             <  2:2627829/148‑1 (MQ=255)
        aCTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCttt                                                                                                                                           <  2:135281/148‑1 (MQ=255)
         cTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTc                                                                                                                                          >  1:3550858/1‑148 (MQ=255)
         cTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTc                                                                                                                                          <  1:556146/148‑1 (MQ=255)
          ttGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCa                                                                                                                                         >  1:3482160/1‑148 (MQ=255)
          ttGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCa                                                                                                                                         >  2:1283507/1‑148 (MQ=255)
          ttGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCa                                                                                                                                         >  2:2149389/1‑148 (MQ=255)
          ttGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCa                                                                                                                                         >  2:3069105/1‑148 (MQ=255)
           tGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAg                                                                                                                                        >  2:816154/1‑148 (MQ=255)
             ccTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGtt                                                                                                                                      >  2:260414/1‑148 (MQ=255)
             ccTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGtt                                                                                                                                      >  2:2060153/1‑148 (MQ=255)
              cTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGtta                                                                                                                                     >  2:2475662/1‑148 (MQ=255)
               tttGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGttat                                                                                                                                    <  2:522319/148‑1 (MQ=255)
               tttGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGttat                                                                                                                                    <  2:3398790/148‑1 (MQ=255)
                ttGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGttatt                                                                                                                                   >  2:152049/1‑148 (MQ=255)
                 tgtgGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGttattt                                                                                                                                  <  2:2501745/148‑1 (MQ=255)
                  gtggtgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGttattta                                                                                                                                 <  2:736481/148‑1 (MQ=255)
                    ggtgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTtatttatt                                                                                                                               <  2:1841307/148‑1 (MQ=255)
                    ggtgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTtatttatt                                                                                                                               <  1:1361248/148‑1 (MQ=255)
                     gtgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATttatta                                                                                                                              <  2:3255906/148‑1 (MQ=255)
                      tgCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTAc                                                                                                                             >  1:1534767/1‑148 (MQ=255)
                       gCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACg                                                                                                                            <  2:1644017/148‑1 (MQ=255)
                       gCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACg                                                                                                                            <  1:2129116/148‑1 (MQ=255)
                        catcatCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACgg                                                                                                                           >  2:235890/1‑148 (MQ=255)
                         atcatcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGc                                                                                                                          <  2:2305820/148‑1 (MQ=255)
                           catcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCtt                                                                                                                        >  2:2405183/1‑148 (MQ=255)
                            atcGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCttt                                                                                                                       <  2:364508/148‑1 (MQ=255)
                              cGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCt                                                                                                                     >  1:2839584/1‑148 (MQ=255)
                                gCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTa                                                                                                                   <  2:2186174/148‑1 (MQ=255)
                                     gcgTCAGCCGGTGATGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAg                                                                                                              >  2:1390713/1‑148 (MQ=255)
                                          aGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGc                                                                                                         <  1:1152243/148‑1 (MQ=255)
                                              ggTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCgtg                                                                                                     >  2:1581355/1‑148 (MQ=255)
                                               gTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCgtgt                                                                                                    <  2:835085/148‑1 (MQ=255)
                                               gTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCgtgt                                                                                                    >  1:1373377/1‑148 (MQ=255)
                                                tgctgcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGtt                                                                                                   <  2:284929/148‑1 (MQ=255)
                                                 gctgcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTg                                                                                                  <  2:19510/148‑1 (MQ=255)
                                                 gctgcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTg                                                                                                  <  2:19399/148‑1 (MQ=255)
                                                    gcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACg                                                                                               <  1:669168/148‑1 (MQ=255)
                                                    gcGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACg                                                                                               <  1:1784383/148‑1 (MQ=255)
                                                      ggTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACgag                                                                                             <  2:2513730/148‑1 (MQ=255)
                                                      ggTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACgag                                                                                             <  1:1321479/148‑1 (MQ=255)
                                                           tggtgTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTaa                                                                                        >  2:1081138/1‑148 (MQ=255)
                                                            ggtgTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAAc                                                                                       <  2:726950/148‑1 (MQ=255)
                                                                 cTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTAc                                                                                  >  2:634376/1‑148 (MQ=255)
                                                                      gCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACt                                                                             >  2:2100487/1‑148 (MQ=255)
                                                                         tGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTa                                                                          >  1:2884688/1‑148 (MQ=255)
                                                                         tGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTa                                                                          >  2:1458236/1‑148 (MQ=255)
                                                                           ggCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGc                                                                        >  2:1865573/1‑148 (MQ=255)
                                                                            gCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTTACTTTAGCg                                                                       >  1:1002015/1‑148 (MQ=255)
                                                                            gCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCg                                                                       >  1:2030260/1‑148 (MQ=255)
                                                                             ccAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGt                                                                      <  1:3221646/148‑1 (MQ=255)
                                                                              cAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTc                                                                     >  2:2925149/1‑148 (MQ=255)
                                                                               aGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCa                                                                    <  1:3069179/148‑1 (MQ=255)
                                                                                ggTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCaa                                                                   <  1:2341916/148‑1 (MQ=255)
                                                                                   tAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTa                                                                <  1:3412215/148‑1 (MQ=255)
                                                                                   tAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTa                                                                <  1:2359906/148‑1 (MQ=255)
                                                                                     gCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAAt                                                              >  2:1547440/1‑148 (MQ=255)
                                                                                        gcgcCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTCGCGTCAATTAATCag                                                           >  1:933534/1‑148 (MQ=255)
                                                                                               gagCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACg                                                    <  1:1916156/148‑1 (MQ=255)
                                                                                               gagCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACg                                                    >  1:2345389/1‑148 (MQ=255)
                                                                                                agCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACgg                                                   >  2:2361862/1‑148 (MQ=255)
                                                                                                    aGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTaaa                                               >  2:2155214/1‑148 (MQ=255)
                                                                                                    aGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTaaa                                               >  1:745364/1‑148 (MQ=255)
                                                                                                     gggTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTaaaa                                              <  2:1197944/148‑1 (MQ=255)
                                                                                                     gggTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTaaaa                                              <  2:2964036/148‑1 (MQ=255)
                                                                                                     gggTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTaaaa                                              <  2:2807526/148‑1 (MQ=255)
                                                                                                      ggTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAAc                                             >  1:2779910/1‑148 (MQ=255)
                                                                                                      ggTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAAc                                             >  2:3208594/1‑148 (MQ=255)
                                                                                                       gTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACa                                            <  1:2149392/148‑1 (MQ=255)
                                                                                                       gTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACa                                            <  1:1092559/148‑1 (MQ=255)
                                                                                                        tAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACaa                                           <  2:3103046/148‑1 (MQ=255)
                                                                                                              caacaaACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAg                                     <  2:3482160/148‑1 (MQ=255)
                                                                                                                  aaacaGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTa                                 <  1:1659148/148‑1 (MQ=255)
                                                                                                                       gTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaaaaaa                            >  2:1378649/1‑148 (MQ=255)
                                                                                                                       gTGCAGTTAATTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTaaaaaa                            <  1:781137/148‑1 (MQ=255)
                                                                                                                        tGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAac                           >  2:2782689/1‑148 (MQ=255)
                                                                                                                         gCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAaca                          >  2:197235/1‑148 (MQ=255)
                                                                                                                           aGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACt                        <  1:1089704/148‑1 (MQ=255)
                                                                                                                            gTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACtt                       >  1:850890/1‑148 (MQ=255)
                                                                                                                                aTTTACGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTgc                   >  2:2403563/1‑148 (MQ=255)
                                                                                                                                  ttACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTgcgc                 <  1:1081138/148‑1 (MQ=255)
                                                                                                                                    aCGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCa               <  1:1042103/148‑1 (MQ=255)
                                                                                                                                    aCGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCa               <  1:1937606/148‑1 (MQ=255)
                                                                                                                                     cGCATATTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCaa              >  2:3075195/1‑148 (MQ=255)
                                                                                                                                       cATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAAtt            <  2:1104153/148‑1 (MQ=255)
                                                                                                                                        atatTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTa           >  1:1426179/1‑148 (MQ=255)
                                                                                                                                         tatTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTAt          <  1:260414/148‑1 (MQ=255)
                                                                                                                                          atTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATg         >  2:315571/1‑148 (MQ=255)
                                                                                                                                          atTCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATg         <  1:489973/148‑1 (MQ=255)
                                                                                                                                            tCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGga       >  1:3036335/1‑148 (MQ=255)
                                                                                                                                            tCTATATACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGga       >  2:27244/1‑148 (MQ=255)
                                                                                                                                             cTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGgag      >  2:2798506/1‑148 (MQ=255)
                                                                                                                                               atatCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCACCGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGaa    >  1:1595912/1‑148 (MQ=255)
                                                                                                                                               atatCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGaa    >  1:3144788/1‑148 (MQ=255)
                                                                                                                                               atatACTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGaa    >  1:2192636/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGGACCGAACTTGCCTTTGTGGTGCATCATCGGCTTCGCGTCAGCCGGTGCTGCGGTAGTGGTGTCTGCGGCATGGGCCAGGTTAGCCGCGCCAAGAGCCAGGGTAGAGGCAACAAACAGTGCAGTTAATTTACGCATATTCTATATCCTTCCTTTCAGTTATTTATTACGGCTTTCTTAAGTAGCGTGCCGTGTTGACGAGATTAACTTTACTGACTTTAGCGTCAATTAATCAGAGCAACGGTAAAACAATGAAAGTGTAAAAAACACTTTTGCGCCAATTATGGAGAAAA  >  CP009273/1819745‑1820037

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: