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.9% 606.6 / 14.7 203I113I (ATC→ATAiapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base A (114/80);  total (118/85)
Fisher's exact test for biased strand distribution p-value = 4.96e-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 207 total aligned reads displayed.

AGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACG  >  CP009273/2870136‑2870422
                                                                                                                                              |                                                                                                                                                
 gCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATCACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATg                                                                                                                                            >  2:2194132/1‑148 (MQ=255)
 gCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATg                                                                                                                                            >  2:1666596/1‑148 (MQ=255)
  cAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATgg                                                                                                                                           >  2:830220/1‑148 (MQ=255)
   aaaTGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGgc                                                                                                                                          >  1:692828/1‑148 (MQ=255)
    aaTGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGgcg                                                                                                                                         >  1:2884535/1‑148 (MQ=255)
      tGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCa                                                                                                                                       >  1:1049178/1‑148 (MQ=255)
        ggTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATc                                                                                                                                     <  2:25782/148‑1 (MQ=255)
          ttATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTg                                                                                                                                   >  2:2018599/1‑148 (MQ=255)
            aTCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGAATCTGGa                                                                                                                                 <  1:1556916/148‑1 (MQ=255)
             tCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGAt                                                                                                                                >  1:160720/1‑148 (MQ=255)
             tCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGAt                                                                                                                                >  2:1417827/1‑148 (MQ=255)
              cGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATa                                                                                                                               <  2:2738923/148‑1 (MQ=255)
                 aGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACtt                                                                                                                            <  1:3320223/148‑1 (MQ=255)
                   tGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTAc                                                                                                                          >  2:473105/1‑148 (MQ=255)
                    gATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACg                                                                                                                         >  2:1840164/1‑148 (MQ=255)
                     atatTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGc                                                                                                                        >  1:2269466/1‑148 (MQ=255)
                     atatTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGc                                                                                                                        >  1:3348937/1‑148 (MQ=255)
                     atatTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGc                                                                                                                        <  2:1459744/148‑1 (MQ=255)
                      tatTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGcc                                                                                                                       >  1:1978509/1‑148 (MQ=255)
                       atTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGccc                                                                                                                      >  2:1584054/1‑148 (MQ=255)
                       atTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGccc                                                                                                                      >  1:3134320/1‑148 (MQ=255)
                       atTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGccc                                                                                                                      >  2:3431004/1‑148 (MQ=255)
                        ttCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGcccc                                                                                                                     >  1:56126/1‑148 (MQ=255)
                             aCATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGa                                                                                                                >  2:91228/1‑148 (MQ=255)
                              cATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAg                                                                                                               <  2:3497257/148‑1 (MQ=255)
                              cATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAg                                                                                                               >  2:3358991/1‑148 (MQ=255)
                              cATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAg                                                                                                               <  2:3293500/148‑1 (MQ=255)
                                 ttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGa                                                                                                            >  2:836771/1‑148 (MQ=255)
                                 ttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGa                                                                                                            >  2:1183230/1‑148 (MQ=255)
                                  tAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGAt                                                                                                           <  1:451250/148‑1 (MQ=255)
                                   aaTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATg                                                                                                          <  2:2852811/148‑1 (MQ=255)
                                    aTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGc                                                                                                         <  2:3319828/148‑1 (MQ=255)
                                    aTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGc                                                                                                         >  2:2979889/1‑148 (MQ=255)
                                    aTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGc                                                                                                         <  2:2422281/148‑1 (MQ=255)
                                     tAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                                                        >  2:1548255/1‑148 (MQ=255)
                                     tAGTCGGCATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                                                        >  1:929081/1‑148 (MQ=255)
                                      aGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTg                                                                                                       <  2:1080177/148‑1 (MQ=255)
                                      aGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTg                                                                                                       <  2:567657/148‑1 (MQ=255)
                                         cGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACg                                                                                                    <  1:2054000/148‑1 (MQ=255)
                                          ggTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGc                                                                                                   <  1:649708/148‑1 (MQ=255)
                                             atatTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGa                                                                                                >  1:2025975/1‑148 (MQ=255)
                                             atatTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGa                                                                                                >  1:3260291/1‑148 (MQ=255)
                                               atTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATg                                                                                              <  1:3179389/148‑1 (MQ=255)
                                                    gaCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAt                                                                                         <  2:2211430/147‑1 (MQ=255)
                                                      ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAtct                                                                                       >  1:595679/1‑148 (MQ=255)
                                                      ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAtct                                                                                       >  2:1041335/1‑148 (MQ=255)
                                                        gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcg                                                                                     <  1:3005446/148‑1 (MQ=255)
                                                           cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggcg                                                                                  <  1:122379/148‑1 (MQ=255)
                                                           cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggcg                                                                                  <  1:1897604/148‑1 (MQ=255)
                                                             cgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCggg                                                                                <  1:1652188/148‑1 (MQ=255)
                                                                tAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                             >  2:1706293/1‑148 (MQ=255)
                                                                tAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                             >  2:3170162/1‑148 (MQ=255)
                                                                tAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                             >  2:647519/1‑148 (MQ=255)
                                                                tAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                             <  1:998080/148‑1 (MQ=255)
                                                                  aTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGAc                                                                           >  2:1722448/1‑148 (MQ=255)
                                                                      tAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTa                                                                       <  1:335413/148‑1 (MQ=255)
                                                                       aaGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTAc                                                                      >  2:3243588/1‑148 (MQ=255)
                                                                       aaGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATACCAATCTCGGCGGGCTGACGTTAc                                                                      >  1:3163750/1‑148 (MQ=255)
                                                                        agagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                     >  1:78702/1‑148 (MQ=255)
                                                                         gagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACaa                                                                    <  1:3206979/148‑1 (MQ=255)
                                                                          agCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAg                                                                   >  2:229354/1‑148 (MQ=255)
                                                                          agCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAg                                                                   >  1:1239155/1‑148 (MQ=255)
                                                                           gCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                                                                  <  1:3210970/148‑1 (MQ=255)
                                                                           gCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                                                                  <  1:1711382/148‑1 (MQ=255)
                                                                            cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                 <  1:993836/148‑1 (MQ=255)
                                                                             tGGCATAACGTGACGGGAAGTTCGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                <  2:3087459/148‑1 (MQ=255)
                                                                               gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                              >  2:119590/1‑148 (MQ=255)
                                                                               gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                              >  2:3460597/1‑148 (MQ=255)
                                                                               gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                              >  1:536970/1‑148 (MQ=255)
                                                                               gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                              <  1:1589388/148‑1 (MQ=255)
                                                                               gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                              >  1:484042/1‑148 (MQ=255)
                                                                               gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                              >  2:2265685/1‑148 (MQ=255)
                                                                               gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                              >  1:2301353/1‑148 (MQ=255)
                                                                                      gTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATa                                                       <  1:1350090/148‑1 (MQ=255)
                                                                                            ggAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccg                                                 >  1:2303380/1‑148 (MQ=255)
                                                                                               aGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAg                                              <  1:2005364/148‑1 (MQ=255)
                                                                                                gTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAgg                                             >  1:225613/1‑148 (MQ=255)
                                                                                                 tACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGt                                            <  1:2193582/148‑1 (MQ=255)
                                                                                                 tACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGt                                            >  1:861436/1‑148 (MQ=255)
                                                                                                  aCGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGtt                                           >  1:168169/1‑148 (MQ=255)
                                                                                                  aCGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGtt                                           <  1:242275/148‑1 (MQ=255)
                                                                                                     gTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCTTAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAg                                        >  1:1445252/1‑148 (MQ=255)
                                                                                                        aTTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                     >  1:1074180/1‑148 (MQ=255)
                                                                                                           gCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCa                                  <  2:876372/148‑1 (MQ=255)
                                                                                                                  aTGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGa                           >  1:1959419/1‑148 (MQ=255)
                                                                                                                  aTGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGa                           <  1:428353/148‑1 (MQ=255)
                                                                                                                   tGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGaa                          <  1:3219496/148‑1 (MQ=255)
                                                                                                                      aGGCAAAGCGCCGCAGCAGATCATAATTAGGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTg                       <  1:549905/148‑1 (MQ=255)
                                                                                                                           aaGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGa                  >  2:2352406/1‑148 (MQ=255)
                                                                                                                               gccgcAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGc              <  2:692828/148‑1 (MQ=255)
                                                                                                                                ccgcAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGcc             <  1:2804892/148‑1 (MQ=255)
                                                                                                                                ccgcAGCAGATCATAATTAGGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGcc             <  1:2567679/148‑1 (MQ=255)
                                                                                                                                    agcagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaa         >  1:842839/1‑148 (MQ=255)
                                                                                                                                      cagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaa       >  1:1009310/1‑148 (MQ=255)
                                                                                                                                      cagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaa       >  1:159005/1‑148 (MQ=255)
                                                                                                                                        gATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAAt     >  1:1494230/1‑148 (MQ=255)
                                                                                                                                        gATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAAt     >  1:3236106/1‑148 (MQ=255)
                                                                                                                                        gATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAAt     >  1:1493882/1‑148 (MQ=255)
                                                                                                                                           cataatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACg  >  1:3441283/1‑148 (MQ=255)
                                                                                                                                           cataatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACg  >  1:2955463/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                
AGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACG  >  CP009273/2870136‑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: