Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,602,131 G→T 15.5% D177Y (GAT→TAT)  tam → trans‑aconitate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,602,1310GT15.5% 891.2 / 99.7 284D177Y (GAT→TAT) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (110/130);  new base T (16/28);  total (126/158)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

AGTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTGGA  >  CP009273/1601983‑1602279
                                                                                                                                                    |                                                                                                                                                    
aGTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGa                                                                                                                                                     >  1:3416131/1‑150 (MQ=255)
aGTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGa                                                                                                                                                     >  2:4140156/1‑150 (MQ=255)
 gTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGat                                                                                                                                                    >  2:2945574/1‑150 (MQ=255)
 gTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGat                                                                                                                                                    <  1:3529437/150‑1 (MQ=255)
     aGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGTGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCt                                                                                                                                                >  2:4395265/1‑150 (MQ=255)
       aTGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTgg                                                                                                                                              >  2:2585894/1‑150 (MQ=255)
        tGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGc                                                                                                                                             <  2:1200776/150‑1 (MQ=255)
         gCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCg                                                                                                                                            >  1:1821051/1‑150 (MQ=255)
          ccAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGa                                                                                                                                           >  2:3380461/1‑150 (MQ=255)
          ccAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGa                                                                                                                                           <  1:4499185/150‑1 (MQ=255)
              aTAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAacca                                                                                                                                       >  1:2045300/1‑150 (MQ=255)
                aaCTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAaccacc                                                                                                                                     <  1:2205198/150‑1 (MQ=255)
                aaCTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAaccacc                                                                                                                                     >  2:4270714/1‑150 (MQ=255)
                   tggctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACctac                                                                                                                                  <  1:4188964/150‑1 (MQ=255)
                   tggctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACctac                                                                                                                                  >  1:3060679/1‑150 (MQ=255)
                      ctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACctactct                                                                                                                               >  2:4179422/1‑148 (MQ=255)
                      ctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTAt                                                                                                                               >  1:2896384/1‑150 (MQ=255)
                      ctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTAt                                                                                                                               >  1:3511183/1‑150 (MQ=255)
                       tggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATc                                                                                                                              >  2:342225/1‑150 (MQ=255)
                        ggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCa                                                                                                                             >  1:1039231/1‑150 (MQ=255)
                         gAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCAc                                                                                                                            <  2:783842/150‑1 (MQ=255)
                            ccGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAg                                                                                                                         <  2:564988/150‑1 (MQ=255)
                             cgtcccATGTGCTCATGGGCGTAGTTGCCTGGGTACAAAACTACCCGGTTCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGTTATTTTGTGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGa                                                                                                                        <  2:1312548/147‑1 (MQ=255)
                             cGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGa                                                                                                                        <  1:592525/150‑1 (MQ=255)
                               aCCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATg                                                                                                                      <  2:2636338/150‑1 (MQ=255)
                                cccATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACGACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGc                                                                                                                     <  2:391151/150‑1 (MQ=255)
                                 ccATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGcc                                                                                                                    <  1:3089582/150‑1 (MQ=255)
                                 ccATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGcc                                                                                                                    >  2:3557683/1‑150 (MQ=255)
                                    tgtgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTc                                                                                                                 <  2:1086458/150‑1 (MQ=255)
                                    tgtgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTc                                                                                                                 <  1:769264/150‑1 (MQ=255)
                                     gtgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTca                                                                                                                <  1:1746962/150‑1 (MQ=255)
                                     gtgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTca                                                                                                                >  1:1514492/1‑150 (MQ=255)
                                       gCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTcaca                                                                                                              <  2:2339019/150‑1 (MQ=255)
                                       gCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTcaca                                                                                                              <  2:4464247/150‑1 (MQ=255)
                                        cTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCacac                                                                                                             <  2:795997/150‑1 (MQ=255)
                                         tCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACAcc                                                                                                            >  1:3587178/1‑150 (MQ=255)
                                          cATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCa                                                                                                           <  1:3772379/150‑1 (MQ=255)
                                                   aGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTAt                                                                                                  <  2:1915951/150‑1 (MQ=255)
                                                   aGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTAt                                                                                                  <  1:2781260/150‑1 (MQ=255)
                                                    gTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATc                                                                                                 <  1:4134946/150‑1 (MQ=255)
                                                    gTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATc                                                                                                 >  1:642997/1‑150 (MQ=255)
                                                    gTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATc                                                                                                 >  1:1242307/1‑150 (MQ=255)
                                                      tGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGa                                                                                               >  2:3593569/1‑150 (MQ=255)
                                                             gAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTg                                                                                        >  1:1912090/1‑150 (MQ=255)
                                                             gAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTg                                                                                        >  2:509577/1‑150 (MQ=255)
                                                              aacaaaacaaCCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATCTCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGa                                                                                       >  1:2946678/1‑150 (MQ=255)
                                                                cAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACt                                                                                     <  2:1702223/150‑1 (MQ=255)
                                                                cAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACt                                                                                     <  2:1156173/150‑1 (MQ=255)
                                                                cAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACt                                                                                     <  2:1155931/150‑1 (MQ=255)
                                                                   aaCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGcc                                                                                  <  2:3313205/150‑1 (MQ=255)
                                                                   aaCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGcc                                                                                  >  2:2823244/1‑150 (MQ=255)
                                                                   aaCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGcc                                                                                  <  2:2600097/150‑1 (MQ=255)
                                                                   aaCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGCTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGcc                                                                                  >  2:2823272/1‑150 (MQ=255)
                                                                       aCCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTg                                                                              >  1:942611/1‑150 (MQ=255)
                                                                         ccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGa                                                                            >  2:2213769/1‑150 (MQ=255)
                                                                         ccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGa                                                                            >  2:2997052/1‑150 (MQ=255)
                                                                          cGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGAt                                                                           >  1:3025146/1‑150 (MQ=255)
                                                                          cGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGAt                                                                           >  1:740734/1‑150 (MQ=255)
                                                                             aTCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTAc                                                                        >  1:805382/1‑150 (MQ=255)
                                                                               cgcgGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGt                                                                      >  2:1560402/1‑150 (MQ=255)
                                                                                   ggCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGt                                                                  <  1:2955251/150‑1 (MQ=255)
                                                                                        gAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTa                                                             <  1:3236096/150‑1 (MQ=255)
                                                                                        gAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTa                                                             <  2:2041605/150‑1 (MQ=255)
                                                                                         aGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGTTATCTGGCGATCCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTAc                                                            <  1:54198/150‑1 (MQ=255)
                                                                                         aGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTAc                                                            >  2:1510243/1‑150 (MQ=255)
                                                                                          gCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACa                                                           <  2:3633952/150‑1 (MQ=255)
                                                                                            cGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAgg                                                         <  2:2563691/150‑1 (MQ=255)
                                                                                               tggctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATc                                                      <  1:3791627/150‑1 (MQ=255)
                                                                                                  ctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGa                                                   >  1:4351749/1‑150 (MQ=255)
                                                                                                   tggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGAc                                                  <  2:3825160/150‑1 (MQ=255)
                                                                                                      cGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCga                                               <  1:489558/150‑1 (MQ=255)
                                                                                                        ttCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCgaga                                             <  2:3220674/150‑1 (MQ=255)
                                                                                                        ttCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCgaga                                             <  2:1780232/150‑1 (MQ=255)
                                                                                                          cATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGc                                           <  2:2452278/150‑1 (MQ=255)
                                                                                                          cATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGACACACCAGGCGATTATCGATTTGGTGACTGCCACTGGATTACGACCGTGGTTACAGGATCTGACCGAGAGc                                           >  2:1658751/1‑150 (MQ=255)
                                                                                                            tGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTTGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGa                                         <  1:3344330/150‑1 (MQ=255)
                                                                                                              cTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                       >  1:1178306/1‑150 (MQ=255)
                                                                                                              cTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                       <  1:997319/150‑1 (MQ=255)
                                                                                                              cTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                       >  2:3177678/1‑150 (MQ=255)
                                                                                                               ttACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAca                                      <  1:2668165/150‑1 (MQ=255)
                                                                                                                tactacGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcag                                     <  1:2262819/150‑1 (MQ=255)
                                                                                                                tactacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcag                                     >  2:1114388/1‑150 (MQ=255)
                                                                                                                   tacGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcagcag                                  <  1:3276431/150‑1 (MQ=255)
                                                                                                                   tacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcagcag                                  <  2:3542378/150‑1 (MQ=255)
                                                                                                                   tacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcagcag                                  <  2:1470385/150‑1 (MQ=255)
                                                                                                                       atatTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCttt                              <  2:3219130/150‑1 (MQ=255)
                                                                                                                        tatTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCtttt                             <  1:2103673/150‑1 (MQ=255)
                                                                                                                        tatTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCtttt                             <  2:1279399/150‑1 (MQ=255)
                                                                                                                          ttttGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTc                           >  1:1613666/1‑150 (MQ=255)
                                                                                                                              gAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTaa                       >  1:1119936/1‑150 (MQ=255)
                                                                                                                              gAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTaa                       <  1:2945572/150‑1 (MQ=255)
                                                                                                                                gCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgc                     <  2:1868988/150‑1 (MQ=255)
                                                                                                                                gCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgc                     <  1:590262/150‑1 (MQ=255)
                                                                                                                                  gAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCga                                               <  2:1550137/122‑1 (MQ=255)
                                                                                                                                  gAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgcgc                   >  1:1361202/1‑150 (MQ=255)
                                                                                                                                     gCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTAc                <  2:3029382/150‑1 (MQ=255)
                                                                                                                                      ccGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTAcc               <  2:3796460/150‑1 (MQ=255)
                                                                                                                                        ggATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCAt             >  2:514531/1‑150 (MQ=255)
                                                                                                                                         gATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATc            <  2:3259529/150‑1 (MQ=255)
                                                                                                                                                   cGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTGGa  <  1:3903731/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTGGA  >  CP009273/1601983‑1602279

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: