Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,602,131 G→T 5.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,1310GT 5.5% 1032.4 / 27.7 271D177Y (GAT→TAT) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (120/135);  new base T (7/8);  total (128/143)
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.07e-01

Only 100 of 278 total aligned reads displayed.

CAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTGGAA  >  CP009273/1601982‑1602280
                                                                                                                                                     |                                                                                                                                                     
 aGTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGa                                                                                                                                                      >  1:3974896/1‑150 (MQ=255)
  gTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGat                                                                                                                                                     >  2:1282012/1‑150 (MQ=255)
    aCAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGatat                                                                                                                                                   <  2:2949245/150‑1 (MQ=255)
    aCAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGatat                                                                                                                                                   <  2:3728336/150‑1 (MQ=255)
     cAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATc                                                                                                                                                  >  2:1857135/1‑150 (MQ=255)
     cAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATc                                                                                                                                                  <  1:4259674/150‑1 (MQ=255)
      aGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCt                                                                                                                                                 >  1:1268384/1‑150 (MQ=255)
      aGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCt                                                                                                                                                 <  2:881147/150‑1 (MQ=255)
        aTGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTgg                                                                                                                                               >  1:957066/1‑150 (MQ=255)
          gCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCg                                                                                                                                             <  1:768732/150‑1 (MQ=255)
           ccAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGa                                                                                                                                            >  2:2887807/1‑150 (MQ=255)
             aGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAac                                                                                                                                          >  2:2147211/1‑150 (MQ=255)
              gATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAacc                                                                                                                                         <  1:3780724/150‑1 (MQ=255)
               aTAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAacca                                                                                                                                        >  2:684672/1‑150 (MQ=255)
               aTAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAacca                                                                                                                                        <  2:3688003/150‑1 (MQ=255)
                tAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAaccac                                                                                                                                       <  2:3350624/150‑1 (MQ=255)
                  aCTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACct                                                                                                                                     >  1:2561389/1‑150 (MQ=255)
                     ggctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACctact                                                                                                                                  <  1:2407466/150‑1 (MQ=255)
                      gctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACctacta                                                                                                                                 <  1:1672746/150‑1 (MQ=255)
                       ctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTAt                                                                                                                                >  2:4303000/1‑150 (MQ=255)
                         ggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCa                                                                                                                              >  2:1330845/1‑150 (MQ=255)
                             ccGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAg                                                                                                                          >  1:2884853/1‑150 (MQ=255)
                                aCCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATg                                                                                                                       <  1:4273548/150‑1 (MQ=255)
                                aCCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATg                                                                                                                       >  1:3640485/1‑150 (MQ=255)
                                  ccATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGcc                                                                                                                     >  1:100671/1‑150 (MQ=255)
                                   cATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCg                                                                                                                    <  2:156182/150‑1 (MQ=255)
                                   cATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCg                                                                                                                    <  2:676113/150‑1 (MQ=255)
                                    aTGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGt                                                                                                                   >  1:3557817/1‑150 (MQ=255)
                                            aTGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAg                                                                                                           >  1:4155167/1‑150 (MQ=255)
                                            aTGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAg                                                                                                           <  2:2011751/150‑1 (MQ=255)
                                            aTGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAg                                                                                                           <  1:2609918/150‑1 (MQ=255)
                                              gcgcGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGc                                                                                                         <  1:60698/150‑1 (MQ=255)
                                                 cgAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGAt                                                                                                      <  2:1085934/150‑1 (MQ=255)
                                                  gAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGAtt                                                                                                     >  2:2486245/1‑150 (MQ=255)
                                                    aGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTAt                                                                                                   <  2:4464201/150‑1 (MQ=255)
                                                     ggtGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATc                                                                                                  <  1:2729103/148‑1 (MQ=255)
                                                      ttGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCg                                                                                                 <  2:2404260/150‑1 (MQ=255)
                                                      ttGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCg                                                                                                 <  1:389178/150‑1 (MQ=255)
                                                      ttGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCg                                                                                                 <  2:4167957/150‑1 (MQ=255)
                                                        gCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGAt                                                                                               <  1:1867434/150‑1 (MQ=255)
                                                          cTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTg                                                                                             <  1:1749407/150‑1 (MQ=255)
                                                            gggAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTggg                                                                                           <  1:1200860/150‑1 (MQ=255)
                                                            gggAACAAAACTACCCGGATCGCGGGCATGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTggg                                                                                           <  1:4468988/150‑1 (MQ=255)
                                                              gAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTg                                                                                         <  1:1040185/150‑1 (MQ=255)
                                                              gAACAAAACTACCCGGATCGCGGGCGTGAGCCATTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTg                                                                                         <  2:1074223/150‑1 (MQ=255)
                                                                 cAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACt                                                                                      <  2:46156/150‑1 (MQ=255)
                                                                 cAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACt                                                                                      <  1:4028665/150‑1 (MQ=255)
                                                                  aaaaCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTCGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTg                                                                                     <  2:2791923/150‑1 (MQ=255)
                                                                     aCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCa                                                                                  >  2:3831819/1‑150 (MQ=255)
                                                                      cTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCAc                                                                                 >  1:3937931/1‑150 (MQ=255)
                                                                        aCCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATAACTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCAcag                                                                               >  1:1822341/1‑148 (MQ=255)
                                                                         cccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCCCTgg                                                                              >  2:3210028/1‑150 (MQ=255)
                                                                         cccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTgg                                                                              >  2:2807344/1‑150 (MQ=255)
                                                                          ccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGa                                                                             >  2:305140/1‑150 (MQ=255)
                                                                          ccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGa                                                                             <  1:2128121/150‑1 (MQ=255)
                                                                          ccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATATGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGCTTGGGTGACTGCCACAGGa                                                                             >  1:1328393/1‑150 (MQ=255)
                                                                           cGGATCTCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATAGATTGGGTGACTGCCACTGGAt                                                                            >  1:2359050/1‑150 (MQ=255)
                                                                            ggATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGAtt                                                                           <  1:301995/150‑1 (MQ=255)
                                                                             gATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTa                                                                          <  1:3670907/150‑1 (MQ=255)
                                                                                        tGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGtt                                                               <  2:3584480/150‑1 (MQ=255)
                                                                                        tGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGtt                                                               <  2:3160888/150‑1 (MQ=255)
                                                                                        tGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAAACACCTACTATCACCAGATGCCGTCACACCAGGCGATTATAGATTTGGTGACTGCCACAGGATTAAGTACGTGGtt                                                               >  1:836145/1‑150 (MQ=255)
                                                                                        tGAGCCGTTGGCTGGCGTTCATGCTTACTAAGATATTTTGAGCGAATCCGGATGTGAGGTCGATATCTTGCGAACCACCTACTATCACCACATGCCGTCACAACAGGCGATTATAGATTTGGTGACTGCCACAGGATTAAGTCCCTGGtt                                                               >  1:2657651/1‑150 (MQ=255)
                                                                                           gCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACa                                                            <  1:807901/150‑1 (MQ=255)
                                                                                              gTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGa                                                         >  1:4400160/1‑150 (MQ=255)
                                                                                                tggctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATc                                                       <  2:2087833/150‑1 (MQ=255)
                                                                                                 ggctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCt                                                      <  2:4408866/150‑1 (MQ=255)
                                                                                                   ctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGa                                                    <  2:3075583/150‑1 (MQ=255)
                                                                                                    tggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGAc                                                   >  1:3289253/1‑150 (MQ=255)
                                                                                                       cGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCga                                                >  1:646387/1‑150 (MQ=255)
                                                                                                       cGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCga                                                >  2:612/1‑150 (MQ=255)
                                                                                                          tCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGagag                                             <  1:882909/150‑1 (MQ=255)
                                                                                                            aTGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCg                                           >  2:713428/1‑150 (MQ=255)
                                                                                                            aTGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCg                                           <  1:2287436/150‑1 (MQ=255)
                                                                                                              gCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGaa                                         <  1:1782247/150‑1 (MQ=255)
                                                                                                               cTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                        <  2:224226/150‑1 (MQ=255)
                                                                                                               cTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                        <  2:3414268/150‑1 (MQ=255)
                                                                                                                ttACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAca                                       <  1:220051/150‑1 (MQ=255)
                                                                                                                 tactacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcag                                      <  1:1486023/150‑1 (MQ=255)
                                                                                                                      cGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCt                                 <  2:257473/150‑1 (MQ=255)
                                                                                                                       gATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCtt                                >  2:1129672/1‑150 (MQ=255)
                                                                                                                       gATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCtt                                <  2:3077105/150‑1 (MQ=255)
                                                                                                                       gATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCtt                                >  2:318243/1‑150 (MQ=255)
                                                                                                                       gATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCtt                                <  2:150071/150‑1 (MQ=255)
                                                                                                                        atatTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCttt                               <  1:4501747/150‑1 (MQ=255)
                                                                                                                        atatTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCttt                               <  1:2807344/150‑1 (MQ=255)
                                                                                                                        atatTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAAGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCttt                               <  2:1712605/150‑1 (MQ=255)
                                                                                                                          atTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCttttt                             <  1:722914/150‑1 (MQ=255)
                                                                                                                          atTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCttttt                             <  1:3709216/150‑1 (MQ=255)
                                                                                                                           ttttGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTc                            <  2:1982853/150‑1 (MQ=255)
                                                                                                                                aGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAg                       <  1:3136510/150‑1 (MQ=255)
                                                                                                                                  cGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgcg                     >  1:3817932/1‑150 (MQ=255)
                                                                                                                                  cGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgcg                     >  1:203908/1‑150 (MQ=255)
                                                                                                                                    aaGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCt                   <  2:2701628/150‑1 (MQ=255)
                                                                                                                                      gCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTAc                 >  2:299783/1‑150 (MQ=255)
                                                                                                                                        cGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCa               >  1:3494382/1‑150 (MQ=255)
                                                                                                                                        cGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCa               >  2:919014/1‑150 (MQ=255)
                                                                                                                                           aTGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCa            >  1:3726517/1‑150 (MQ=255)
                                                                                                                                              tgAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGAt         <  2:2618800/150‑1 (MQ=255)
                                                                                                                                               gAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATg        >  1:1938992/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTGGAA  >  CP009273/1601982‑1602280

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: