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.6% 616.0 / 18.1 206I113I (ATC→ATAiapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base C (9/0);  new base A (99/98);  total (108/98)
Fisher's exact test for biased strand distribution p-value = 3.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01
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 219 total aligned reads displayed.

TTCAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACG  >  CP009273/2870133‑2870422
                                                                                                                                                 |                                                                                                                                                
ttCAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATcataat                                                                                                                                                <  2:1427866/148‑1 (MQ=255)
     cAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATgg                                                                                                                                           >  2:2671687/1‑148 (MQ=255)
         tGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCa                                                                                                                                       >  1:724380/1‑148 (MQ=255)
         tGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCa                                                                                                                                       >  1:2507080/1‑148 (MQ=255)
           ggTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGTCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATc                                                                                                                                     >  1:2475428/1‑148 (MQ=255)
            gTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCt                                                                                                                                    >  2:511708/1‑148 (MQ=255)
             ttATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTg                                                                                                                                   >  1:2543419/1‑148 (MQ=255)
             ttATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTg                                                                                                                                   >  1:1941414/1‑148 (MQ=255)
             ttATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTg                                                                                                                                   <  2:127223/148‑1 (MQ=255)
              tATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTgg                                                                                                                                  <  2:1340643/148‑1 (MQ=255)
               aTCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGa                                                                                                                                 <  2:1866826/148‑1 (MQ=255)
               aTCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGa                                                                                                                                 <  1:1205658/148‑1 (MQ=255)
               aTCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGa                                                                                                                                 <  2:772117/148‑1 (MQ=255)
                tCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGAt                                                                                                                                <  2:1012352/148‑1 (MQ=255)
                   cAGTGATATTCGGAAATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACt                                                                                                                             <  1:1044757/148‑1 (MQ=255)
                    aGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACtt                                                                                                                            <  1:1293903/148‑1 (MQ=255)
                     gTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGCTCATAATTATGGCGCATCTGGATCCTTa                                                                                                                           >  2:2359534/1‑148 (MQ=255)
                     gTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTa                                                                                                                           <  2:2913664/148‑1 (MQ=255)
                       gATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACg                                                                                                                         >  1:377916/1‑148 (MQ=255)
                       gATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACg                                                                                                                         >  1:1193744/1‑148 (MQ=255)
                        atatTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGc                                                                                                                        >  2:2999804/1‑148 (MQ=255)
                        atatTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGc                                                                                                                        >  2:3000110/1‑148 (MQ=255)
                            tCGGACATTTAATTGTCGGTATTTTTATCCCGCCCGCTATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTAATGATGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCg                                                                                                                    <  1:2145486/148‑1 (MQ=255)
                            tCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCg                                                                                                                    <  1:2276883/148‑1 (MQ=255)
                            tCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCg                                                                                                                    >  2:839999/1‑148 (MQ=255)
                            tCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCg                                                                                                                    >  1:3281320/1‑148 (MQ=255)
                             cGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGc                                                                                                                   >  2:526915/1‑148 (MQ=255)
                               gACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTg                                                                                                                 <  1:2950487/148‑1 (MQ=255)
                                aCATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGa                                                                                                                >  1:1064328/1‑148 (MQ=255)
                                  aTTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGc                                                                                                              >  2:3447638/1‑148 (MQ=255)
                                   tttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                                                             >  2:162787/1‑148 (MQ=255)
                                      aaTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATg                                                                                                          >  2:676148/1‑148 (MQ=255)
                                      aaTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATg                                                                                                          <  2:3035795/148‑1 (MQ=255)
                                      aaTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATg                                                                                                          <  2:1072951/148‑1 (MQ=255)
                                        tAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                                                        >  1:1371055/1‑148 (MQ=255)
                                          gTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGa                                                                                                      >  2:113239/1‑148 (MQ=255)
                                           tCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGAc                                                                                                     >  2:1166493/1‑148 (MQ=255)
                                              gTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGcc                                                                                                  >  1:2818151/1‑148 (MQ=255)
                                               tataTTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCg                                                                                                 <  1:3387619/148‑1 (MQ=255)
                                               tataTTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCg                                                                                                 <  1:2241361/148‑1 (MQ=255)
                                               tataTTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCg                                                                                                 <  2:2607872/148‑1 (MQ=255)
                                                   tttATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGc                                                                                             <  1:108672/148‑1 (MQ=255)
                                                    ttATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGcc                                                                                            >  2:1931847/1‑148 (MQ=255)
                                                    ttATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGcc                                                                                            >  2:604374/1‑148 (MQ=255)
                                                     tataCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCa                                                                                           >  1:2328519/1‑148 (MQ=255)
                                                      ataCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCaa                                                                                          >  2:1297935/1‑148 (MQ=255)
                                                       taCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAt                                                                                         >  1:3523448/1‑148 (MQ=255)
                                                                  aTAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                              >  2:2282523/1‑148 (MQ=255)
                                                                  aTAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                              >  2:2451966/1‑148 (MQ=255)
                                                                    aaTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGa                                                                            >  1:384773/1‑148 (MQ=255)
                                                                      tCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                          <  2:294844/148‑1 (MQ=255)
                                                                          aaGAGCTGGCATAACGTGACGGGAAGTACGGTTATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTAc                                                                      >  1:2773135/1‑148 (MQ=255)
                                                                           agagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                     <  2:348724/148‑1 (MQ=255)
                                                                           agagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                     <  2:2064526/148‑1 (MQ=255)
                                                                           agagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                     <  2:2405786/148‑1 (MQ=255)
                                                                              gCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                                                                  <  2:3074065/148‑1 (MQ=255)
                                                                              gCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACGTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                                                                  <  2:1302708/148‑1 (MQ=255)
                                                                               cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                 <  2:1674413/148‑1 (MQ=255)
                                                                               cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                 <  1:868685/148‑1 (MQ=255)
                                                                               cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                 <  2:2400588/148‑1 (MQ=255)
                                                                                 ggCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAAt                                                               >  2:93643/1‑148 (MQ=255)
                                                                                   cATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                                                             <  2:1941414/148‑1 (MQ=255)
                                                                                   cATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                                                             >  2:1198467/1‑148 (MQ=255)
                                                                                    aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGga                                                            <  1:3052963/148‑1 (MQ=255)
                                                                                    aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGga                                                            <  2:3173598/148‑1 (MQ=255)
                                                                                     tAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgat                                                           >  2:1617776/1‑148 (MQ=255)
                                                                                       aCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatga                                                         <  2:1471585/148‑1 (MQ=255)
                                                                                       aCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatga                                                         >  1:2146675/1‑148 (MQ=255)
                                                                                       aCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatga                                                         <  2:3415272/148‑1 (MQ=255)
                                                                                             cGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgc                                                   >  1:839038/1‑148 (MQ=255)
                                                                                             cGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgc                                                   <  1:1489433/148‑1 (MQ=255)
                                                                                               ggAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccg                                                 <  2:1165058/148‑1 (MQ=255)
                                                                                                   gTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAgg                                             >  2:833984/1‑148 (MQ=255)
                                                                                                   gTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAgg                                             <  2:131839/148‑1 (MQ=255)
                                                                                                      cGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGttt                                          <  2:2806441/148‑1 (MQ=255)
                                                                                                         tGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAgg                                       >  2:2189760/1‑148 (MQ=255)
                                                                                                         tGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAgg                                       >  1:1112643/1‑148 (MQ=255)
                                                                                                         tGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAgg                                       >  2:2898342/1‑148 (MQ=255)
                                                                                                         tGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAgg                                       >  1:3553668/1‑148 (MQ=255)
                                                                                                         tGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAgg                                       <  2:3061015/148‑1 (MQ=255)
                                                                                                            tagCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtgt                                    >  1:2500039/3‑148 (MQ=255)
                                                                                                                cGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                <  2:618653/148‑1 (MQ=255)
                                                                                                                 gCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGc                               <  2:2853649/148‑1 (MQ=255)
                                                                                                                  cTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCt                              >  2:1084056/1‑148 (MQ=255)
                                                                                                                      tGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGaa                          <  2:2534449/148‑1 (MQ=255)
                                                                                                                       gAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAAt                         <  2:420691/148‑1 (MQ=255)
                                                                                                                          ggCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTgg                      <  2:839038/148‑1 (MQ=255)
                                                                                                                           gCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGc                     >  1:197687/1‑148 (MQ=255)
                                                                                                                           gCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGc                     <  2:45400/148‑1 (MQ=255)
                                                                                                                           gCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGc                     <  1:1168998/148‑1 (MQ=255)
                                                                                                                            cAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCa                    >  2:1104595/1‑148 (MQ=255)
                                                                                                                              aaGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGa                  <  2:2791185/148‑1 (MQ=255)
                                                                                                                                gcgcCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAAc                >  2:334346/1‑148 (MQ=255)
                                                                                                                                  gccgcAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGc              <  2:890653/148‑1 (MQ=255)
                                                                                                                                  gccgcAGCAGATCATAATTAGGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGc              <  1:1867748/148‑1 (MQ=255)
                                                                                                                                         cagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaa       >  1:2414978/1‑148 (MQ=255)
                                                                                                                                         cagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaa       <  1:890698/148‑1 (MQ=255)
                                                                                                                                         cagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaa       <  1:2671687/148‑1 (MQ=255)
                                                                                                                                          agATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaaa      <  1:113239/148‑1 (MQ=255)
                                                                                                                                            atcataATTAGGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATa    <  1:608308/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: