Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,943,414 C→A 7.4% D456Y (GAT→TAT)  aspS ← aspartyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,943,4140CA 7.4% 1094.8 / 41.6 285D456Y (GAT→TAT) aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (147/117);  new base A (9/12);  total (156/129)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 290 total aligned reads displayed.

GATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGGCGCAGTGCACCCA  >  CP009273/1943265‑1943563
                                                                                                                                                     |                                                                                                                                                     
 aTACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCt                                                                                                                                                      <  2:1896920/150‑1 (MQ=255)
 aTACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCt                                                                                                                                                      >  1:2491724/1‑150 (MQ=255)
  tACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCtt                                                                                                                                                     >  2:555176/1‑150 (MQ=255)
   aCCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCttt                                                                                                                                                    >  2:86758/1‑150 (MQ=255)
     ccAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATATTTCg                                                                                                                                                  <  1:3587859/150‑1 (MQ=255)
      cAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCgg                                                                                                                                                 >  1:1890786/1‑150 (MQ=255)
       aGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCggt                                                                                                                                                >  2:2318366/1‑150 (MQ=255)
         aaTACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCggtga                                                                                                                                              >  1:1438907/1‑150 (MQ=255)
          aTACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCggtgag                                                                                                                                             >  1:3786619/1‑150 (MQ=255)
              cAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCggtgaggtga                                                                                                                                         >  2:3003070/1‑150 (MQ=255)
              cAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCggtgaggtga                                                                                                                                         >  2:406903/1‑150 (MQ=255)
                  cacCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACgg                                                                                                                                     >  1:3450037/1‑150 (MQ=255)
                  cacCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACgg                                                                                                                                     >  1:4228897/1‑150 (MQ=255)
                  cacCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACgg                                                                                                                                     >  1:1643359/1‑150 (MQ=255)
                    ccGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGAt                                                                                                                                   >  2:1087313/1‑150 (MQ=255)
                    ccGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGAt                                                                                                                                   <  1:4921839/150‑1 (MQ=255)
                    ccGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATATTTCGGTGAGGTGAACGGAt                                                                                                                                   <  2:1658699/150‑1 (MQ=255)
                     cGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATg                                                                                                                                  >  2:1295711/1‑150 (MQ=255)
                     cGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATg                                                                                                                                  >  2:3649990/1‑150 (MQ=255)
                     cGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATg                                                                                                                                  >  1:4038110/1‑150 (MQ=255)
                         tgctgcATATCACCATTATGGATACGTACTGAACCACCGCACACATCTTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTACCGCTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGc                                                                                                                              >  1:3624529/1‑150 (MQ=255)
                           ctgcATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCAt                                                                                                                            <  2:2390812/150‑1 (MQ=255)
                           ctgcATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCAt                                                                                                                            <  2:4689693/150‑1 (MQ=255)
                            tgcATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCAtt                                                                                                                           >  1:1745019/1‑150 (MQ=255)
                             gcATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTg                                                                                                                          >  1:3912649/1‑150 (MQ=255)
                              cATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGc                                                                                                                         >  2:1219644/1‑150 (MQ=255)
                              cATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGc                                                                                                                         <  2:2288973/150‑1 (MQ=255)
                                tatCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCg                                                                                                                       >  2:3789212/1‑150 (MQ=255)
                                 atCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGt                                                                                                                      >  2:1396768/1‑150 (MQ=255)
                                    aCCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAg                                                                                                                   >  1:2921684/1‑150 (MQ=255)
                                       aTTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGgcc                                                                                                                <  2:296961/150‑1 (MQ=255)
                                       aTTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGgcc                                                                                                                <  1:4811310/150‑1 (MQ=255)
                                        ttATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGgccg                                                                                                               <  2:1991325/150‑1 (MQ=255)
                                         tATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGgccgc                                                                                                              >  1:3345820/1‑150 (MQ=255)
                                           tGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATCTCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCt                                                                                                            <  1:1242831/150‑1 (MQ=255)
                                           tGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATATTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCt                                                                                                            >  1:905979/1‑150 (MQ=255)
                                             gATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATATTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTc                                                                                                          >  2:1588684/1‑150 (MQ=255)
                                                 cGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACcg                                                                                                      >  2:4656782/1‑150 (MQ=255)
                                                   tACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACcgtc                                                                                                    <  1:1318450/150‑1 (MQ=255)
                                                   tACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACcgtc                                                                                                    <  1:305484/150‑1 (MQ=255)
                                                     cTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATATTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACcgtcgt                                                                                                  <  1:2847912/150‑1 (MQ=255)
                                                      tGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCgtcgtc                                                                                                 <  2:736589/150‑1 (MQ=255)
                                                       gAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCt                                                                                                >  1:5082853/1‑150 (MQ=255)
                                                        aaCCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCtt                                                                                               >  2:2641372/1‑150 (MQ=255)
                                                         accaccGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTc                                                                                              <  2:318587/150‑1 (MQ=255)
                                                         accaccGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTc                                                                                              <  1:1949645/150‑1 (MQ=255)
                                                          ccaccGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCa                                                                                             <  2:2558917/150‑1 (MQ=255)
                                                               gCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCCAACAt                                                                                        >  2:17207/1‑150 (MQ=255)
                                                                cccACTTCGTAACCATTGATGACCAGATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATc                                                                                       <  1:3712290/150‑1 (MQ=255)
                                                                   aCTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCggg                                                                                    >  2:492182/1‑150 (MQ=255)
                                                                   aCTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCggg                                                                                    >  2:2066000/1‑150 (MQ=255)
                                                                        gTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTc                                                                               >  1:2261818/1‑150 (MQ=255)
                                                                           aCCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGAt                                                                            <  1:4875922/150‑1 (MQ=255)
                                                                           aCCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGAt                                                                            >  2:3112836/1‑150 (MQ=255)
                                                                            ccATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATa                                                                           <  2:4639288/150‑1 (MQ=255)
                                                                            ccATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATa                                                                           <  1:4327434/150‑1 (MQ=255)
                                                                                tgatgaCCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAAccc                                                                       <  2:3510200/150‑1 (MQ=255)
                                                                                 gatgaCCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCca                                                                      <  2:2796332/150‑1 (MQ=255)
                                                                                     aCCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGc                                                                  <  2:1009744/150‑1 (MQ=255)
                                                                                      ccATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCg                                                                 <  2:3856000/150‑1 (MQ=255)
                                                                                          atCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGc                                                             >  1:2673277/1‑150 (MQ=255)
                                                                                           tCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGcc                                                            >  2:4186316/1‑150 (MQ=255)
                                                                                              tAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCAt                                                         >  2:2708000/1‑150 (MQ=255)
                                                                                                    ttCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCtt                                                   <  1:2318366/150‑1 (MQ=255)
                                                                                                          aCCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTtcgtcg                                             <  1:2966153/150‑1 (MQ=255)
                                                                                                              cATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCa                                         <  1:1559146/150‑1 (MQ=255)
                                                                                                              cATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCa                                         <  1:4370083/150‑1 (MQ=255)
                                                                                                              cATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCa                                         >  2:3996127/1‑150 (MQ=255)
                                                                                                                ttttCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGa                                       <  2:2501889/150‑1 (MQ=255)
                                                                                                                 tttCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGAc                                      >  2:3599497/1‑150 (MQ=255)
                                                                                                                 tatCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTACAAACATCGGGAAGTCGATAACCCACAGCCGTGCCCATTAGCTTTCGTCGGTCAGAc                                      >  1:4385708/3‑150 (MQ=255)
                                                                                                                     cGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAg                                  <  2:497616/150‑1 (MQ=255)
                                                                                                                     cGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAg                                  >  1:2555697/1‑150 (MQ=255)
                                                                                                                     cGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAg                                  <  2:497315/150‑1 (MQ=255)
                                                                                                                     cGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAg                                  >  1:2044675/1‑150 (MQ=255)
                                                                                                                        tGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTc                               <  1:1430262/150‑1 (MQ=255)
                                                                                                                         gCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGATCGGATGGTGCATTGCCGTCAGGCCGCCTTCAACGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCACATTTGCTTTCGTCGGTCAGACCAAGGTCt                              <  2:3268086/150‑1 (MQ=255)
                                                                                                                          cAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCtt                             >  2:3344490/1‑150 (MQ=255)
                                                                                                                           aGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGACGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCttt                            >  1:2802958/1‑150 (MQ=255)
                                                                                                                           aGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCttt                            >  2:3021775/1‑150 (MQ=255)
                                                                                                                             cTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTAc                          <  1:2355653/150‑1 (MQ=255)
                                                                                                                              tttCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTAcc                         >  1:2867402/1‑150 (MQ=255)
                                                                                                                               ttCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTAccc                        <  1:1386529/150‑1 (MQ=255)
                                                                                                                                 cAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCAc                      >  2:1500991/1‑150 (MQ=255)
                                                                                                                                  aGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACt                     >  1:4070222/1‑150 (MQ=255)
                                                                                                                                  aGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACt                     >  1:4638482/1‑150 (MQ=255)
                                                                                                                                  aGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACt                     >  1:1667893/1‑150 (MQ=255)
                                                                                                                                   gTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACtt                    <  2:4778604/150‑1 (MQ=255)
                                                                                                                                     tCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTc                  <  2:4453667/150‑1 (MQ=255)
                                                                                                                                        gCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAgg               <  1:3649926/150‑1 (MQ=255)
                                                                                                                                         cAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGgc              >  1:1035315/1‑150 (MQ=255)
                                                                                                                                          aGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGgcg             >  2:4163793/1‑150 (MQ=255)
                                                                                                                                           gCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGgcgc            >  1:1806142/1‑150 (MQ=255)
                                                                                                                                             cGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGGCGCAg          >  2:1148784/1‑150 (MQ=255)
                                                                                                                                              gTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGGCGCAGt         <  2:3912649/150‑1 (MQ=255)
                                                                                                                                              gTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGGCGCAGt         <  1:3582404/150‑1 (MQ=255)
                                                                                                                                              gTCATATATTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGGCGCAGt         >  1:1145828/1‑150 (MQ=255)
                                                                                                                                               tCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGGCGCAGTg        <  1:2190924/150‑1 (MQ=255)
                                                                                                                                                 atatCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGGCGCAGTGCa      <  1:849397/150‑1 (MQ=255)
                                                                                                                                                     cTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGGCGCAGTGCACCCa  <  2:2414239/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGTCTTCAAACATCGGGAAGTCGATAACCCACAGCGGTGCCCATTTGCTTTCGTCGGTCAGACCAAGGTCTTTACCCACTTTCAGGCGCAGTGCACCCA  >  CP009273/1943265‑1943563

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: