Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,880,857 G→T 16.3% M252I (ATG→ATT yeaX → putative YeaWX dioxygenase beta subunit, reductase component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,880,8570GT16.3% 863.8 / 99.4 268M252I (ATG→ATTyeaXputative YeaWX dioxygenase beta subunit, reductase component
Reads supporting (aligned to +/- strand):  ref base G (108/116);  new base T (19/25);  total (127/141)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

ATTGAAGCGGTAAGAAGTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAAGAGCGTGCCAGCCA  >  CP009273/1880708‑1881006
                                                                                                                                                     |                                                                                                                                                     
aTTGAAGCGGTAAGAAGTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATg                                                                                                                                                       <  1:952131/150‑1 (MQ=255)
aTTGAAGCGGTAAGAAGTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTGGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATg                                                                                                                                                       >  2:3661021/1‑150 (MQ=255)
 ttGAAGCGGTAAGAAGTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGa                                                                                                                                                      <  1:39478/150‑1 (MQ=255)
     aGCGGTAAGAAGTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTAt                                                                                                                                                  <  2:3111048/150‑1 (MQ=255)
      gCGGTAAGAAGTGAAGCTTCGCGTCTGGACATCGCCGCCGATACGCTTCACTTTGAGCAATTTTCTATCTACTACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATTACTAtt                                                                                                                                                 <  1:1376413/150‑1 (MQ=255)
      gCGGTAAGAAGTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTAtt                                                                                                                                                 <  1:552350/150‑1 (MQ=255)
        ggTAAGAAGTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTAtttt                                                                                                                                               <  2:519628/150‑1 (MQ=255)
        ggTAAGAAGTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTAtttt                                                                                                                                               >  2:1165297/1‑150 (MQ=255)
         gTAAGAAGTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTg                                                                                                                                              <  1:733374/150‑1 (MQ=255)
          tAAGAAGTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGc                                                                                                                                             <  1:4263331/150‑1 (MQ=255)
              aagTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGt                                                                                                                                         >  2:1975706/1‑150 (MQ=255)
               agTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATTACTATTTTGCAGGtt                                                                                                                                        <  2:2615070/150‑1 (MQ=255)
                gTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTa                                                                                                                                       <  1:4597958/150‑1 (MQ=255)
                gTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTa                                                                                                                                       <  2:219145/150‑1 (MQ=255)
                 tGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTAt                                                                                                                                      <  2:941869/150‑1 (MQ=255)
                  gAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTAtt                                                                                                                                     <  2:2453423/150‑1 (MQ=255)
                  gAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTAtt                                                                                                                                     >  1:3340237/1‑150 (MQ=255)
                  gAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTAtt                                                                                                                                     <  1:3292781/150‑1 (MQ=255)
                  gAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTAtt                                                                                                                                     >  1:3159667/1‑150 (MQ=255)
                              cTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATaaa                                                                                                                         >  2:1478094/1‑150 (MQ=255)
                               tGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATTACTATTTTGCAGGTTATTGAAAATAATAAAg                                                                                                                        <  2:149424/150‑1 (MQ=255)
                                ggACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGc                                                                                                                       >  1:1997079/1‑150 (MQ=255)
                                ggACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGc                                                                                                                       >  2:3010730/1‑150 (MQ=255)
                                ggACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGc                                                                                                                       <  2:614501/150‑1 (MQ=255)
                                 gACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGcc                                                                                                                      <  2:2926808/150‑1 (MQ=255)
                                 gACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAAAAAAGcc                                                                                                                      >  2:1802548/1‑150 (MQ=255)
                                  aCATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCcg                                                                                                                     <  1:2355915/150‑1 (MQ=255)
                                        ccgccgATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAg                                                                                                               >  1:2685164/1‑150 (MQ=255)
                                        ccgccgATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAg                                                                                                               >  2:1707589/1‑150 (MQ=255)
                                           ccgATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTgg                                                                                                            >  2:1122683/1‑150 (MQ=255)
                                            cgATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATTACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGa                                                                                                           >  1:4626922/1‑150 (MQ=255)
                                              aTACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAAt                                                                                                         <  2:2400836/150‑1 (MQ=255)
                                              aTACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAAt                                                                                                         <  1:153171/150‑1 (MQ=255)
                                                  gCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGttt                                                                                                     <  1:371224/150‑1 (MQ=255)
                                                         tttGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGt                                                                                              >  1:2433034/1‑150 (MQ=255)
                                                         tttGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGt                                                                                              <  1:905576/150‑1 (MQ=255)
                                                          ttGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATTACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTg                                                                                             <  1:2598854/150‑1 (MQ=255)
                                                          ttGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTg                                                                                             <  1:429862/150‑1 (MQ=255)
                                                             aGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAAGGTTTATGTCGTGAAg                                                                                          >  1:2952142/1‑150 (MQ=255)
                                                               cAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATTACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAggg                                                                                        <  1:4009155/150‑1 (MQ=255)
                                                                aaTTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAgggg                                                                                       >  1:1476149/1‑150 (MQ=255)
                                                                  tttGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTa                                                                                     <  1:1826642/150‑1 (MQ=255)
                                                                    tGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATg                                                                                   >  2:744321/1‑150 (MQ=255)
                                                                    tGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATg                                                                                   >  2:1121973/1‑150 (MQ=255)
                                                                     gCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGc                                                                                  <  2:411025/150‑1 (MQ=255)
                                                                      cTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCg                                                                                 >  2:3532512/1‑150 (MQ=255)
                                                                       tATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAATGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCgg                                                                                >  1:179546/1‑150 (MQ=255)
                                                                        aTCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGa                                                                               <  1:1313848/150‑1 (MQ=255)
                                                                        aTCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGa                                                                               <  1:4190139/150‑1 (MQ=255)
                                                                        aTCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGa                                                                               <  2:3945898/150‑1 (MQ=255)
                                                                        aTCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGa                                                                               <  2:3907909/150‑1 (MQ=255)
                                                                             aGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTg                                                                          <  2:1236290/150‑1 (MQ=255)
                                                                             aGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTg                                                                          >  1:3334823/1‑150 (MQ=255)
                                                                              gACAAAACCGGCGATGCATTTACCCGGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTgt                                                                         <  2:4489581/150‑2 (MQ=255)
                                                                               aCAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCg                                                                        <  2:434181/150‑1 (MQ=255)
                                                                                cAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGa                                                                       >  2:3196683/1‑150 (MQ=255)
                                                                                 aaaaCCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATTACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGaa                                                                      >  1:48759/1‑150 (MQ=255)
                                                                                   aaCCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAAc                                                                    <  1:3291629/150‑1 (MQ=255)
                                                                                    aCCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATTACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACa                                                                   <  1:4639922/150‑1 (MQ=255)
                                                                                     ccGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAg                                                                  >  1:282268/1‑150 (MQ=255)
                                                                                      cGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGc                                                                 <  2:480070/150‑1 (MQ=255)
                                                                                       ggCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCa                                                                >  2:3902/1‑150 (MQ=255)
                                                                                       ggCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCa                                                                >  1:382116/1‑150 (MQ=255)
                                                                                           aTGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTCTTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATAc                                                            >  2:1372276/1‑150 (MQ=255)
                                                                                            tGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACt                                                           <  1:3452100/150‑1 (MQ=255)
                                                                                              cATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTgg                                                         >  2:3776277/1‑150 (MQ=255)
                                                                                                 ttACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAg                                                      <  1:4534346/150‑1 (MQ=255)
                                                                                                   aCCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGt                                                    <  1:2781919/150‑1 (MQ=255)
                                                                                                         gTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCt                                              <  2:3861174/150‑1 (MQ=255)
                                                                                                         gTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCt                                              >  2:2002667/1‑150 (MQ=255)
                                                                                                          tGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTg                                             >  1:1658192/1‑150 (MQ=255)
                                                                                                               gCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCAt                                        >  2:4426022/1‑150 (MQ=255)
                                                                                                                cccGTTCCGGATAAGAGTTTGTGGTGCCGGTAGAGATGACTATTTTGCATGTTATTGAATATTATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACATCAATACTGGAAGGTGAAGCTGACCATc                                       <  1:1757915/150‑1 (MQ=255)
                                                                                                                cccGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATc                                       <  2:3013605/150‑1 (MQ=255)
                                                                                                                 ccGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCg                                      <  2:234886/150‑1 (MQ=255)
                                                                                                                 ccGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCg                                      <  1:3925204/150‑1 (MQ=255)
                                                                                                                 ccGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCg                                      <  2:1920050/150‑1 (MQ=255)
                                                                                                                   gTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATTACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCggg                                    >  2:3482590/1‑150 (MQ=255)
                                                                                                                    ttCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGa                                   <  2:4606598/150‑1 (MQ=255)
                                                                                                                    ttCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGa                                   >  1:1040846/1‑150 (MQ=255)
                                                                                                                       cGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCa                                <  2:2516963/150‑1 (MQ=255)
                                                                                                                       cGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCa                                <  1:1009390/150‑1 (MQ=255)
                                                                                                                       cGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCa                                >  1:584660/1‑150 (MQ=255)
                                                                                                                        ggAAAAGAGTTTGTGGTGCCGGAAGAGATTACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCaa                               >  1:2408558/1‑150 (MQ=255)
                                                                                                                           aaaGAGTTTGTGGTGCCGGAAGAGATTACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAatat                            >  2:4084761/1‑150 (MQ=255)
                                                                                                                           aaaGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAatat                            >  1:3403711/1‑150 (MQ=255)
                                                                                                                           aaaGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAatat                            <  1:1358072/150‑1 (MQ=255)
                                                                                                                           aaaGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAatat                            <  1:138158/150‑1 (MQ=255)
                                                                                                                                       gtgCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGaa                <  2:1763201/150‑1 (MQ=255)
                                                                                                                                        tgCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAag               <  2:3469043/150‑1 (MQ=255)
                                                                                                                                         gCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAaga              >  2:3754969/1‑150 (MQ=255)
                                                                                                                                            ggAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAAGAGCg           <  1:1122683/150‑1 (MQ=255)
                                                                                                                                                  gaTTACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAAGAGCGTgccag     >  2:3085103/1‑150 (MQ=255)
                                                                                                                                                  gaTGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAAGAGCGTgccag     >  2:1414872/1‑150 (MQ=255)
                                                                                                                                                  gaTGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAAGAGCGTgccag     >  1:4588847/1‑150 (MQ=255)
                                                                                                                                                  gaTGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAAGAGCGTgccag     >  2:2278389/1‑150 (MQ=255)
                                                                                                                                                   aTGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAAGAGCGTgccagc    >  1:3128919/1‑150 (MQ=255)
                                                                                                                                                    tGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAAGAGCGTgccagcc   >  2:2671358/1‑150 (MQ=255)
                                                                                                                                                     gACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAAGAGCGTgccagcca  <  2:3340236/150‑1 (MQ=255)
                                                                                                                                                     gACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAAGAGCGTgccagcca  <  1:2165641/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTGAAGCGGTAAGAAGTGAAGCTGCGCGTCTGGACATCGCCGCCGATACGCTGCACTTTGAGCAATTTGCTATCGAAGACAAAACCGGCGATGCATTTACCCTGGTGCTTGCCCGTTCCGGAAAAGAGTTTGTGGTGCCGGAAGAGATGACTATTTTGCAGGTTATTGAAAATAATAAAGCCGCGAAAGTGGAATGTTTATGTCGTGAAGGGGTATGCGGAACCTGCGAAACAGCAATACTGGAAGGTGAAGCTGACCATCGGGATCAATATTTTAGCGATGAAGAGCGTGCCAGCCA  >  CP009273/1880708‑1881006

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: