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,2780CA98.4% 639.3 / 0.9 199I113I (ATC→ATAiapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (100/95);  total (102/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 205 total aligned reads displayed.

CAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATAC  >  CP009273/2870135‑2870421
                                                                                                                                               |                                                                                                                                               
cAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTa                                                                                                                                             <  1:1176891/148‑1 (MQ=255)
cAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTa                                                                                                                                             <  2:2622646/148‑1 (MQ=255)
  gCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATg                                                                                                                                           >  1:1721535/1‑148 (MQ=255)
      aTGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGgcgc                                                                                                                                       >  1:1026020/1‑148 (MQ=255)
        gggTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCAt                                                                                                                                     <  2:2120372/148‑1 (MQ=255)
        gggTTATCGCAGTGATATTCGGACATGTAATAGTCGGTATATTTATACCGCCCGCGAAAATCGTAAGAGCTTGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGAGCCGCAGCAGTTCATAATTATGGCGCAt                                                                                                                                     >  2:949494/1‑148 (MQ=255)
         ggTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATc                                                                                                                                    <  2:211135/148‑1 (MQ=255)
            tATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTgg                                                                                                                                 >  1:170362/1‑148 (MQ=255)
            tATCGCAGTGATATTCGGACATTTAATAGTCGGGATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTgg                                                                                                                                 >  1:1630989/1‑148 (MQ=255)
             aTCGCAGTTATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGAGGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGa                                                                                                                                <  1:1952868/148‑1 (MQ=255)
             aTCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGa                                                                                                                                <  2:3367663/148‑1 (MQ=255)
             aTCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGa                                                                                                                                <  1:538433/148‑1 (MQ=255)
              tCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGTTAATCGTAAGAGCTGGCATAACTTGTCGGGAAGTACGGTGATTGCAGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGAt                                                                                                                               <  1:294511/148‑1 (MQ=255)
                     gATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACg                                                                                                                        >  2:2384176/1‑148 (MQ=255)
                      atatTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGc                                                                                                                       >  1:2289090/1‑148 (MQ=255)
                        atTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGccc                                                                                                                     >  1:2779647/1‑148 (MQ=255)
                        atTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGccc                                                                                                                     >  1:1491013/1‑148 (MQ=255)
                        atTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGccc                                                                                                                     >  1:1227731/1‑148 (MQ=255)
                          tCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCg                                                                                                                   >  1:1235552/1‑148 (MQ=255)
                          tCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCg                                                                                                                   <  1:1468132/148‑1 (MQ=255)
                               cATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAg                                                                                                              >  1:1909886/1‑148 (MQ=255)
                                 tttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                                                            >  2:276884/1‑148 (MQ=255)
                                 tttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                                                            >  2:1946708/1‑148 (MQ=255)
                                 tttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                                                            <  2:3088743/148‑1 (MQ=255)
                                 tttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                                                            >  2:879254/1‑148 (MQ=255)
                                  ttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGa                                                                                                           <  2:2337503/148‑1 (MQ=255)
                                  ttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGa                                                                                                           <  2:1192909/148‑1 (MQ=255)
                                  ttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGa                                                                                                           >  2:1025216/1‑148 (MQ=255)
                                    aaTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATg                                                                                                         <  1:3588019/148‑1 (MQ=255)
                                        gTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGa                                                                                                     <  2:3470359/148‑1 (MQ=255)
                                         tCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGAc                                                                                                    >  1:1407119/1‑148 (MQ=255)
                                            gTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGcc                                                                                                 >  2:2947464/1‑148 (MQ=255)
                                             tataTTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCg                                                                                                >  1:3496143/1‑148 (MQ=255)
                                              atatTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGa                                                                                               >  1:972267/1‑148 (MQ=255)
                                              atatTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGa                                                                                               >  2:3719957/1‑148 (MQ=255)
                                                 tttATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGc                                                                                            <  1:2750577/148‑1 (MQ=255)
                                                      aCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAtc                                                                                       >  2:3171974/1‑148 (MQ=255)
                                                        cgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATCCTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAtctc                                                                                     >  2:1364153/1‑148 (MQ=255)
                                                         gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcg                                                                                    <  1:3270207/148‑1 (MQ=255)
                                                         gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcg                                                                                    <  1:2384176/148‑1 (MQ=255)
                                                            cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggcg                                                                                 >  2:284371/1‑148 (MQ=255)
                                                            cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggcg                                                                                 >  2:2162432/1‑148 (MQ=255)
                                                              cgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCggg                                                                               >  2:1129284/1‑148 (MQ=255)
                                                                aTAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                             >  1:938147/1‑148 (MQ=255)
                                                                     cGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGt                                                                        <  2:1102879/148‑1 (MQ=255)
                                                                         agagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                    <  1:3754561/148‑1 (MQ=255)
                                                                             cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                <  2:1902818/148‑1 (MQ=255)
                                                                             cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                <  2:3202564/148‑1 (MQ=255)
                                                                             cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                <  1:2806837/148‑1 (MQ=255)
                                                                             cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                <  1:459321/148‑1 (MQ=255)
                                                                              tGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                               <  1:2855122/148‑1 (MQ=255)
                                                                                gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                             >  2:382079/1‑148 (MQ=255)
                                                                                 cATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                                                            >  2:228962/1‑148 (MQ=255)
                                                                                 cATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                                                            >  1:1177382/1‑148 (MQ=255)
                                                                                 cATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                                                            <  2:2370295/148‑1 (MQ=255)
                                                                                 cATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                                                            <  2:36457/148‑1 (MQ=255)
                                                                                  aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGga                                                           <  2:3062588/148‑1 (MQ=255)
                                                                                  aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGga                                                           <  2:3524873/148‑1 (MQ=255)
                                                                                  aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGga                                                           <  2:3524830/148‑1 (MQ=255)
                                                                                    aaCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatg                                                         <  2:254078/148‑1 (MQ=255)
                                                                                      cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                       >  1:868095/1‑148 (MQ=255)
                                                                                      cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                       >  1:1204521/1‑148 (MQ=255)
                                                                                       gTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTTACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATa                                                      <  2:2014229/148‑1 (MQ=255)
                                                                                        tGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATaa                                                     >  2:2621238/1‑148 (MQ=255)
                                                                                                  tACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGt                                           >  1:2601898/1‑148 (MQ=255)
                                                                                                     ggTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTa                                        <  2:3797461/148‑1 (MQ=255)
                                                                                                       tGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAgg                                      <  1:2294748/148‑1 (MQ=255)
                                                                                                         aTTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                    >  2:1460670/1‑148 (MQ=255)
                                                                                                         aTTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                    >  1:827690/1‑148 (MQ=255)
                                                                                                           tGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTc                                  <  2:849075/148‑1 (MQ=255)
                                                                                                           tGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTc                                  <  1:2722710/148‑1 (MQ=255)
                                                                                                            gCCGCTCATGAAGGCAAAGCGCCGAAGCAGTTCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGATTGGATGATAACGCCGCAGGTTTAGGTGTCa                                 <  1:2631537/148‑1 (MQ=255)
                                                                                                                   aTGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGa                          <  1:2104081/148‑1 (MQ=255)
                                                                                                                   aTGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGa                          <  2:1815098/148‑1 (MQ=255)
                                                                                                                   aTGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGa                          >  1:3399885/1‑148 (MQ=255)
                                                                                                                   aTGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGa                          >  1:3562562/1‑148 (MQ=255)
                                                                                                                       aGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTg                      >  1:2017744/1‑148 (MQ=255)
                                                                                                                       aGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTg                      >  1:1196301/1‑148 (MQ=255)
                                                                                                                        ggCAACGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTgg                     <  1:336580/148‑1 (MQ=255)
                                                                                                                        ggCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTgg                     <  1:2075028/148‑1 (MQ=255)
                                                                                                                        ggCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTgg                     >  1:3566148/1‑148 (MQ=255)
                                                                                                                         gCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGc                    <  2:2809136/148‑1 (MQ=255)
                                                                                                                         gCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGc                    >  2:1816589/1‑148 (MQ=255)
                                                                                                                         gCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGc                    >  2:2446347/1‑148 (MQ=255)
                                                                                                                          cAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCa                   >  1:133708/1‑148 (MQ=255)
                                                                                                                             aGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                >  2:689245/1‑148 (MQ=255)
                                                                                                                             aGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                >  2:689496/1‑148 (MQ=255)
                                                                                                                             aGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                <  2:170362/148‑1 (MQ=255)
                                                                                                                             aGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                >  2:2742653/1‑148 (MQ=255)
                                                                                                                              gcgcCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAAc               >  1:538025/1‑148 (MQ=255)
                                                                                                                              gcgcCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAAc               >  2:2689526/1‑148 (MQ=255)
                                                                                                                                 ccgcAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGcc            <  2:425929/148‑1 (MQ=255)
                                                                                                                                      gcagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaa       >  2:2562917/1‑148 (MQ=255)
                                                                                                                                      gcagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaa       >  1:2234036/1‑148 (MQ=255)
                                                                                                                                       cagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaa      <  2:338944/148‑1 (MQ=255)
                                                                                                                                       cagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaa      <  2:1491004/148‑1 (MQ=255)
                                                                                                                                          atcataATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATa   >  1:539056/1‑148 (MQ=255)
                                                                                                                                          atcataATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATa   <  1:808945/148‑1 (MQ=255)
                                                                                                                                           tcataaTTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATAc  >  1:3325021/1‑148 (MQ=255)
                                                                                                                                           tcataaTTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATAc  >  2:3027442/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATAC  >  CP009273/2870135‑2870421

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: