Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,219,858 G→T 6.8% V71V (GTC→GTA pepE ← (alpha)‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,219,8580GT 6.8% 824.2 / 21.7 220V71V (GTC→GTApepE(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (116/89);  new base T (7/8);  total (123/97)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 221 total aligned reads displayed.

CTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGC  >  CP009273/4219712‑4220005
                                                                                                                                                  |                                                                                                                                                   
  ccAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACa                                                                                                                                                  <  2:2821447/148‑1 (MQ=255)
   cAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACaa                                                                                                                                                 >  2:906152/1‑148 (MQ=255)
       cAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCt                                                                                                                                             <  1:367610/148‑1 (MQ=255)
        aaTATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTa                                                                                                                                            >  1:823030/1‑148 (MQ=255)
          tataCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATg                                                                                                                                          <  2:135175/148‑1 (MQ=255)
          tataCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATg                                                                                                                                          >  2:1228818/1‑148 (MQ=255)
             aCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAAt                                                                                                                                       >  2:2884944/1‑148 (MQ=255)
                  gcgcCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCgg                                                                                                                                  <  2:2550253/148‑1 (MQ=255)
                     ccACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGa                                                                                                                               >  2:2511596/1‑148 (MQ=255)
                      cACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGAc                                                                                                                              >  2:465812/1‑148 (MQ=255)
                         gTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGa                                                                                                                           <  1:2020410/148‑1 (MQ=255)
                         gTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGa                                                                                                                           <  1:2020493/148‑1 (MQ=255)
                          tttAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGaa                                                                                                                          >  1:1412753/1‑148 (MQ=255)
                           ttAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGaaa                                                                                                                         >  1:2057103/1‑148 (MQ=255)
                           ttAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGaaa                                                                                                                         >  1:3121146/1‑148 (MQ=255)
                              aCCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAacac                                                                                                                      <  2:139925/148‑1 (MQ=255)
                               ccACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACAcc                                                                                                                     <  2:102574/148‑1 (MQ=255)
                               ccACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACAcc                                                                                                                     <  1:2219985/148‑1 (MQ=255)
                               ccACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACAcc                                                                                                                     <  1:1640838/148‑1 (MQ=255)
                               ccACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACAcc                                                                                                                     >  2:472067/1‑148 (MQ=255)
                                cACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACAccc                                                                                                                    >  1:1968750/1‑148 (MQ=255)
                                cACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACAccc                                                                                                                    >  1:3042981/1‑148 (MQ=255)
                                 aCGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCa                                                                                                                   >  2:1170558/1‑148 (MQ=255)
                                 aCGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCa                                                                                                                   >  2:1129469/1‑148 (MQ=255)
                                     cAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCgg                                                                                                               >  2:2691980/1‑148 (MQ=255)
                                      aGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCACCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGa                                                                                                              <  2:1454487/148‑1 (MQ=255)
                                       gTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAg                                                                                                             >  1:2350029/1‑148 (MQ=255)
                                        tAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGc                                                                                                            <  2:3396316/148‑1 (MQ=255)
                                         aaTTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGGAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCg                                                                                                           <  2:675644/148‑1 (MQ=255)
                                           ttGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCgag                                                                                                         >  2:2394504/1‑148 (MQ=255)
                                            tGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCgaga                                                                                                        <  2:2190151/148‑1 (MQ=255)
                                            tGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCgaga                                                                                                        <  2:2364520/148‑1 (MQ=255)
                                              gTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAAc                                                                                                      <  2:1510151/148‑1 (MQ=255)
                                                   agcagcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAg                                                                                                 >  2:3424430/1‑148 (MQ=255)
                                                     cagcCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGcc                                                                                               >  2:2222885/1‑148 (MQ=255)
                                                       gcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGt                                                                                             <  1:2536872/148‑1 (MQ=255)
                                                              cTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCt                                                                                      <  2:3100969/148‑1 (MQ=255)
                                                                 gcgGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCtgtg                                                                                   <  1:2186952/148‑1 (MQ=255)
                                                                   ggCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTa                                                                                 <  2:2024023/148‑1 (MQ=255)
                                                                        tGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatcat                                                                            >  1:473996/1‑148 (MQ=255)
                                                                          tttCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATcc                                                                          <  2:1179086/148‑1 (MQ=255)
                                                                          tttCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATcc                                                                          >  1:469361/1‑148 (MQ=255)
                                                                          tttCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATcc                                                                          <  1:1133621/148‑1 (MQ=255)
                                                                          tttCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATcc                                                                          <  1:3533019/148‑1 (MQ=255)
                                                                           ttCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATccc                                                                         <  1:2103966/148‑1 (MQ=255)
                                                                            tCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCa                                                                        <  2:1968750/148‑1 (MQ=255)
                                                                             cAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAg                                                                       <  1:44808/148‑1 (MQ=255)
                                                                                cAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTc                                                                    >  2:1359005/1‑148 (MQ=255)
                                                                                  aCTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTg                                                                  >  1:532784/1‑148 (MQ=255)
                                                                                   cTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGc                                                                 <  2:3044308/148‑1 (MQ=255)
                                                                                    tGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCg                                                                >  2:2597832/1‑148 (MQ=255)
                                                                                     ggAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGt                                                               >  1:1378599/1‑148 (MQ=255)
                                                                                         aGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACg                                                           >  1:1299832/1‑148 (MQ=255)
                                                                                          gTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGc                                                          >  2:751985/1‑148 (MQ=255)
                                                                                          gTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGc                                                          >  2:1556652/1‑148 (MQ=255)
                                                                                           tATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGcc                                                         >  1:148053/1‑148 (MQ=255)
                                                                                            aTTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCa                                                        >  2:2161878/1‑148 (MQ=255)
                                                                                             ttCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAg                                                       >  2:924520/1‑148 (MQ=255)
                                                                                              tCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGc                                                      >  1:3375644/1‑148 (MQ=255)
                                                                                                ccgccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGa                                                    >  2:333747/1‑148 (MQ=255)
                                                                                                 cgccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGaa                                                   >  1:1056909/1‑148 (MQ=255)
                                                                                                 cgccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGaa                                                   <  1:2172326/148‑1 (MQ=255)
                                                                                                 cgccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGaa                                                   >  2:194347/1‑148 (MQ=255)
                                                                                                  gccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGaaa                                                  <  1:1503576/148‑1 (MQ=255)
                                                                                                  gccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGaaa                                                  <  1:465183/148‑1 (MQ=255)
                                                                                                     gccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAggg                                               <  2:52070/148‑1 (MQ=255)
                                                                                                     gccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAggg                                               <  2:2787565/148‑1 (MQ=255)
                                                                                                     gccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAggg                                               <  1:2703712/148‑1 (MQ=255)
                                                                                                         aCAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATaa                                           >  1:499019/1‑148 (MQ=255)
                                                                                                          cAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGGATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATaaa                                          <  1:1102499/148‑1 (MQ=255)
                                                                                                           aaTCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAac                                         <  2:737077/148‑1 (MQ=255)
                                                                                                            aTCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAaca                                        >  1:2418599/1‑148 (MQ=255)
                                                                                                            aTCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAaca                                        <  1:906159/148‑1 (MQ=255)
                                                                                                             tCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATAATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAacac                                       >  1:2395210/1‑148 (MQ=255)
                                                                                                                cGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGc                                    >  1:1890462/1‑148 (MQ=255)
                                                                                                                   tctcAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGa                                 >  1:2350400/1‑148 (MQ=255)
                                                                                                                   tctcAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGa                                 >  1:291000/1‑148 (MQ=255)
                                                                                                                   tctcAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGa                                 >  1:3388314/1‑148 (MQ=255)
                                                                                                                    ctcAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTGTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAg                                >  1:692651/1‑148 (MQ=255)
                                                                                                                       aGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcgg                             <  1:290026/148‑1 (MQ=255)
                                                                                                                         cATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcggcg                           >  1:510199/1‑148 (MQ=255)
                                                                                                                           ttttCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGAc                         >  1:606445/1‑148 (MQ=255)
                                                                                                                           ttttCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGAc                         >  1:1656872/1‑148 (MQ=255)
                                                                                                                           ttttCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCCCTGAGCGGCGAc                         >  1:321643/1‑148 (MQ=255)
                                                                                                                            tttCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGAcc                        <  2:294634/148‑1 (MQ=255)
                                                                                                                            tttCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGAcc                        >  1:432403/1‑148 (MQ=255)
                                                                                                                            tttCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAAAACCGCTGAGCGGCGAcc                        >  1:2108265/1‑148 (MQ=255)
                                                                                                                             ttCAATCGCGGCAACGGGATCGACAACGCTATGAATAACGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTGTTCGCAGTTTAATCATCCCAGGTCATCCGTACCCCAGCGAAACCGATAAACACCCCAGAGCCGCGAccc                       >  1:1088002/1‑148 (MQ=255)
                                                                                                                              tCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACAGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCt                      <  1:3164167/148‑1 (MQ=255)
                                                                                                                                 aTCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCa                   >  2:1843615/1‑148 (MQ=255)
                                                                                                                                  tCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCaa                  >  2:3156329/1‑148 (MQ=255)
                                                                                                                                     cgGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTg               <  1:96061/148‑1 (MQ=255)
                                                                                                                                       gCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGtt             >  1:1171077/1‑148 (MQ=255)
                                                                                                                                        cAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTc            >  1:2929989/1‑148 (MQ=255)
                                                                                                                                        cAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTc            >  1:3397563/1‑148 (MQ=255)
                                                                                                                                        cAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATAATCCCAGGTCTGCGTTACGCCCCCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCCACTGTTc            >  1:1785061/1‑148 (MQ=255)
                                                                                                                                             ggATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCaa       >  1:3250563/1‑148 (MQ=255)
                                                                                                                                             ggATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCaa       >  1:107205/1‑148 (MQ=255)
                                                                                                                                                tCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTa    >  2:2976079/1‑148 (MQ=255)
                                                                                                                                                  gACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGc  <  1:1353433/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGC  >  CP009273/4219712‑4220005

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: