Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,870,278 C→A 88.4% 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,2780CA88.4% 440.9 / 57.8 181I113I (ATC→ATAiapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base C (14/7);  new base A (89/71);  total (103/78)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 192 total aligned reads displayed.

TCAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACG  >  CP009273/2870134‑2870422
                                                                                                                                                |                                                                                                                                                
tCAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAAtt                                                                                                                                               <  1:438086/148‑1 (MQ=255)
 cAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTa                                                                                                                                              <  2:2299341/148‑1 (MQ=255)
 cAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTa                                                                                                                                              <  2:3050620/148‑1 (MQ=255)
    cAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATgg                                                                                                                                           >  2:1827509/1‑148 (MQ=255)
    cAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATgg                                                                                                                                           <  2:42835/148‑1 (MQ=255)
       aTGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGgcgc                                                                                                                                        >  1:3074649/1‑148 (MQ=255)
         gggTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCAt                                                                                                                                      <  1:2582112/148‑1 (MQ=255)
         gggTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCAt                                                                                                                                      <  1:2582152/148‑1 (MQ=255)
             tATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTgg                                                                                                                                  >  1:834717/1‑148 (MQ=255)
             tATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTgg                                                                                                                                  <  1:1363707/148‑1 (MQ=255)
                cGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATa                                                                                                                               <  2:1708113/148‑1 (MQ=255)
                  cAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACt                                                                                                                             >  2:2383881/1‑148 (MQ=255)
                   aGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACtt                                                                                                                            <  2:338920/148‑1 (MQ=255)
                    gTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTa                                                                                                                           >  1:2214430/1‑148 (MQ=255)
                     tGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTAc                                                                                                                          >  1:345101/1‑148 (MQ=255)
                     tGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTAc                                                                                                                          >  1:2010765/1‑148 (MQ=255)
                      gATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACg                                                                                                                         <  1:1733160/148‑1 (MQ=255)
                           tCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCg                                                                                                                    <  2:274328/148‑1 (MQ=255)
                           tCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCg                                                                                                                    >  2:1226794/1‑148 (MQ=255)
                              gACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTg                                                                                                                 >  1:68194/1‑148 (MQ=255)
                               aCATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGa                                                                                                                >  1:1456949/1‑148 (MQ=255)
                               aCATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGa                                                                                                                >  2:2274912/1‑148 (MQ=255)
                                cATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAg                                                                                                               <  2:147671/148‑1 (MQ=255)
                                cATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAg                                                                                                               >  1:241418/1‑148 (MQ=255)
                                  tttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                                                             >  2:1760250/1‑148 (MQ=255)
                                  tttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                                                             >  2:1271049/1‑148 (MQ=255)
                                  tttAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                                                             >  1:1408011/1‑148 (MQ=255)
                                     aaTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATg                                                                                                          <  2:946415/148‑1 (MQ=255)
                                      aTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGc                                                                                                         <  1:1027818/148‑1 (MQ=255)
                                      aTAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGc                                                                                                         <  1:423356/148‑1 (MQ=255)
                                       tAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                                                        <  1:3100445/148‑1 (MQ=255)
                                          tCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGAc                                                                                                     >  2:1286128/1‑148 (MQ=255)
                                               atatTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGa                                                                                                >  1:2675369/1‑148 (MQ=255)
                                                   ttATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGcc                                                                                            >  1:1164829/1‑148 (MQ=255)
                                                    tataCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCa                                                                                           >  2:2607104/1‑148 (MQ=255)
                                                    tataCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCa                                                                                           <  1:2470755/148‑1 (MQ=255)
                                                     ataCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCaa                                                                                          <  2:2021483/148‑1 (MQ=255)
                                                      taCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAt                                                                                         >  2:1396458/1‑148 (MQ=255)
                                                      taCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAt                                                                                         >  1:539647/1‑148 (MQ=255)
                                                       aCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAtc                                                                                        >  2:1349830/1‑148 (MQ=255)
                                                       aCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAtc                                                                                        >  1:3049228/1‑148 (MQ=255)
                                                       aCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAtc                                                                                        >  2:2997270/1‑148 (MQ=255)
                                                        ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAtct                                                                                       >  2:722309/1‑148 (MQ=255)
                                                         cgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAtctc                                                                                      >  2:1303315/1‑148 (MQ=255)
                                                          gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcg                                                                                     <  2:317397/148‑1 (MQ=255)
                                                          gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcg                                                                                     <  2:317564/148‑1 (MQ=255)
                                                             cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggcg                                                                                  <  1:1801687/148‑1 (MQ=255)
                                                             cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggcg                                                                                  >  2:1083436/1‑148 (MQ=255)
                                                             cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggcg                                                                                  >  2:2651407/1‑148 (MQ=255)
                                                                 aTAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                              >  1:960978/1‑148 (MQ=255)
                                                                 aTAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                              >  2:3088911/1‑148 (MQ=255)
                                                                  tAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                             <  2:517301/148‑1 (MQ=255)
                                                                     tCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                          >  2:82738/1‑148 (MQ=255)
                                                                         aaGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTAc                                                                      >  1:1797903/1‑148 (MQ=255)
                                                                          agagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                     >  1:188182/1‑148 (MQ=255)
                                                                          agagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                     >  2:224155/1‑148 (MQ=255)
                                                                           gagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACaa                                                                    <  1:1531739/148‑1 (MQ=255)
                                                                            agCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAg                                                                   >  2:2282053/1‑148 (MQ=255)
                                                                             gCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGATCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                                                                  <  1:3114831/148‑1 (MQ=255)
                                                                               tGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                <  1:579528/148‑1 (MQ=255)
                                                                                 gCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                              >  2:1307413/1‑148 (MQ=255)
                                                                                   aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGga                                                            >  1:1247927/1‑148 (MQ=255)
                                                                                    tAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgat                                                           >  1:2391814/1‑148 (MQ=255)
                                                                                     aaCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatg                                                          <  1:1894286/148‑1 (MQ=255)
                                                                                      aCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatga                                                         <  2:3074649/148‑1 (MQ=255)
                                                                                       cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                        <  2:2876303/148‑1 (MQ=255)
                                                                                       cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                                                        >  1:1054484/1‑148 (MQ=255)
                                                                                        gTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATa                                                       <  2:2851833/148‑1 (MQ=255)
                                                                                              ggAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccg                                                 <  1:2282053/148‑1 (MQ=255)
                                                                                                   tACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGt                                            >  1:1663200/1‑148 (MQ=255)
                                                                                                     cGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGttt                                          <  2:2191672/148‑1 (MQ=255)
                                                                                                     cGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGttt                                          <  2:15441/148‑1 (MQ=255)
                                                                                                      ggTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTa                                         >  2:1836412/1‑148 (MQ=255)
                                                                                                         gATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgt                                      <  2:600738/148‑1 (MQ=255)
                                                                                                         gATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgt                                      <  2:899973/148‑1 (MQ=255)
                                                                                                         gATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgt                                      <  2:3209115/148‑1 (MQ=255)
                                                                                                          aTTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                     <  2:68194/148‑1 (MQ=255)
                                                                                                          aTTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                     >  2:2809792/1‑148 (MQ=255)
                                                                                                             gCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCa                                  <  1:1307937/148‑1 (MQ=255)
                                                                                                             gCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCa                                  <  1:459892/148‑1 (MQ=255)
                                                                                                               ctctaaTGAAGGCAAAGCGCCGCAGCATATCATATTTTTGGCGAATCTGGATACTTACGCCCCGCTGCGCGATGATGACGCAGATGCCAATCTGGGCGGGCTGACGTTACATGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                <  1:1175614/143‑1 (MQ=255)
                                                                                                               cGCTCATGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                <  1:2606002/148‑1 (MQ=255)
                                                                                                               cGCTCATGAAGGAAAAGCGCCGCAGCAGATCATAATTAGGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                <  1:2775593/148‑1 (MQ=255)
                                                                                                                    aTGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGa                           <  2:1797898/148‑1 (MQ=255)
                                                                                                                     tGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGaa                          >  1:1674293/1‑148 (MQ=255)
                                                                                                                     tGAAGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGaa                          >  1:3170188/1‑148 (MQ=255)
                                                                                                                       aaGGCAAAGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAAtt                        >  1:2779770/1‑148 (MQ=255)
                                                                                                                             aaGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGACTGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGa                  >  2:1348707/1‑148 (MQ=255)
                                                                                                                              aGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                 >  1:1773610/1‑148 (MQ=255)
                                                                                                                              aGCGCCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                 >  2:3285074/1‑148 (MQ=255)
                                                                                                                               gcgcCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAAc                >  1:715298/1‑148 (MQ=255)
                                                                                                                               gcgcCGCAGCAGATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAAc                >  2:2490670/1‑148 (MQ=255)
                                                                                                                                     cagcagATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGa          >  2:2759190/1‑148 (MQ=255)
                                                                                                                                         agATCATAATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaaaa      <  1:2464235/148‑1 (MQ=255)
                                                                                                                                           atcataATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATa    >  2:1839922/1‑148 (MQ=255)
                                                                                                                                            tcataaTTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATAc   >  2:319199/1‑148 (MQ=255)
                                                                                                                                             cataatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGTCAGAACGCCTGAAAAATAct  >  2:1946016/1‑147 (MQ=255)
                                                                                                                                             cataatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACg  >  2:3275846/1‑148 (MQ=255)
                                                                                                                                             cataatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACg  >  2:484957/1‑148 (MQ=255)
                                                                                                                                             cataatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACg  >  2:1637895/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TCAGCAAATGGGTTATCGCAGTGATATTCGGACATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACG  >  CP009273/2870134‑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: