Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,049,645 G→A 33.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,6450GA33.6% 411.7 / 302.1 292E16K (GAA→AAA) fau5‑formyltetrahydrofolate cyclo‑ligase family protein
Reads supporting (aligned to +/- strand):  ref base G (101/93);  new base A (59/39);  total (160/132)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 295 total aligned reads displayed.

TACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTT  >  CP009273/3049496‑3049794
                                                                                                                                                     |                                                                                                                                                     
tACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAgg                                                                                                                                                       >  1:1441396/1‑150 (MQ=255)
 aCAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGa                                                                                                                                                      <  2:2892402/150‑1 (MQ=255)
 aCAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGa                                                                                                                                                      >  2:580725/1‑150 (MQ=255)
 aCAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGa                                                                                                                                                      >  2:3422773/1‑150 (MQ=255)
    gCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAAt                                                                                                                                                   >  2:1971144/1‑150 (MQ=255)
          ggcggcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTca                                                                                                                                             >  2:869266/1‑150 (MQ=255)
          ggcggcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTca                                                                                                                                             <  2:727546/150‑1 (MQ=255)
          ggcgacATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTca                                                                                                                                             <  1:2247741/150‑1 (MQ=255)
            cggcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTcaac                                                                                                                                           >  1:2372607/1‑150 (MQ=255)
             ggcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTcaaca                                                                                                                                          >  2:1393057/1‑150 (MQ=255)
             ggcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTcaaca                                                                                                                                          >  1:2187416/1‑150 (MQ=255)
              gcATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAACAGGAAATGGGTcaacaa                                                                                                                                         <  1:3368651/150‑1 (MQ=255)
               cATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAg                                                                                                                                        <  2:3041161/150‑1 (MQ=255)
                aTCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGc                                                                                                                                       <  1:4084478/150‑1 (MQ=255)
                aTCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGc                                                                                                                                       <  2:482928/150‑1 (MQ=255)
                aTCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGc                                                                                                                                       >  2:3443148/1‑150 (MQ=255)
                aTCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGc                                                                                                                                       >  1:2573164/1‑150 (MQ=255)
                 tctcGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGcc                                                                                                                                      <  2:3354686/150‑1 (MQ=255)
                   tcGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGc                                                                                                                                    <  2:2862013/150‑1 (MQ=255)
                     ggAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTa                                                                                                                                  <  2:2992516/150‑1 (MQ=255)
                     ggAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTa                                                                                                                                  <  1:532713/150‑1 (MQ=255)
                     ggAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTa                                                                                                                                  <  1:443519/150‑1 (MQ=255)
                       agaTTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTAcc                                                                                                                                >  2:456965/1‑150 (MQ=255)
                         aTTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCg                                                                                                                              >  1:1174687/1‑150 (MQ=255)
                         aTTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCg                                                                                                                              <  2:1637/150‑1 (MQ=255)
                            cccTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGgat                                                                                                                           >  1:3550698/1‑150 (MQ=255)
                            cccTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGgat                                                                                                                           >  2:275776/1‑150 (MQ=255)
                             ccTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGgatg                                                                                                                          >  1:1771612/1‑150 (MQ=255)
                                tcttATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGatgatg                                                                                                                       <  2:2636015/150‑1 (MQ=255)
                                tcttATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGatgatg                                                                                                                       >  2:2240621/1‑150 (MQ=255)
                                   tATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACt                                                                                                                    >  2:2971697/1‑150 (MQ=255)
                                       tGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATc                                                                                                                >  1:1132679/1‑150 (MQ=255)
                                       tGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATc                                                                                                                >  1:1505259/1‑150 (MQ=255)
                                        ggCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATcc                                                                                                               >  1:2721705/1‑150 (MQ=255)
                                         gCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATccc                                                                                                              >  1:987854/1‑150 (MQ=255)
                                              aGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggt                                                                                                         >  1:3262865/1‑150 (MQ=255)
                                               gCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggtg                                                                                                        >  2:2078770/1‑150 (MQ=255)
                                               gCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCggtg                                                                                                        <  1:2058657/150‑1 (MQ=255)
                                                  aTGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGgtggtg                                                                                                     >  1:3841449/1‑150 (MQ=255)
                                                   tGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGa                                                                                                    <  2:3589502/150‑1 (MQ=255)
                                                    gACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGAt                                                                                                   >  2:2153408/1‑150 (MQ=255)
                                                     aCGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATg                                                                                                  <  2:1669111/150‑1 (MQ=255)
                                                      cGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATgg                                                                                                 >  1:1319811/1‑150 (MQ=255)
                                                       gCAACGACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGc                                                                                                <  2:322910/150‑1 (MQ=255)
                                                               cAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACgg                                                                                        >  2:1619663/1‑150 (MQ=255)
                                                                 gAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTc                                                                                      >  1:3171975/1‑150 (MQ=255)
                                                                 gAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTc                                                                                      <  2:2481375/150‑1 (MQ=255)
                                                                  aaCTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCg                                                                                     <  1:3434489/150‑1 (MQ=255)
                                                                           tGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTcc                                                                            <  1:2537127/150‑1 (MQ=255)
                                                                            gACATTATCCCGACAAGAAATCCGCAAAATGATTCGTCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTAGCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCct                                                                           <  2:3977727/150‑1 (MQ=255)
                                                                               aTTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCctctc                                                                        >  1:3313889/1‑150 (MQ=255)
                                                                                 tATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCtt                                                                      <  2:1224359/150‑1 (MQ=255)
                                                                                    cccGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTg                                                                   <  2:3536611/150‑1 (MQ=255)
                                                                                      cGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGAt                                                                 >  2:3726852/1‑150 (MQ=255)
                                                                                      cGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGAt                                                                 >  2:3726828/1‑150 (MQ=255)
                                                                                        aCAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATgg                                                               >  2:4086971/1‑150 (MQ=255)
                                                                                        aCAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATgg                                                               <  2:2441573/150‑1 (MQ=255)
                                                                                          aaGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCg                                                             >  2:1434184/1‑150 (MQ=255)
                                                                                           aGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATGCCTCTCTTTTGATGGCGa                                                            >  2:1097134/1‑150 (MQ=255)
                                                                                            gAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGaa                                                           <  2:1882808/150‑1 (MQ=255)
                                                                                             aaaTCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTAGTCCTCTCTTTTGATGGCGAAc                                                          >  2:4116876/1‑150 (MQ=255)
                                                                                               aTCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTc                                                        >  1:2496820/1‑150 (MQ=255)
                                                                                               aTCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTc                                                        <  1:3336976/150‑1 (MQ=255)
                                                                                                   gCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGaca                                                    <  2:3349351/150‑1 (MQ=255)
                                                                                                   gCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTATTTTGAAGGCGAAATCGaca                                                    >  1:1125008/1‑150 (MQ=255)
                                                                                                    cAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGacac                                                   >  2:1715644/1‑150 (MQ=255)
                                                                                                     aaaaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCGACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACAcc                                                  <  2:2920880/150‑1 (MQ=255)
                                                                                                       aaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCTTTTGATGGCGAACTCGACAcccc                                                >  1:3500024/1‑149 (MQ=255)
                                                                                                             tCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGTACTCGACACCCAGCCACt                                          >  2:3276037/1‑150 (MQ=255)
                                                                                                               ggCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCa                                        <  1:580725/150‑1 (MQ=255)
                                                                                                                  aaCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAg                                     >  1:3788233/1‑150 (MQ=255)
                                                                                                                  aaCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAg                                     >  2:1166357/1‑150 (MQ=255)
                                                                                                                  aaCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAg                                     >  2:3165168/1‑150 (MQ=255)
                                                                                                                    cgtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCGCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaa                                   >  1:2511257/1‑150 (MQ=255)
                                                                                                                     gtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCAAAGaac                                  >  1:768773/1‑150 (MQ=255)
                                                                                                                     gtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaac                                  >  1:859737/1‑150 (MQ=255)
                                                                                                                      tcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaaca                                 <  2:3006006/150‑1 (MQ=255)
                                                                                                                       cgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCACTCTTTTGATGGCGAACTCGACACCCAGCCACACATAGaacaa                                >  1:3939878/1‑150 (MQ=255)
                                                                                                                        gtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaacaac                               >  1:1683076/1‑150 (MQ=255)
                                                                                                                           gtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAActct                            >  2:1750479/1‑150 (MQ=255)
                                                                                                                            tGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTg                           <  2:1305875/150‑1 (MQ=255)
                                                                                                                              cGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGgc                         <  1:3691276/150‑1 (MQ=255)
                                                                                                                               gTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGgcg                        >  1:1231491/1‑150 (MQ=255)
                                                                                                                                  aaCGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTTATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGgcgcgc                     >  1:300425/1‑150 (MQ=255)
                                                                                                                                    cGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCg                   <  1:1852675/150‑1 (MQ=255)
                                                                                                                                     gCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCgg                  <  2:653324/150‑1 (MQ=255)
                                                                                                                                      ccGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGt                 <  2:1922939/150‑1 (MQ=255)
                                                                                                                                       cGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTa                <  1:2971697/150‑1 (MQ=255)
                                                                                                                                        ggAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTaa               <  1:3094534/150‑1 (MQ=255)
                                                                                                                                         gAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAg              <  2:3914957/150‑1 (MQ=255)
                                                                                                                                         gAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACACATAGAACAACTCTGGCGCGCCGGTAAg              >  1:138896/1‑150 (MQ=255)
                                                                                                                                          aacaacAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgc             <  1:904563/150‑1 (MQ=255)
                                                                                                                                          aacaacAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgc             <  1:2381690/150‑1 (MQ=255)
                                                                                                                                           acaactGCAGGAAATGGGTCTACAAGGCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTGGCTGTTTTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgcg            <  1:2603444/150‑1 (MQ=255)
                                                                                                                                               cagcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGta        <  2:3960748/150‑1 (MQ=255)
                                                                                                                                               cagcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGta        >  1:2321626/1‑150 (MQ=255)
                                                                                                                                               cagcagAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGta        <  2:1800032/150‑1 (MQ=255)
                                                                                                                                               cagcagAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCGGGCGCGCCGGTACGCGCGta        >  2:385524/1‑150 (MQ=255)
                                                                                                                                                agcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGtat       >  1:1036781/1‑150 (MQ=255)
                                                                                                                                                agcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAAAGCGACACCCAGCCACTAAAAGAACAACACTGGCGCGCCGGTAATCGAGtaa       >  1:1036506/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACAGCCTGCGGCGGCATCTCGGAGATTCCCTTCTTATCTGGCACCAGCCATGACGCAACTACCAGAACTCCCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTT  >  CP009273/3049496‑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: