Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,070,543 G→T 100% D57Y (GAT→TAT)  dtd → D‑tyr‑tRNA(Tyr) deacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,070,5430GT96.1% 542.4 / 12.2 180D57Y (GAT→TAT) dtdD‑tyr‑tRNA(Tyr) deacylase
Reads supporting (aligned to +/- strand):  ref base G (4/3);  new base T (91/82);  total (95/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 180 total aligned reads displayed.

ACCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGT  >  CP009273/4070399‑4070687
                                                                                                                                                |                                                                                                                                                
aCCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATg                                                                                                                                               >  2:3477331/1‑148 (MQ=255)
    gTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGa                                                                                                                                           >  1:715811/1‑148 (MQ=255)
    gTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGa                                                                                                                                           <  2:3246175/148‑1 (MQ=255)
     tGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGaa                                                                                                                                          <  2:2581915/148‑1 (MQ=255)
      gCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAg                                                                                                                                         >  2:1334595/1‑148 (MQ=255)
       ccAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCATACCGTCTGTGCGAGCGGGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAgg                                                                                                                                        <  1:235089/148‑1 (MQ=255)
       ccAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAgg                                                                                                                                        <  1:3308201/148‑1 (MQ=255)
          gCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCaa                                                                                                                                     >  1:481838/1‑148 (MQ=255)
             tCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGAt                                                                                                                                  >  1:2671130/1‑148 (MQ=255)
              cACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATg                                                                                                                                 >  1:2996633/1‑148 (MQ=255)
              cACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATg                                                                                                                                 <  1:376263/148‑1 (MQ=255)
               aCCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTTGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGa                                                                                                                                >  1:688993/1‑148 (MQ=255)
               aCCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGa                                                                                                                                <  2:1488823/148‑1 (MQ=255)
                  gTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAAtc                                                                                                                             <  2:233893/148‑1 (MQ=255)
                   tGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAAtct                                                                                                                            <  1:3139565/148‑1 (MQ=255)
                      agggagAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCaa                                                                                                                         <  2:2883095/148‑1 (MQ=255)
                         gagaAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGt                                                                                                                      >  1:1554500/1‑148 (MQ=255)
                          agaAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTg                                                                                                                     >  1:3480500/1‑148 (MQ=255)
                          agaAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTg                                                                                                                     <  1:1666417/148‑1 (MQ=255)
                            aaGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCa                                                                                                                   >  1:220302/1‑148 (MQ=255)
                             aGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCaa                                                                                                                  >  2:2117615/1‑148 (MQ=255)
                             aGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCaa                                                                                                                  >  1:3077090/1‑148 (MQ=255)
                                     gCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAggcggg                                                                                                          <  2:2540446/148‑1 (MQ=255)
                                        aaaTTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGgcgggcgg                                                                                                       <  1:3444133/148‑1 (MQ=255)
                                         aaTTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGggcggc                                                                                                      <  1:2955345/148‑1 (MQ=255)
                                         aaTTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGggcggc                                                                                                      >  1:1780474/1‑148 (MQ=255)
                                           ttGGCGCGGGACTTTTGGTGTTATTGGTTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAg                                                                                                    >  1:1496465/1‑148 (MQ=255)
                                           ttGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAg                                                                                                    >  2:505841/1‑148 (MQ=255)
                                           ttGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAg                                                                                                    >  2:3438739/1‑148 (MQ=255)
                                           ttGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAg                                                                                                    >  1:2580867/1‑148 (MQ=255)
                                           ttGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAg                                                                                                    <  2:3426604/148‑1 (MQ=255)
                                            tGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAgt                                                                                                   <  1:1835185/148‑1 (MQ=255)
                                             ggCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAgtg                                                                                                  >  2:3048341/1‑148 (MQ=255)
                                                 cgGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCt                                                                                              <  1:749109/148‑1 (MQ=255)
                                                   ggACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCtgg                                                                                            >  1:2409177/1‑148 (MQ=255)
                                                    gACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCtggt                                                                                           >  1:1086478/1‑148 (MQ=255)
                                                    gACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCtggt                                                                                           <  2:1850299/148‑1 (MQ=255)
                                                     aCTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCtggtg                                                                                          >  2:947766/1‑148 (MQ=255)
                                                         ttGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGttt                                                                                      <  2:1833830/148‑1 (MQ=255)
                                                         ttGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGttt                                                                                      <  2:2458188/148‑1 (MQ=255)
                                                         ttGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGttt                                                                                      <  2:577842/148‑1 (MQ=255)
                                                             tgtTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCa                                                                                  <  2:334740/148‑1 (MQ=255)
                                                             tgtTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCa                                                                                  >  2:1782492/1‑148 (MQ=255)
                                                                   tGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTAc                                                                            >  2:2055864/1‑148 (MQ=255)
                                                                   tGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTAc                                                                            <  1:1233111/148‑1 (MQ=255)
                                                                    gggTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTAcc                                                                           <  2:415024/148‑1 (MQ=255)
                                                                    gggTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGATACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTAcc                                                                           <  1:1186986/148‑1 (MQ=255)
                                                                     ggTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTAccc                                                                          <  1:2696338/148‑1 (MQ=255)
                                                                     ggTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTAccc                                                                          >  2:1574055/1‑148 (MQ=255)
                                                                      gtgtCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGGGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGTTGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCt                                                                         <  2:565973/148‑1 (MQ=255)
                                                                         tCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTcgc                                                                      <  2:2074761/148‑1 (MQ=255)
                                                                         tCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTcgc                                                                      <  2:1105919/148‑1 (MQ=255)
                                                                             aaaGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCa                                                                  >  1:3537152/1‑148 (MQ=255)
                                                                               aGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGa                                                                <  2:615761/148‑1 (MQ=255)
                                                                                ggATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGAt                                                               >  2:695696/1‑148 (MQ=255)
                                                                                 gATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATa                                                              >  2:2918739/1‑148 (MQ=255)
                                                                                 gATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATa                                                              <  1:2262204/148‑1 (MQ=255)
                                                                                  aTGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATAc                                                             >  1:3561289/1‑148 (MQ=255)
                                                                                    gacgacGATCAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCg                                                           >  2:797106/1‑148 (MQ=255)
                                                                                    gacgacGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCg                                                           >  2:3137124/1‑148 (MQ=255)
                                                                                      cgacgaACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCCGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGaa                                                         <  2:436019/148‑1 (MQ=255)
                                                                                       gacgaACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAAc                                                        <  1:2049237/148‑1 (MQ=255)
                                                                                        acgaACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACg                                                       >  1:707101/1‑148 (MQ=255)
                                                                                         cgaACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACgg                                                      <  1:1385619/148‑1 (MQ=255)
                                                                                         cgaACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACgg                                                      <  1:1382298/148‑1 (MQ=255)
                                                                                         cgaACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACgg                                                      <  1:3157268/148‑1 (MQ=255)
                                                                                           aaCAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACgggg                                                    <  1:471610/148‑1 (MQ=255)
                                                                                             cAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGa                                                  >  1:1483231/1‑148 (MQ=255)
                                                                                             cAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGa                                                  <  1:558133/148‑1 (MQ=255)
                                                                                              aGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGAt                                                 <  1:359678/148‑1 (MQ=255)
                                                                                               gAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATg                                                >  1:3039622/1‑148 (MQ=255)
                                                                                                 aaGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATgcg                                              <  1:804157/148‑1 (MQ=255)
                                                                                                   gCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGcc                                            <  2:2671130/148‑1 (MQ=255)
                                                                                                   gCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGcc                                            >  2:883796/1‑148 (MQ=255)
                                                                                                   gCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGcc                                            <  2:1518449/148‑1 (MQ=255)
                                                                                                     aaaCCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCa                                          >  2:787178/1‑148 (MQ=255)
                                                                                                      aaCCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCaa                                         <  2:1924464/148‑1 (MQ=255)
                                                                                                       aCCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAg                                        >  1:331315/1‑148 (MQ=255)
                                                                                                       aCCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAg                                        >  1:330660/1‑148 (MQ=255)
                                                                                                         cGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGtt                                      >  2:2254050/1‑148 (MQ=255)
                                                                                                         cGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGtt                                      >  1:2787733/1‑148 (MQ=255)
                                                                                                           tCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTtc                                    >  1:2599178/1‑148 (MQ=255)
                                                                                                                gCGAGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCaa                               <  2:1015341/148‑1 (MQ=255)
                                                                                                                   aGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGaa                                                         <  2:845991/119‑1 (MQ=255)
                                                                                                                   aGCGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGaa                                                         >  1:845991/1‑119 (MQ=255)
                                                                                                                     cGTGTGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTg                          >  2:1736056/1‑148 (MQ=255)
                                                                                                                      gtgtGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGc                         <  1:2348406/148‑1 (MQ=255)
                                                                                                                      gtgtGCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGc                         <  1:3015818/148‑1 (MQ=255)
                                                                                                                        gtgCTCGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCAt                       >  1:1648976/1‑148 (MQ=255)
                                                                                                                             cGGCTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCg                  <  2:481838/148‑1 (MQ=255)
                                                                                                                                cTACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGAt               >  1:2313725/1‑148 (MQ=255)
                                                                                                                                 tACCGCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATc              >  2:879265/1‑148 (MQ=255)
                                                                                                                                     gCATCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCgcgc          >  1:2842322/1‑148 (MQ=255)
                                                                                                                                       aTCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCag        >  1:3297671/1‑148 (MQ=255)
                                                                                                                                        tCTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCaga       <  2:3561289/148‑1 (MQ=255)
                                                                                                                                         cTTTAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCagag      <  2:3010419/148‑1 (MQ=255)
                                                                                                                                           ttAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGc    >  2:1753099/1‑148 (MQ=255)
                                                                                                                                            tAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCg   <  2:2038495/148‑1 (MQ=255)
                                                                                                                                            tAGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCg   >  1:2641208/1‑148 (MQ=255)
                                                                                                                                             aGCTATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGt  <  1:947992/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ACCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTCGCCGCAGATACCGAACGGGGGATGCGCCCAAGTTTCTCCAAAGGTGCATCACCGGATCGCGCAGAGGCGT  >  CP009273/4070399‑4070687

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: