Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,284,857 G→T 100% D53Y (GAT→TAT)  rssA → putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,284,8570GT98.6% 657.5 / ‑0.8 209D53Y (GAT→TAT) rssAputative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base T (113/93);  total (113/96)
Fisher's exact test for biased strand distribution p-value = 9.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 216 total aligned reads displayed.

AAAATAGGGCTGGCGCTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCGATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAG  >  CP009273/1284713‑1285001
                                                                                                                                                |                                                                                                                                                
 aaaTAGGGCTGGCGCTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGctatcg                                                                                                                                              <  1:372155/148‑1 (MQ=255)
 aaaTAGGGCTGGCGCTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGctatcg                                                                                                                                              <  1:1368840/148‑1 (MQ=255)
   aTAGGGCTGGCGCTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGctatcgat                                                                                                                                            >  2:1683516/1‑148 (MQ=255)
      gggCTGGCGCTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTAt                                                                                                                                         <  2:1786340/148‑1 (MQ=255)
       ggCTGGCGCTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATc                                                                                                                                        <  2:1853541/148‑1 (MQ=255)
       ggCTGGCGCTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATc                                                                                                                                        >  2:362896/1‑148 (MQ=255)
        gCTGGCGCTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCt                                                                                                                                       >  2:2341386/1‑148 (MQ=255)
         cTGGCGCTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTg                                                                                                                                      <  2:2407937/148‑1 (MQ=255)
           ggCGCTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTgcg                                                                                                                                    <  1:405075/148‑1 (MQ=255)
            gcgcTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTgcgc                                                                                                                                   >  1:1667899/1‑148 (MQ=255)
                               cgAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACctc                                                                                                                >  1:2861386/1‑148 (MQ=255)
                               cgAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACctc                                                                                                                <  1:1805086/148‑1 (MQ=255)
                                gagaGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACctct                                                                                                               <  1:3318579/148‑1 (MQ=255)
                                gagaGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGAGATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACctct                                                                                                               <  2:840333/148‑1 (MQ=255)
                                  gagGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCttt                                                                                                             >  1:2243882/1‑148 (MQ=255)
                                   agGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTc                                                                                                            >  1:489168/1‑148 (MQ=255)
                                    ggTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCa                                                                                                           <  2:3009515/148‑1 (MQ=255)
                                       tGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGtt                                                                                                        >  2:935178/1‑148 (MQ=255)
                                       tGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGtt                                                                                                        <  2:1580380/148‑1 (MQ=255)
                                       tGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGtt                                                                                                        >  1:3508893/1‑148 (MQ=255)
                                        ggTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTa                                                                                                       <  1:2594130/148‑1 (MQ=255)
                                        ggTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTa                                                                                                       <  2:2265150/148‑1 (MQ=255)
                                          tCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTAtt                                                                                                     >  2:316825/1‑148 (MQ=255)
                                           cGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTg                                                                                                    <  1:2440122/148‑1 (MQ=255)
                                               tatTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGAt                                                                                                >  2:169634/1‑148 (MQ=255)
                                               tatTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGAt                                                                                                >  2:798893/1‑148 (MQ=255)
                                                 ttGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGt                                                                                              >  1:558274/1‑148 (MQ=255)
                                                 ttGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGt                                                                                              <  1:692622/148‑1 (MQ=255)
                                                  tGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTTCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGAATATCTGCGCTGGAAGATTGGGTGACCTCTTTAAGTTATTGGGATGtt                                                                                             >  1:2435682/1‑148 (MQ=255)
                                                  tGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTTGGATGtt                                                                                             >  2:1622541/1‑148 (MQ=255)
                                                  tGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGtt                                                                                             <  1:2737833/148‑1 (MQ=255)
                                                     cGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTa                                                                                          <  1:2913080/148‑1 (MQ=255)
                                                      gTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTAc                                                                                         >  2:1450968/1‑148 (MQ=255)
                                                         attaATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGcc                                                                                      >  2:2949048/1‑148 (MQ=255)
                                                          ttaATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCt                                                                                     <  2:2751115/148‑1 (MQ=255)
                                                          ttaATGAGCTAAAAAAAGTGGGTATTGAAATTGATATCGGTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCt                                                                                     <  2:2131171/148‑1 (MQ=255)
                                                           taATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTg                                                                                    >  1:706550/1‑148 (MQ=255)
                                                           taATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTg                                                                                    <  2:1760480/148‑1 (MQ=255)
                                                           taATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTg                                                                                    >  2:2769538/1‑148 (MQ=255)
                                                           taATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTg                                                                                    >  1:1036322/1‑148 (MQ=255)
                                                           taATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTg                                                                                    >  1:1081921/1‑148 (MQ=255)
                                                            aaTGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGa                                                                                   <  1:2137293/148‑1 (MQ=255)
                                                            aaTGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGa                                                                                   <  1:321009/148‑1 (MQ=255)
                                                             aTGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGAt                                                                                  >  1:2207489/1‑148 (MQ=255)
                                                             aTGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGAt                                                                                  >  1:118610/1‑148 (MQ=255)
                                                             aTGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGAt                                                                                  <  1:2990281/148‑1 (MQ=255)
                                                             aTGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGAt                                                                                  >  1:1170890/1‑148 (MQ=255)
                                                               gcgcTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATgg                                                                                <  1:1450484/148‑1 (MQ=255)
                                                                cgcTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGa                                                                               <  2:543241/148‑1 (MQ=255)
                                                                    aaaaaaaGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGAtctc                                                                           >  2:2202071/1‑148 (MQ=255)
                                                                    aaaaaaaGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGAtctc                                                                           >  2:434951/1‑148 (MQ=255)
                                                                        aaaGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCt                                                                       <  1:2823239/148‑1 (MQ=255)
                                                                            tgtgTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCa                                                                   >  1:3368138/4‑148 (MQ=255)
                                                                            tGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCa                                                                   >  1:1426438/1‑148 (MQ=255)
                                                                            tGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTAAGCCTGATGGATCTCTCCTGGCa                                                                   >  1:465679/1‑148 (MQ=255)
                                                                             gggTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAg                                                                  >  1:2229303/1‑148 (MQ=255)
                                                                              ggTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAgc                                                                 >  2:1310564/1‑148 (MQ=255)
                                                                              ggTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAgc                                                                 >  1:3246650/1‑148 (MQ=255)
                                                                                     aaaTTGATATCGGTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTgg                                                          <  1:816166/148‑1 (MQ=255)
                                                                                      aaTTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTggg                                                         >  1:2117061/1‑148 (MQ=255)
                                                                                        ttGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGtt                                                       >  2:696587/1‑148 (MQ=255)
                                                                                        ttGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGtt                                                       >  2:1065009/1‑148 (MQ=255)
                                                                                         tGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTa                                                      <  1:1310564/148‑1 (MQ=255)
                                                                                          gATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTAc                                                     >  1:2765890/1‑148 (MQ=255)
                                                                                          gATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTAc                                                     <  2:487732/148‑1 (MQ=255)
                                                                                          gATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTAc                                                     >  1:2934094/1‑148 (MQ=255)
                                                                                             atCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTgc                                                  >  1:2693849/1‑148 (MQ=255)
                                                                                             atCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTgc                                                  >  1:3153881/1‑148 (MQ=255)
                                                                                             atCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTgc                                                  >  2:1379042/1‑148 (MQ=255)
                                                                                               cGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTgcgc                                                >  1:1162597/1‑148 (MQ=255)
                                                                                               cGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTgcgc                                                >  1:3485290/1‑148 (MQ=255)
                                                                                                 ttGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCgcgg                                              >  2:518578/1‑148 (MQ=255)
                                                                                                    cAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGa                                           >  2:2657774/1‑148 (MQ=255)
                                                                                                      ggATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGc                                         >  2:3026065/1‑148 (MQ=255)
                                                                                                        aTGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCgt                                       >  2:3484833/1‑148 (MQ=255)
                                                                                                          gTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCgtgt                                     <  2:1852381/148‑1 (MQ=255)
                                                                                                          gTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCgtgt                                     >  2:3226897/1‑148 (MQ=255)
                                                                                                           ttCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTc                                    <  1:698316/148‑1 (MQ=255)
                                                                                                               aTTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTtca                                >  1:2425360/1‑148 (MQ=255)
                                                                                                               aTTGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTtca                                >  1:2425566/1‑148 (MQ=255)
                                                                                                                ttGGTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACTCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTtcaa                               >  2:191253/1‑148 (MQ=255)
                                                                                                                   gTTCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTtcaatca                            <  2:2826779/148‑1 (MQ=255)
                                                                                                                    ttCGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGa                                           >  2:757036/1‑132 (MQ=255)
                                                                                                                      cGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATa                         <  2:1488295/148‑1 (MQ=255)
                                                                                                                      cGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATa                         <  1:1125335/148‑1 (MQ=255)
                                                                                                                      cGCTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATa                         <  1:1683516/148‑1 (MQ=255)
                                                                                                                        cTGGTGGGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATAcc                       >  1:967000/1‑148 (MQ=255)
                                                                                                                          ggtggGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACcgc                     >  1:2747533/1‑148 (MQ=255)
                                                                                                                          ggtggGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACcgc                     >  1:545824/1‑148 (MQ=255)
                                                                                                                            tggGCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGa                   <  2:1176552/148‑1 (MQ=255)
                                                                                                                             gggCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGaa                  <  1:3229957/148‑1 (MQ=255)
                                                                                                                              ggCGCTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGaaa                 <  2:1023158/148‑1 (MQ=255)
                                                                                                                                 gcTGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAaataa              <  1:1011586/148‑1 (MQ=255)
                                                                                                                                   tGCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATg            <  2:2016175/148‑1 (MQ=255)
                                                                                                                                    gCCTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGc           <  2:1882521/148‑1 (MQ=255)
                                                                                                                                     ccTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGcc          <  2:3368042/148‑1 (MQ=255)
                                                                                                                                     ccTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGcc          <  2:172536/148‑1 (MQ=255)
                                                                                                                                     ccTATGCATGCTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGcc          >  1:864782/1‑148 (MQ=255)
                                                                                                                                          gcatgcTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGaaa     >  2:713133/1‑148 (MQ=255)
                                                                                                                                           catgcTATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAAc    <  1:798893/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AAAATAGGGCTGGCGCTGGGATCTGGCGCGGCGAGAGGTTGGTCGCATATTGGCGTTATTAATGCGCTAAAAAAAGTGGGTATTGAAATTGATATCGTTGCAGGATGTTCAATTGGTTCGCTGGTGGGCGCTGCCTATGCATGCGATCGATTATCTGCGCTGGAAGATTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAG  >  CP009273/1284713‑1285001

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: