Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,818,010 C→A 66.3% H46N (CAT→AAT)  trmH → tRNA mG18‑2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,818,0100CA66.3% 319.6 / inf 331H46N (CAT→AAT) trmHtRNA mG18‑2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (55/56);  new base A (123/97);  total (178/153)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

Only 100 of 331 total aligned reads displayed.

AAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGC  >  CP009273/3817862‑3818158
                                                                                                                                                    |                                                                                                                                                    
 aaTTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAAt                                                                                                                                                    <  1:3787707/150‑1 (MQ=255)
    tAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGaa                                                                                                                                                 >  1:982114/1‑150 (MQ=255)
    tAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGaa                                                                                                                                                 >  1:31665/1‑150 (MQ=255)
    tAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGaa                                                                                                                                                 >  1:3648207/1‑150 (MQ=255)
    tAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGaa                                                                                                                                                 <  1:4263309/150‑1 (MQ=255)
     aaTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAg                                                                                                                                                >  2:1365444/1‑150 (MQ=255)
         ttttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCa                                                                                                                                            <  2:3475963/150‑1 (MQ=255)
            tATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGc                                                                                                                                         >  1:4051755/1‑150 (MQ=255)
             aTGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGcc                                                                                                                                        <  1:698918/150‑1 (MQ=255)
             aTGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGcc                                                                                                                                        <  1:997055/150‑1 (MQ=255)
               gAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGt                                                                                                                                      >  2:4739088/1‑150 (MQ=255)
               gAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGt                                                                                                                                      >  1:2921199/1‑150 (MQ=255)
                 aCCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGAGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTACTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTGCACGCCGTCt                                                                                                                                    >  2:3967429/1‑150 (MQ=255)
                   ccAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTgg                                                                                                                                  <  2:1185375/150‑1 (MQ=255)
                     aaCACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGcc                                                                                                                                >  1:1791085/1‑150 (MQ=255)
                     aaCACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGcc                                                                                                                                >  1:2489996/1‑150 (MQ=255)
                      acacGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCt                                                                                                                               <  2:801597/150‑1 (MQ=255)
                      acacGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCt                                                                                                                               <  1:214288/150‑1 (MQ=255)
                      acacGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCt                                                                                                                               <  1:5101019/150‑1 (MQ=255)
                       cacGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTg                                                                                                                              <  2:1404244/150‑1 (MQ=255)
                        acGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTgg                                                                                                                             <  2:4266612/150‑1 (MQ=255)
                          gTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTa                                                                                                                           >  2:3699451/1‑150 (MQ=255)
                          gTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTa                                                                                                                           >  2:5021299/1‑150 (MQ=255)
                           ttATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAg                                                                                                                          >  2:2550673/1‑150 (MQ=255)
                           ttATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAg                                                                                                                          >  2:4477213/1‑150 (MQ=255)
                           ttATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAg                                                                                                                          >  1:1705065/1‑150 (MQ=255)
                              tGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTCGCCg                                                                                                                       >  2:4861602/1‑150 (MQ=255)
                              tGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCg                                                                                                                       >  2:2038974/1‑150 (MQ=255)
                              tGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCAAAACGGTTCTGCGAGTAATCTTACCCCAGATGCCGATGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGACg                                                                                                                       >  2:4643230/1‑150 (MQ=255)
                                cACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCa                                                                                                                     >  1:4859223/1‑150 (MQ=255)
                                 aCGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCAt                                                                                                                    >  2:324776/1‑150 (MQ=255)
                                   gCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATgc                                                                                                                  <  1:2858947/150‑1 (MQ=255)
                                   gCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATgc                                                                                                                  >  2:1062670/1‑150 (MQ=255)
                                    cATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATgcg                                                                                                                 >  1:544220/1‑150 (MQ=255)
                                     aTCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATgcgc                                                                                                                >  1:1621389/1‑150 (MQ=255)
                                      tCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCa                                                                                                               >  1:4301089/1‑150 (MQ=255)
                                           gAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCCTAACGTTTCTGCGATTTTTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATg                                                                                                          <  2:416593/150‑1 (MQ=255)
                                           gAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATg                                                                                                          <  2:415861/150‑1 (MQ=255)
                                            aaaTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATCACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATgg                                                                                                         <  2:959089/150‑1 (MQ=255)
                                                gCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTc                                                                                                     <  1:4713605/150‑1 (MQ=255)
                                                   cGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGc                                                                                                  <  1:4995996/150‑1 (MQ=255)
                                                    gCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCa                                                                                                 <  1:1597691/150‑1 (MQ=255)
                                                    gCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCa                                                                                                 <  1:3914956/150‑1 (MQ=255)
                                                      cAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgc                                                                                               >  2:4123743/1‑150 (MQ=255)
                                                          cggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgcggcg                                                                                           <  2:1465873/150‑1 (MQ=255)
                                                           ggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGcggcgg                                                                                          >  1:1008301/1‑150 (MQ=255)
                                                              aGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgta                                                                                       >  1:5095093/1‑150 (MQ=255)
                                                                 cTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgtagta                                                                                    >  2:2261366/1‑150 (MQ=255)
                                                                  tGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCATCGGCGGGTAGTaa                                                                                   <  2:158998/150‑1 (MQ=255)
                                                                      cTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGc                                                                               >  2:2197822/1‑150 (MQ=255)
                                                                      cTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGc                                                                               >  2:1278208/1‑150 (MQ=255)
                                                                      cTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGc                                                                               >  2:2298347/1‑150 (MQ=255)
                                                                        gACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTg                                                                             <  1:3104551/150‑1 (MQ=255)
                                                                        gACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTg                                                                             >  2:4936683/1‑150 (MQ=255)
                                                                          ccGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTggg                                                                           >  1:1461567/1‑150 (MQ=255)
                                                                              cTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACa                                                                       <  2:941484/150‑1 (MQ=255)
                                                                                 cATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGt                                                                    >  2:3392245/1‑150 (MQ=255)
                                                                                      aGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGaaaa                                                               <  2:3254184/150‑1 (MQ=255)
                                                                                      aGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGaaaa                                                               >  2:2573647/1‑150 (MQ=255)
                                                                                        cAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAaca                                                             >  2:2993793/1‑150 (MQ=255)
                                                                                         aGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAacac                                                            <  1:2975030/150‑1 (MQ=255)
                                                                                         aGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAacac                                                            >  2:61722/1‑150 (MQ=255)
                                                                                            tCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACAcc                                                         >  1:855486/1‑150 (MQ=255)
                                                                                             ccACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCg                                                        >  2:5103946/1‑150 (MQ=255)
                                                                                             ccACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCg                                                        <  1:1319846/150‑1 (MQ=255)
                                                                                              cacaAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGc                                                       >  1:3287532/1‑150 (MQ=255)
                                                                                                caAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCAc                                                     <  1:244518/150‑1 (MQ=255)
                                                                                                   aCCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCAt                                                  >  1:2784249/1‑150 (MQ=255)
                                                                                                     cTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTg                                                <  1:459891/150‑1 (MQ=255)
                                                                                                      tCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTgg                                               >  1:3469150/1‑150 (MQ=255)
                                                                                                      tCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCAGTgg                                               >  2:3097398/1‑150 (MQ=255)
                                                                                                      tCATAAAGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGACGCATGCGCACCATGGATTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTgg                                               <  1:2533067/150‑1 (MQ=255)
                                                                                                        aTAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCg                                             >  1:2713225/1‑150 (MQ=255)
                                                                                                          aaCGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGAt                                           >  2:3631480/1‑150 (MQ=255)
                                                                                                            cGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGc                                         <  1:2344917/150‑1 (MQ=255)
                                                                                                               ttCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGt                                      >  2:3047228/1‑150 (MQ=255)
                                                                                                                tCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTc                                     <  1:4477213/150‑1 (MQ=255)
                                                                                                                  tGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGc                                   >  2:1216929/1‑150 (MQ=255)
                                                                                                                        tattCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAtct                             <  1:5205533/150‑1 (MQ=255)
                                                                                                                         attCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAtctc                            <  1:2109583/150‑1 (MQ=255)
                                                                                                                           tCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCaa                          >  2:1096048/1‑150 (MQ=255)
                                                                                                                            cGTACCGCAGATTCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCtaa                         >  1:4900974/1‑150 (MQ=255)
                                                                                                                              tACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAgg                       >  1:4108083/1‑150 (MQ=255)
                                                                                                                              tACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAgg                       >  2:4963697/1‑150 (MQ=255)
                                                                                                                                ccGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGcc                     <  1:1904100/150‑1 (MQ=255)
                                                                                                                                ccGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGcc                     <  2:5197068/150‑1 (MQ=255)
                                                                                                                                ccGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGcc                     <  1:4400849/150‑1 (MQ=255)
                                                                                                                                  gCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAg                   >  2:2414307/1‑150 (MQ=255)
                                                                                                                                        gCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATg             <  2:4788423/150‑1 (MQ=255)
                                                                                                                                        gCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATg             <  2:956251/150‑1 (MQ=255)
                                                                                                                                           gTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAg          <  1:2722041/150‑1 (MQ=255)
                                                                                                                                           gTTGGCGTACATGAAGGTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAg          <  1:407373/150‑1 (MQ=255)
                                                                                                                                           gTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAg          >  1:1868363/1‑150 (MQ=255)
                                                                                                                                             tGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCTCCATTGGCGATGCCGTCGCTCATCTCTAAGTCCAGGGCCTGCAGAt        >  2:2821447/1‑150 (MQ=255)
                                                                                                                                             tGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGAt        >  2:4016400/1‑150 (MQ=255)
                                                                                                                                             tGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGAt        >  2:1726196/1‑150 (MQ=255)
                                                                                                                                               gCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTc      <  1:1536467/150‑1 (MQ=255)
                                                                                                                                               gCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTc      <  1:1241352/150‑1 (MQ=255)
                                                                                                                                                 gTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTg    <  2:577216/150‑1 (MQ=255)
                                                                                                                                                   aCATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGc  <  2:19215/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGC  >  CP009273/3817862‑3818158

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: