Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 185,393 A→T 87.7% F201I (TTC→ATC)  map ← methionine aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273185,3930AT87.7% 576.4 / 56.4 198F201I (TTC→ATC) mapmethionine aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (11/14);  new base T (84/89);  total (95/103)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

GTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGAATCGC  >  CP009273/185248‑185540
                                                                                                                                                 |                                                                                                                                                   
gTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCg                                                                                                                                                   <  1:2790196/148‑1 (MQ=255)
 tCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGt                                                                                                                                                  >  2:716392/1‑148 (MQ=255)
  cAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTc                                                                                                                                                 <  2:1762507/148‑1 (MQ=255)
  cAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTc                                                                                                                                                 <  1:522554/148‑1 (MQ=255)
     aaTTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATc                                                                                                                                              <  1:1594476/148‑1 (MQ=255)
       tttCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATccc                                                                                                                                            >  2:450410/1‑148 (MQ=255)
        ttCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCa                                                                                                                                           >  2:178466/1‑148 (MQ=255)
             aGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGttt                                                                                                                                      <  1:1836514/148‑1 (MQ=255)
              gCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTc                                                                                                                                     >  1:1685787/1‑148 (MQ=255)
               ccGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCa                                                                                                                                    >  2:1826677/1‑148 (MQ=255)
                  ttATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTa                                                                                                                                 >  1:392008/1‑148 (MQ=255)
                   tATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTac                                                                                                                                >  2:2622898/1‑148 (MQ=255)
                    aTCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTacg                                                                                                                               <  2:2193985/148‑1 (MQ=255)
                         tCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGtt                                                                                                                          >  2:1528944/1‑148 (MQ=255)
                         tCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGtt                                                                                                                          >  1:461040/1‑148 (MQ=255)
                          caccacAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTg                                                                                                                         >  1:1234357/1‑148 (MQ=255)
                          caccacAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTg                                                                                                                         >  1:3061717/1‑148 (MQ=255)
                           accacAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTgg                                                                                                                        <  1:2069566/148‑1 (MQ=255)
                           accacAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTgg                                                                                                                        >  1:2310852/1‑148 (MQ=255)
                            ccacAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGt                                                                                                                       <  1:3281622/148‑1 (MQ=255)
                             cacaATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGtt                                                                                                                      <  2:755504/148‑1 (MQ=255)
                              acaATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGttt                                                                                                                     <  1:374907/148‑1 (MQ=255)
                                  tAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACg                                                                                                                 >  1:1961637/1‑148 (MQ=255)
                                       tGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGt                                                                                                            <  1:2840609/148‑1 (MQ=255)
                                          tCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACTACGTTGGTTTCACGGGAGTCAt                                                                                                         <  2:1638838/148‑1 (MQ=255)
                                          tCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCAt                                                                                                         >  1:932044/1‑148 (MQ=255)
                                             tatTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGt                                                                                                      >  1:2471223/1‑148 (MQ=255)
                                              atTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTg                                                                                                     <  1:2699017/148‑1 (MQ=255)
                                               ttGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTgc                                                                                                    >  1:126399/1‑148 (MQ=255)
                                                   gCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTgcagca                                                                                                <  2:1904935/148‑1 (MQ=255)
                                                    cAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCAc                                                                                               >  2:294702/1‑148 (MQ=255)
                                                       aCAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTg                                                                                            >  2:491122/1‑148 (MQ=255)
                                                          aGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCgg                                                                                         >  1:2346549/1‑148 (MQ=255)
                                                            cTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGtt                                                                                       >  1:1344655/1‑148 (MQ=255)
                                                                 aTCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCa                                                                                  <  2:945336/148‑1 (MQ=255)
                                                                 aTCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCa                                                                                  >  1:2723876/1‑148 (MQ=255)
                                                                 aTCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTGTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCa                                                                                  <  2:2174372/148‑1 (MQ=255)
                                                                    tttGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATgg                                                                               >  2:2902571/1‑148 (MQ=255)
                                                                    tttGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATgg                                                                               >  2:2769982/1‑148 (MQ=255)
                                                                     gtGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGa                                                                              <  2:2477191/147‑1 (MQ=255)
                                                                      tGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGCTTACAGAACGACGTTGGTTTCACGGGAGTCATAGTGAAGAACCAGCGGTTCTTCATGGaa                                                                             >  2:1211691/1‑148 (MQ=255)
                                                                       ggTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAg                                                                            <  1:2241952/148‑1 (MQ=255)
                                                                        gTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGc                                                                           <  1:2546768/148‑1 (MQ=255)
                                                                        gTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGc                                                                           <  2:1026538/148‑1 (MQ=255)
                                                                         tttttACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGcc                                                                          >  1:3006759/1‑148 (MQ=255)
                                                                          ttttACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCcg                                                                         >  2:254912/1‑148 (MQ=255)
                                                                             tACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGa                                                                      <  2:2114515/148‑1 (MQ=255)
                                                                                  tCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAAt                                                                 <  2:177335/148‑1 (MQ=255)
                                                                                   ccAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATa                                                                <  2:1751123/148‑1 (MQ=255)
                                                                                      gCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCg                                                             <  1:3107367/148‑1 (MQ=255)
                                                                                      gCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCg                                                             <  2:125028/148‑1 (MQ=255)
                                                                                       ccATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCgt                                                            <  2:2928522/148‑1 (MQ=255)
                                                                                       ccATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCgt                                                            >  2:634973/1‑148 (MQ=255)
                                                                                        cATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCgtg                                                           >  1:2911903/1‑148 (MQ=255)
                                                                                          tCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTc                                                         >  2:579807/1‑148 (MQ=255)
                                                                                            tttCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCg                                                       >  2:1682826/1‑148 (MQ=255)
                                                                                              tCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCa                                                     <  1:2104233/148‑1 (MQ=255)
                                                                                              tCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCa                                                     <  1:2459115/148‑1 (MQ=255)
                                                                                                   gTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATAtt                                                <  1:1479091/148‑1 (MQ=255)
                                                                                                     gCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCa                                              <  1:657284/148‑1 (MQ=255)
                                                                                                        gATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacga                                           <  1:878947/148‑1 (MQ=255)
                                                                                                        gATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacga                                           <  2:1865718/148‑1 (MQ=255)
                                                                                                        gATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacga                                           <  1:575038/148‑1 (MQ=255)
                                                                                                        gATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacga                                           >  1:1287511/1‑148 (MQ=255)
                                                                                                         aTCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCACTATTCacgaa                                          >  2:3142981/1‑148 (MQ=255)
                                                                                                         aTCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacgaa                                          <  1:1488021/148‑1 (MQ=255)
                                                                                                         aTCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacgaa                                          >  2:1244539/1‑148 (MQ=255)
                                                                                                           ctctTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacgaacg                                        <  2:317254/148‑1 (MQ=255)
                                                                                                           ctctTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacgaacg                                        <  2:666390/148‑1 (MQ=255)
                                                                                                             ctTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAacgac                                      <  1:683374/148‑1 (MQ=255)
                                                                                                             ctTTTTTACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTCTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAacgac                                      <  1:450410/148‑1 (MQ=255)
                                                                                                              ttttttACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCACTATTCACGAacgacg                                     >  2:2834313/1‑148 (MQ=255)
                                                                                                              ttttttACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAacgacg                                     <  2:609508/148‑1 (MQ=255)
                                                                                                               tttttACCCGCGTTTACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGTTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACgg                                    <  1:1978860/148‑1 (MQ=255)
                                                                                                               tttttACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACgg                                    <  1:1890122/148‑1 (MQ=255)
                                                                                                                 tttACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGgag                                  >  2:1122111/1‑148 (MQ=255)
                                                                                                                 tttACCCGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGgag                                  <  1:293694/148‑1 (MQ=255)
                                                                                                                     cccGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGc                              >  1:2078972/1‑148 (MQ=255)
                                                                                                                     cccGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGc                              <  1:265506/148‑1 (MQ=255)
                                                                                                                      ccGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGcc                             <  1:2902571/148‑1 (MQ=255)
                                                                                                                      ccGCGTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGcc                             <  1:3237534/148‑1 (MQ=255)
                                                                                                                        gcgTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCtt                           <  1:2277434/148‑1 (MQ=255)
                                                                                                                        gcgTTGACCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTAATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCtt                           <  1:79121/148‑1 (MQ=255)
                                                                                                                           ttGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAAGACCGGGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTg                        <  2:2029996/148‑1 (MQ=255)
                                                                                                                              aCCATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCtt                     <  2:782983/148‑1 (MQ=255)
                                                                                                                               ccATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTc                    >  1:1018072/1‑148 (MQ=255)
                                                                                                                                cATTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCg                   >  2:1342330/1‑148 (MQ=255)
                                                                                                                                 aTTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGa                  <  2:216554/148‑1 (MQ=255)
                                                                                                                                 aTTGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGa                  >  2:1033871/1‑148 (MQ=255)
                                                                                                                                  ttGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGAc                 >  2:1720208/1‑148 (MQ=255)
                                                                                                                                  ttGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGAc                 <  2:89566/148‑1 (MQ=255)
                                                                                                                                  ttGGCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGAc                 <  2:89814/148‑1 (MQ=255)
                                                                                                                                    ggCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACaa               <  1:3191795/148‑1 (MQ=255)
                                                                                                                                     gCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACaaa              >  2:11342/1‑148 (MQ=255)
                                                                                                                                     gCTCGATGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACaaa              >  1:1213607/1‑148 (MQ=255)
                                                                                                                                          aTGGTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCt         >  1:616832/1‑148 (MQ=255)
                                                                                                                                            ggTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGa       >  2:3019237/1‑148 (MQ=255)
                                                                                                                                            ggTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGa       >  2:548174/1‑148 (MQ=255)
                                                                                                                                            ggTGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGa       <  2:1464821/148‑1 (MQ=255)
                                                                                                                                              tGATCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGAAt     >  1:174190/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCCAGCCATCTTTCATGGTGCGGATCTCTTTTTTACCCGCGTTGACCATTGGCTCGATGGTGAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATTTCTGAATCGC  >  CP009273/185248‑185540

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: ‑