Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,644,550 C→A 24.5% pseudogene (216/765 nt) ydfE → Qin prophage; pseudogene;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,644,5500CA24.5% 586.1 / 136.3 248pseudogene (216/765 nt)ydfEQin prophage; pseudogene;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (98/89);  new base A (29/32);  total (127/121)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

AGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTC  >  CP009273/1644406‑1644697
                                                                                                                                                |                                                                                                                                                   
 gCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGAtt                                                                                                                                                 <  2:196988/148‑1 (MQ=255)
  ccACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTc                                                                                                                                                <  2:1411605/148‑1 (MQ=255)
    acaGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCAt                                                                                                                                              >  1:554355/1‑148 (MQ=255)
    acaGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCAt                                                                                                                                              <  1:473312/148‑1 (MQ=255)
    acaGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCAt                                                                                                                                              >  1:2081058/1‑148 (MQ=255)
      aGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATgg                                                                                                                                            <  1:29025/148‑1 (MQ=255)
       gCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATggg                                                                                                                                           <  1:132753/148‑1 (MQ=255)
       gCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATggg                                                                                                                                           >  1:2970710/1‑148 (MQ=255)
           ttGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGc                                                                                                                                       <  1:430247/148‑1 (MQ=255)
             gAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTg                                                                                                                                     <  2:3141304/148‑1 (MQ=255)
                gCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAat                                                                                                                                  >  2:3188256/1‑148 (MQ=255)
                gCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAat                                                                                                                                  <  1:446994/148‑1 (MQ=255)
                 cTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAata                                                                                                                                 <  2:2451174/148‑1 (MQ=255)
                  tGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAatat                                                                                                                                <  1:2074703/148‑1 (MQ=255)
                   ggCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAatata                                                                                                                               >  1:3130717/1‑148 (MQ=255)
                   ggCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAatata                                                                                                                               >  1:338829/1‑148 (MQ=255)
                   ggCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAatata                                                                                                                               <  1:2232234/148‑1 (MQ=255)
                   ggCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAatata                                                                                                                               <  2:1146036/148‑1 (MQ=255)
                    gCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAatatat                                                                                                                              >  2:2014201/1‑148 (MQ=255)
                    gCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAatatat                                                                                                                              >  1:1370184/1‑148 (MQ=255)
                       acgaACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCa                                                                                                                           <  1:179030/148‑1 (MQ=255)
                          aaCGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGa                                                                                                                        <  2:1880963/148‑1 (MQ=255)
                                gTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCg                                                                                                                  <  1:706018/148‑1 (MQ=255)
                                  gCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCt                                                                                                                <  1:3039980/148‑1 (MQ=255)
                                      gAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGa                                                                                                            <  2:1332174/148‑1 (MQ=255)
                                         cAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATga                                                                                                         <  1:1881752/148‑1 (MQ=255)
                                            aaTTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATgagga                                                                                                      <  2:337888/148‑1 (MQ=255)
                                              tttCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGt                                                                                                    <  1:1100529/148‑1 (MQ=255)
                                              tttCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGt                                                                                                    >  2:111847/1‑148 (MQ=255)
                                               ttCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTc                                                                                                   <  1:2090771/148‑1 (MQ=255)
                                               ttCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTc                                                                                                   >  2:2564858/1‑148 (MQ=255)
                                               ttCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTc                                                                                                   <  2:2717425/148‑1 (MQ=255)
                                                 cACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTg                                                                                                 <  2:3540523/148‑1 (MQ=255)
                                                  aCTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGc                                                                                                <  2:2596897/148‑1 (MQ=255)
                                                    tGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGctt                                                                                              >  2:2346443/1‑148 (MQ=255)
                                                      gCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttct                                                                                            <  1:1298283/148‑1 (MQ=255)
                                                      gCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttct                                                                                            >  1:142062/1‑148 (MQ=255)
                                                      gCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTACAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCGGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttct                                                                                            <  1:2314748/148‑1 (MQ=255)
                                                       cAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttctt                                                                                           <  2:982669/148‑1 (MQ=255)
                                                       cAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttctt                                                                                           >  2:1436298/1‑148 (MQ=255)
                                                        aaTTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCttt                                                                                          <  1:2380836/148‑1 (MQ=255)
                                                          ttCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTg                                                                                        >  2:929168/1‑148 (MQ=255)
                                                          ttCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTg                                                                                        >  1:1260419/1‑148 (MQ=255)
                                                              tGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTccc                                                                                    >  1:2746533/1‑148 (MQ=255)
                                                                gAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCgg                                                                                  >  2:1209587/1‑148 (MQ=255)
                                                                 aaCATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGt                                                                                 <  1:845366/148‑1 (MQ=255)
                                                                  aCATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGtt                                                                                >  1:2950972/1‑148 (MQ=255)
                                                                   cATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTg                                                                               <  1:3245635/148‑1 (MQ=255)
                                                                   cATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTg                                                                               >  2:2321786/1‑148 (MQ=255)
                                                                   cATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTg                                                                               >  1:2586636/1‑148 (MQ=255)
                                                                   cATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTg                                                                               >  1:2283981/1‑148 (MQ=255)
                                                                      tATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGGCCGCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTACCGGTTgaga                                                                            >  2:1669435/1‑148 (MQ=255)
                                                                      tATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTgaga                                                                            >  2:3432480/1‑148 (MQ=255)
                                                                       aTCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGAc                                                                           >  1:911255/1‑148 (MQ=255)
                                                                         ccGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTg                                                                         >  2:2826798/1‑148 (MQ=255)
                                                                          cGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGa                                                                        >  2:3139757/1‑148 (MQ=255)
                                                                          cGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGa                                                                        >  2:1958425/1‑148 (MQ=255)
                                                                           ggCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGaa                                                                       >  2:2357971/1‑148 (MQ=255)
                                                                                ctCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatccga                                                                  <  2:122865/148‑1 (MQ=255)
                                                                                ctCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatccga                                                                  >  2:910247/1‑148 (MQ=255)
                                                                                  cAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatccgatc                                                                <  1:996928/148‑1 (MQ=255)
                                                                                       cTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATg                                                           >  2:2287687/1‑148 (MQ=255)
                                                                                       cTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATg                                                           <  1:1060956/148‑1 (MQ=255)
                                                                                           gcTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACg                                                       <  1:3513783/148‑1 (MQ=255)
                                                                                             ttATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTc                                                     <  1:748246/148‑1 (MQ=255)
                                                                                              tataAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCa                                                    >  2:105242/1‑148 (MQ=255)
                                                                                               ataAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCAc                                                   <  1:1605009/148‑1 (MQ=255)
                                                                                               ataAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCAc                                                   <  1:1604897/148‑1 (MQ=255)
                                                                                                   aTTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTCAATCCGATCCGATGAACGACACtttt                                               >  1:3360359/1‑148 (MQ=255)
                                                                                                    tttATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTg                                              >  2:1736663/1‑148 (MQ=255)
                                                                                                    tttATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTg                                              >  1:2037635/1‑148 (MQ=255)
                                                                                                     ttATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGa                                             >  1:3195776/1‑148 (MQ=255)
                                                                                                        ttGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACaa                                          >  2:3240848/1‑148 (MQ=255)
                                                                                                        ttGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACaa                                          >  2:1618571/1‑148 (MQ=255)
                                                                                                              gCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGc                                    <  1:2626025/148‑1 (MQ=255)
                                                                                                              gCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGc                                    >  2:841006/1‑148 (MQ=255)
                                                                                                              gCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGc                                    >  2:2627081/1‑148 (MQ=255)
                                                                                                                gAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGccc                                  <  1:254992/148‑1 (MQ=255)
                                                                                                                 aGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGcccc                                 >  1:782617/1‑148 (MQ=255)
                                                                                                                      aaaCCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAg                            >  1:3533260/1‑148 (MQ=255)
                                                                                                                      aaaCCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAg                            >  1:820432/1‑148 (MQ=255)
                                                                                                                      aaaCCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAg                            >  2:2179441/1‑148 (MQ=255)
                                                                                                                        aCCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTa                          >  1:80294/1‑148 (MQ=255)
                                                                                                                        aCCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTa                          <  2:669377/148‑1 (MQ=255)
                                                                                                                           tGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAg                       >  1:1786774/1‑148 (MQ=255)
                                                                                                                           tGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAg                       >  2:2336327/1‑148 (MQ=255)
                                                                                                                           tGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAg                       >  1:2508792/1‑148 (MQ=255)
                                                                                                                            ggCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGa                      <  1:111847/148‑1 (MQ=255)
                                                                                                                             gCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGaa                     >  1:1271364/1‑148 (MQ=255)
                                                                                                                               aTACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACg                   <  2:1288875/148‑1 (MQ=255)
                                                                                                                                  ccccGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTg                >  2:553536/1‑148 (MQ=255)
                                                                                                                                  ccccGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTg                >  2:2498539/1‑148 (MQ=255)
                                                                                                                                     cGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCa             >  2:3420780/1‑148 (MQ=255)
                                                                                                                                     cGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCa             <  1:2348767/148‑1 (MQ=255)
                                                                                                                                     cGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCa             <  2:2089643/148‑1 (MQ=255)
                                                                                                                                     cGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCa             <  1:1219922/148‑1 (MQ=255)
                                                                                                                                     cGCCCTGCCCTCGATACATGGGATGCTGAATATAGGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGGTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCa             <  1:1219914/148‑1 (MQ=255)
                                                                                                                                       ccctgccctCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATg           >  2:699185/1‑148 (MQ=255)
                                                                                                                                            ccctAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTaa      <  2:2970710/148‑1 (MQ=255)
                                                                                                                                                cGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTc  <  2:478166/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTC  >  CP009273/1644406‑1644697

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: