Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,870,278 C→A 100% I113I (ATC→ATA iap → aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,870,2780CA95.8% 528.7 / 12.8 179I113I (ATC→ATAiapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base C (6/2);  new base A (84/87);  total (90/89)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 187 total aligned reads displayed.

TTCAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACG  >  CP009273/2870133‑2870422
                                                                                                                                                 |                                                                                                                                                
ttCAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATcataat                                                                                                                                                <  1:2245369/148‑1 (MQ=255)
ttCAGCAAATGGGTTATCGCAGTGATAGTCGGACATTTATTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTTTTTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATcataat                                                                                                                                                <  1:1942788/148‑1 (MQ=255)
   aGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTAt                                                                                                                                             >  1:2736187/1‑148 (MQ=255)
     cAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATgg                                                                                                                                           <  2:2965316/148‑1 (MQ=255)
     cAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATgg                                                                                                                                           <  2:1799609/148‑1 (MQ=255)
        aTGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGgcgc                                                                                                                                        >  2:2462851/1‑148 (MQ=255)
        aTGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGgcgc                                                                                                                                        >  1:2272703/1‑148 (MQ=255)
               aTCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGa                                                                                                                                 <  2:124536/148‑1 (MQ=255)
               aTCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGa                                                                                                                                 <  2:2547786/148‑1 (MQ=255)
                 cGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATACTTATGGCGCATCTGGATa                                                                                                                               >  2:2057174/1‑148 (MQ=255)
                   cAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACt                                                                                                                             <  2:570350/148‑1 (MQ=255)
                      tGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTAc                                                                                                                          >  2:1852994/1‑148 (MQ=255)
                       gATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACg                                                                                                                         >  1:1672344/1‑148 (MQ=255)
                       gATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACg                                                                                                                         >  2:2613618/1‑148 (MQ=255)
                         tatTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGcc                                                                                                                       <  2:491025/148‑1 (MQ=255)
                           ttCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGcccc                                                                                                                     >  2:3412546/1‑148 (MQ=255)
                           ttCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGcccc                                                                                                                     >  2:1647916/1‑148 (MQ=255)
                            tCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCg                                                                                                                    >  2:1685001/1‑148 (MQ=255)
                              ggACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCt                                                                                                                  <  2:3293762/148‑1 (MQ=255)
                               gACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTg                                                                                                                 <  1:765116/148‑1 (MQ=255)
                               gACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTg                                                                                                                 >  1:826534/1‑148 (MQ=255)
                                  aTTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGc                                                                                                              >  1:2928910/1‑148 (MQ=255)
                                   tttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                                                             >  2:2446604/1‑148 (MQ=255)
                                    ttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGa                                                                                                            <  2:3451559/148‑1 (MQ=255)
                                       aTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGc                                                                                                         >  2:3298515/1‑148 (MQ=255)
                                        tAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                                                        >  1:222776/1‑148 (MQ=255)
                                         aGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTg                                                                                                       <  2:1910243/148‑1 (MQ=255)
                                           tCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGAc                                                                                                     >  1:1148161/1‑148 (MQ=255)
                                           tCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGAc                                                                                                     >  1:1148167/1‑148 (MQ=255)
                                             ggTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGc                                                                                                   <  2:3532846/148‑1 (MQ=255)
                                              gTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGcc                                                                                                  >  1:2622784/1‑148 (MQ=255)
                                                atatTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGa                                                                                                <  2:753208/148‑1 (MQ=255)
                                                 tatTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGAt                                                                                               >  1:3551364/1‑148 (MQ=255)
                                                       taCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAt                                                                                         <  1:3336275/148‑1 (MQ=255)
                                                         ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAtct                                                                                       <  2:1610397/148‑1 (MQ=255)
                                                            cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcgg                                                                                    <  1:2018441/148‑1 (MQ=255)
                                                             ccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggc                                                                                   >  2:533447/1‑148 (MQ=255)
                                                             ccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggc                                                                                   >  2:2167585/1‑148 (MQ=255)
                                                              cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggcg                                                                                  <  1:1984851/148‑1 (MQ=255)
                                                                  aTAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                              >  2:2903017/1‑148 (MQ=255)
                                                                     aTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGAc                                                                           >  2:2148702/1‑148 (MQ=255)
                                                                     aTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGAc                                                                           >  2:985686/1‑148 (MQ=255)
                                                                      tCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                          >  1:1006741/1‑148 (MQ=255)
                                                                       cGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGt                                                                         <  1:1825102/148‑1 (MQ=255)
                                                                         tAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                       <  2:127077/148‑1 (MQ=255)
                                                                         tAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                       >  1:316084/1‑148 (MQ=255)
                                                                                tGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                <  2:3210463/148‑1 (MQ=255)
                                                                                tGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                <  1:2942207/148‑1 (MQ=255)
                                                                                  gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGACTg                                                              >  1:1039828/1‑148 (MQ=255)
                                                                                  gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                              >  2:2857326/1‑148 (MQ=255)
                                                                                  gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                              >  2:1109493/1‑148 (MQ=255)
                                                                                   cATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                                                             <  2:1830634/148‑1 (MQ=255)
                                                                                    aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGga                                                            <  1:2542401/148‑1 (MQ=255)
                                                                                     tAACGTGACGGGAAGTACGGTGATTGCCGCTCATGTAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgat                                                           >  2:2467805/1‑148 (MQ=255)
                                                                                      aaCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatg                                                          >  2:3299219/1‑148 (MQ=255)
                                                                                        cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                        >  2:3031893/1‑148 (MQ=255)
                                                                                        cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                        <  2:3328741/148‑1 (MQ=255)
                                                                                        cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                        <  2:3328799/148‑1 (MQ=255)
                                                                                           gACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAc                                                     >  1:2569193/1‑148 (MQ=255)
                                                                                            aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcg                                                    >  1:335450/1‑148 (MQ=255)
                                                                                            aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcg                                                    <  1:988376/148‑1 (MQ=255)
                                                                                                gAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccgc                                                >  1:3492886/1‑148 (MQ=255)
                                                                                                 aaGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCa                                               <  2:177347/148‑1 (MQ=255)
                                                                                                    tACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGt                                            >  1:1280798/1‑148 (MQ=255)
                                                                                                          gATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgt                                      <  1:2998686/148‑1 (MQ=255)
                                                                                                          gATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgt                                      <  1:2507617/148‑1 (MQ=255)
                                                                                                           aTTGCCGCTCATGAAGGAAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGGTGACGCCTATGCCAAGCTCGGCGGGCTGACGTTACCAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                     <  1:2823134/148‑1 (MQ=255)
                                                                                                             tGCCGCTCATGACGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTc                                   <  1:756720/148‑1 (MQ=255)
                                                                                                                cGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                <  2:3138610/148‑1 (MQ=255)
                                                                                                                 gCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGc                               <  2:3335184/148‑1 (MQ=255)
                                                                                                                     aTGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGa                           >  2:659889/1‑148 (MQ=255)
                                                                                                                     aTGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGa                           <  1:2206801/148‑1 (MQ=255)
                                                                                                                      tGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGaa                          >  1:2332162/1‑148 (MQ=255)
                                                                                                                      tGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGaa                          <  2:3054415/148‑1 (MQ=255)
                                                                                                                      tGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGaa                          <  1:2057174/148‑1 (MQ=255)
                                                                                                                       gAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAAt                         <  2:2348494/148‑1 (MQ=255)
                                                                                                                        aaGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAAtt                        <  2:1474142/148‑1 (MQ=255)
                                                                                                                          ggCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTgg                      >  1:1731533/1‑148 (MQ=255)
                                                                                                                          ggCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTgg                      >  2:1050415/1‑148 (MQ=255)
                                                                                                                           gCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGc                     <  1:1685001/148‑1 (MQ=255)
                                                                                                                           gCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGc                     >  1:3061998/1‑148 (MQ=255)
                                                                                                                              aaGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGa                  >  1:2433614/1‑148 (MQ=255)
                                                                                                                               aGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                 >  2:2511060/1‑148 (MQ=255)
                                                                                                                               aGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                 <  2:187781/148‑1 (MQ=255)
                                                                                                                                 cgccgcAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACg               >  1:150773/1‑148 (MQ=255)
                                                                                                                                 cgccgcAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACg               <  2:1633678/148‑1 (MQ=255)
                                                                                                                                  gccgcAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGc              <  2:230056/148‑1 (MQ=255)
                                                                                                                                  gccgcAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGc              <  1:3327181/148‑1 (MQ=255)
                                                                                                                                   ccgcAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGcc             <  2:417445/148‑1 (MQ=255)
                                                                                                                                    cgcAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCt            <  2:2422067/148‑1 (MQ=255)
                                                                                                                                      cagcagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGa          >  1:3541106/1‑148 (MQ=255)
                                                                                                                                       agcagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaa         >  1:3278201/1‑148 (MQ=255)
                                                                                                                                         cagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACTCCGCAGGTTTAGGTGTCATCCTGGAATTGGCAGAACGCCTGaaaa       <  2:1148161/148‑1 (MQ=255)
                                                                                                                                         cagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaa       <  2:1148163/148‑1 (MQ=255)
                                                                                                                                          agATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGACTGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaaa      >  2:1385786/1‑148 (MQ=255)
                                                                                                                                          agATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaaa      <  1:2418348/148‑1 (MQ=255)
                                                                                                                                           gATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAAt     <  2:1658801/148‑1 (MQ=255)
                                                                                                                                            atcataATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATa    >  1:1403002/1‑148 (MQ=255)
                                                                                                                                              cataatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACg  <  2:3551364/148‑1 (MQ=255)
                                                                                                                                              cataatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACg  <  2:1154456/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACG  >  CP009273/2870133‑2870422

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: