Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,049,645 G→A 11.6% E16K (GAA→AAA)  fau → 5‑formyltetrahydrofolate cyclo‑ligase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,049,6450GA11.6% 858.0 / 84.3 266E16K (GAA→AAA) fau5‑formyltetrahydrofolate cyclo‑ligase family protein
Reads supporting (aligned to +/- strand):  ref base G (124/111);  new base A (17/14);  total (141/125)
Fisher's exact test for biased strand distribution p-value = 8.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 270 total aligned reads displayed.

CAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTT  >  CP009273/3049498‑3049794
                                                                                                                                                   |                                                                                                                                                     
    cTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATgg                                                                                                                                                 >  1:2256716/1‑150 (MQ=255)
        ggcggcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTca                                                                                                                                             >  1:894982/1‑150 (MQ=255)
        ggcggcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATGATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTca                                                                                                                                             >  2:2731051/1‑150 (MQ=255)
          cggcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTcaac                                                                                                                                           <  2:2285324/150‑1 (MQ=255)
          cggcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTcaac                                                                                                                                           <  2:288002/150‑1 (MQ=255)
            gcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTcaacaa                                                                                                                                         <  1:644019/150‑1 (MQ=255)
             cATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAg                                                                                                                                        >  1:2367494/1‑150 (MQ=255)
               tctcGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGcc                                                                                                                                      >  2:4159166/1‑150 (MQ=255)
                  cGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCt                                                                                                                                   >  2:319435/1‑150 (MQ=255)
                  cGGAGATTCCCTTCATATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCt                                                                                                                                   >  1:2830961/1‑150 (MQ=255)
                    gagaTTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAc                                                                                                                                 >  2:4317257/1‑150 (MQ=255)
                     agaTTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAcc                                                                                                                                >  2:5117668/1‑150 (MQ=255)
                     agaTTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATGATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAcc                                                                                                                                >  1:190486/1‑150 (MQ=255)
                     agaTTCCCTTCTTATCTGGCACCAGCCATAACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGACATGGGTCAACAAGCCGCTAcc                                                                                                                                >  2:546202/1‑150 (MQ=255)
                        ttCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCgg                                                                                                                             <  2:3144074/150‑1 (MQ=255)
                              tcttATCTGGCACCAGCCTTGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGatgatg                                                                                                                       <  2:2103950/150‑1 (MQ=255)
                               cttATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGAtggtga                                                                                                                      >  1:3943508/1‑150 (MQ=255)
                                    cTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTAt                                                                                                                 >  2:176873/1‑150 (MQ=255)
                                     tGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATc                                                                                                                >  2:2404748/1‑150 (MQ=255)
                                        cACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATcccc                                                                                                             <  1:918799/150‑1 (MQ=255)
                                         aCCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATccccc                                                                                                            <  1:774226/150‑1 (MQ=255)
                                         aCCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATccccc                                                                                                            <  1:3351713/150‑1 (MQ=255)
                                         aCCAGCCATGACGCAACTACCAGAACTCCCACTGACATGATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATccccc                                                                                                            <  1:2078940/150‑1 (MQ=255)
                                          ccAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCg                                                                                                           <  1:2283520/150‑1 (MQ=255)
                                              ccATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggtgg                                                                                                       <  1:4982034/150‑1 (MQ=255)
                                                  gACGCAACTACCAGAACTCCCACTGACATGATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGAt                                                                                                   >  1:2868880/1‑150 (MQ=255)
                                                  gACGCAACTACCAGAACTCCCACTGACATGATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGAt                                                                                                   <  2:1590672/150‑1 (MQ=255)
                                                     gCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGc                                                                                                <  1:1160209/150‑1 (MQ=255)
                                                      cAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGca                                                                                               >  2:830983/1‑150 (MQ=255)
                                                        aCTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGcaca                                                                                             >  1:540860/1‑150 (MQ=255)
                                                         cTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACAt                                                                                            <  2:508982/150‑1 (MQ=255)
                                                          tACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATa                                                                                           <  2:3800278/150‑1 (MQ=255)
                                                          tACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATa                                                                                           <  2:462402/150‑1 (MQ=255)
                                                           aCCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATAc                                                                                          >  2:3476284/1‑150 (MQ=255)
                                                            ccAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACg                                                                                         <  1:5112251/150‑1 (MQ=255)
                                                            ccAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACg                                                                                         >  1:1012259/1‑150 (MQ=255)
                                                               gAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTc                                                                                      <  1:1714184/150‑1 (MQ=255)
                                                                  cTCCCACTGACATGATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCt                                                                                   <  1:3819106/150‑1 (MQ=255)
                                                                   tCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGAGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTCTCCCCCGGTGGTGATGGCACATACGGTCGCTg                                                                                  >  2:4253559/1‑150 (MQ=255)
                                                                    cccACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGt                                                                                 <  1:3311962/150‑1 (MQ=255)
                                                                    cccACTGACATGATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGt                                                                                 >  1:4361914/1‑150 (MQ=255)
                                                                      cACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTAt                                                                               >  2:4764191/1‑150 (MQ=255)
                                                                          gACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCct                                                                           <  1:3198149/150‑1 (MQ=255)
                                                                          gACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCct                                                                           >  2:457578/1‑150 (MQ=255)
                                                                            cATGATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCctct                                                                         <  1:4952814/150‑1 (MQ=255)
                                                                              ttATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCctctct                                                                       >  2:1103443/1‑150 (MQ=255)
                                                                               tATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCtt                                                                      <  1:4808786/150‑1 (MQ=255)
                                                                                  cccGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTg                                                                   >  1:3263728/1‑150 (MQ=255)
                                                                                    cGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGAt                                                                 >  1:4244291/1‑150 (MQ=255)
                                                                                      aCAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTGAttg                                                               >  1:384427/1‑148 (MQ=255)
                                                                                       cAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGc                                                              >  1:1463205/1‑150 (MQ=255)
                                                                                           aaaTCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAAc                                                          >  2:2060348/1‑150 (MQ=255)
                                                                                           aaaTCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAAc                                                          <  2:4864488/150‑1 (MQ=255)
                                                                                            aaTCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACt                                                         >  1:270815/1‑150 (MQ=255)
                                                                                            aaTCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACt                                                         <  1:75317/150‑1 (MQ=255)
                                                                                            aaTCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACt                                                         <  1:2309171/150‑1 (MQ=255)
                                                                                             aTCCGCATAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTTGTGATGGCACATACGGTCGCTGTTTTCCTCTCTTTTGATGGCGAACTc                                                        <  2:3772677/150‑1 (MQ=255)
                                                                                              tCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCg                                                       <  2:492127/150‑1 (MQ=255)
                                                                                                cGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGac                                                     <  2:4473667/150‑1 (MQ=255)
                                                                                                 gCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGaca                                                    >  2:20323/1‑150 (MQ=255)
                                                                                                  cAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGacac                                                   >  2:3535353/1‑150 (MQ=255)
                                                                                                  cAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGacac                                                   >  2:3535932/1‑150 (MQ=255)
                                                                                                   aaaaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACAcc                                                  >  2:3865263/1‑150 (MQ=255)
                                                                                                         aTTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCa                                            >  1:5078254/1‑150 (MQ=255)
                                                                                                         aTTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCa                                            >  1:1005980/1‑150 (MQ=255)
                                                                                                          ttCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCAc                                           <  2:54257/150‑1 (MQ=255)
                                                                                                           tCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACt                                          >  1:3548302/1‑150 (MQ=255)
                                                                                                            cGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTc                                         >  2:4097376/1‑150 (MQ=255)
                                                                                                              gCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCAt                                       >  1:3999821/1‑150 (MQ=255)
                                                                                                               cAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATa                                      <  1:1321467/150‑1 (MQ=255)
                                                                                                                aaCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAg                                     >  2:2478242/1‑150 (MQ=255)
                                                                                                                  cgtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaa                                   >  2:5036417/1‑150 (MQ=255)
                                                                                                                  cgtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTATTTTGATTGCGAAATCGACACCCACCCCCTCAGAGaa                                   >  1:1098095/1‑150 (MQ=255)
                                                                                                                   gtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaac                                  >  2:1625423/1‑150 (MQ=255)
                                                                                                                   gtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaac                                  >  1:408426/1‑150 (MQ=255)
                                                                                                                    tcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaaca                                 <  2:1670162/150‑1 (MQ=255)
                                                                                                                    tcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGGATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaaca                                 >  2:3671646/1‑150 (MQ=255)
                                                                                                                     cgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaacaa                                >  1:4163564/1‑150 (MQ=255)
                                                                                                                      gtcgtGCGTTAACGCCGGAACAACATCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaacaac                               <  2:604854/150‑1 (MQ=255)
                                                                                                                       tcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAAct                              <  2:960012/150‑1 (MQ=255)
                                                                                                                       tcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAAct                              <  1:2451947/150‑1 (MQ=255)
                                                                                                                          tGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTg                           <  2:1026220/150‑1 (MQ=255)
                                                                                                                             gTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGgcg                        <  1:1220297/150‑1 (MQ=255)
                                                                                                                              ttAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGgcgc                       <  1:4519963/150‑1 (MQ=255)
                                                                                                                              ttAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACACGACACCCAGCCACTCATAGAACAACTCTGgcgc                       >  1:1801835/1‑150 (MQ=255)
                                                                                                                                  cGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCg                   <  2:3341929/150‑1 (MQ=255)
                                                                                                                                   gCCGGAACAACAGCAGGAACTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTAGTCCTCTCTTTTGATGGCGCACTCGACACCCAGCCACTCATAGAACAACGCTGGCGCGCCgg                  >  1:4475159/1‑150 (MQ=255)
                                                                                                                                   gCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCgg                  <  2:4893720/150‑1 (MQ=255)
                                                                                                                                   gCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCgg                  >  1:3772246/1‑150 (MQ=255)
                                                                                                                                    ccGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGt                 >  2:1437085/1‑150 (MQ=255)
                                                                                                                                    ccGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGt                 >  1:4708228/1‑150 (MQ=255)
                                                                                                                                     ctgAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTa                >  2:1058214/3‑150 (MQ=255)
                                                                                                                                      ggAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTaa               >  2:3358779/1‑150 (MQ=255)
                                                                                                                                        aacaacAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgc             <  1:3995519/150‑1 (MQ=255)
                                                                                                                                        aacaacAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgc             <  1:4390393/150‑1 (MQ=255)
                                                                                                                                        aacaacAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGGAAgc             >  2:4375095/1‑150 (MQ=255)
                                                                                                                                          caacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCCCTCATAGAACAACTCTGGCGCGCCGGTAAgcgc           >  1:1394096/1‑150 (MQ=255)
                                                                                                                                              agcagGACATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGtat       <  1:4902365/150‑1 (MQ=255)
                                                                                                                                                 agGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATc    >  2:4956865/1‑150 (MQ=255)
                                                                                                                                                   gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCtt  >  1:694519/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTT  >  CP009273/3049498‑3049794

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: