Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,602,131 G→T 24.8% 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,1310GT24.8% 727.7 / 201.2 302D177Y (GAT→TAT) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (107/120);  new base T (30/45);  total (137/165)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

CAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTGGAA  >  CP009273/1601982‑1602280
                                                                                                                                                     |                                                                                                                                                     
  gTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGat                                                                                                                                                     <  1:3225359/150‑1 (MQ=255)
   tACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCtata                                                                                                                                                    >  2:4191736/1‑150 (MQ=255)
   tACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCtata                                                                                                                                                    >  2:4191959/1‑150 (MQ=255)
   tACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCtata                                                                                                                                                    >  2:4617545/1‑150 (MQ=255)
      aGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCt                                                                                                                                                 >  2:4891020/1‑150 (MQ=255)
      aGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCt                                                                                                                                                 >  2:2973295/1‑150 (MQ=255)
      aGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCt                                                                                                                                                 <  1:1118875/150‑1 (MQ=255)
         tGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGc                                                                                                                                              <  1:1526368/150‑1 (MQ=255)
             aGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAac                                                                                                                                          >  1:2623908/1‑150 (MQ=255)
                         ggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCa                                                                                                                              <  1:4032557/150‑1 (MQ=255)
                          gAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCAc                                                                                                                             <  1:1561251/150‑1 (MQ=255)
                          gAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCAc                                                                                                                             >  1:4859650/1‑150 (MQ=255)
                          gAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCAc                                                                                                                             <  2:673309/150‑1 (MQ=255)
                            gCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTATGAGCGAAGCCGGAGGTGAGGTCGATATCTGGCGAACCAACTACTATCACCa                                                                                                                           >  2:3952390/1‑150 (MQ=255)
                               gACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGAt                                                                                                                        <  1:3636169/150‑1 (MQ=255)
                               gACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGAt                                                                                                                        >  1:1901698/1‑150 (MQ=255)
                               gACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTAATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGAt                                                                                                                        <  1:5113773/150‑1 (MQ=255)
                                  ccATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGcc                                                                                                                     <  1:4301947/150‑1 (MQ=255)
                                      gtgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTca                                                                                                                 <  1:4123719/150‑1 (MQ=255)
                                      gtgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTca                                                                                                                 <  1:3659501/150‑1 (MQ=255)
                                       tgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTcac                                                                                                                >  2:1628616/1‑150 (MQ=255)
                                           cATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCa                                                                                                            <  2:2384782/150‑1 (MQ=255)
                                            aTGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAg                                                                                                           >  2:4973966/1‑150 (MQ=255)
                                            aTGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAg                                                                                                           <  1:3836995/150‑1 (MQ=255)
                                             tGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAgg                                                                                                          <  2:2452606/150‑1 (MQ=255)
                                             tGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAgg                                                                                                          <  1:2988550/150‑1 (MQ=255)
                                                gcgAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGa                                                                                                       <  1:3438940/150‑1 (MQ=255)
                                                    aGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTAt                                                                                                   <  1:4822366/150‑1 (MQ=255)
                                                     gTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATc                                                                                                  <  2:643910/150‑1 (MQ=255)
                                                      ttGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCg                                                                                                 <  2:3014469/150‑1 (MQ=255)
                                                      ttGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCg                                                                                                 >  2:1711498/1‑150 (MQ=255)
                                                      ttGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCg                                                                                                 >  1:604097/1‑150 (MQ=255)
                                                        gCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGAt                                                                                               <  2:4874504/150‑1 (MQ=255)
                                                         ccTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGAtt                                                                                              <  2:2702905/150‑1 (MQ=255)
                                                            gggAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTggg                                                                                           <  2:4706253/150‑1 (MQ=255)
                                                            gggAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTggg                                                                                           <  1:2373836/150‑1 (MQ=255)
                                                            gggAACAAAACGCCCCGGAGCGCGGGCGTGAGCCGTTGGCTGGCGGTCATGCTTACTACGAGATTTTGCGCGAAGCCGGATGTGAGGTCGATATCCGGCGCACCACCTACTCTCACCAGATGCCGTCACACCAGGCGATTATCGATTggg                                                                                           <  2:3669303/150‑1 (MQ=255)
                                                             ggAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGt                                                                                          <  2:5045164/150‑1 (MQ=255)
                                                              gAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTg                                                                                         <  2:3455925/150‑1 (MQ=255)
                                                                 cAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACt                                                                                      >  2:1930760/1‑150 (MQ=255)
                                                                     aCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCa                                                                                  >  1:783207/1‑150 (MQ=255)
                                                                       tACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACt                                                                                <  2:714446/150‑1 (MQ=255)
                                                                           cGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGAt                                                                            >  1:2426954/1‑150 (MQ=255)
                                                                            ggATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGAGATTTTGAGCGACGCCGGATGGGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGagt                                                                           >  2:1578693/1‑148 (MQ=255)
                                                                              aTCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTAc                                                                         >  1:2872934/1‑150 (MQ=255)
                                                                              aTCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTAc                                                                         >  1:2819184/1‑150 (MQ=255)
                                                                              aTCGCGGGCGTGAGCCGTTGGCTGGAGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTAc                                                                         >  2:1328245/1‑150 (MQ=255)
                                                                                cgcgGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGt                                                                       >  1:1318230/1‑150 (MQ=255)
                                                                                   gggCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCg                                                                    >  1:2910985/1‑150 (MQ=255)
                                                                                      cGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTgg                                                                 <  1:696101/150‑1 (MQ=255)
                                                                                       gTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGt                                                                <  1:1233217/150‑1 (MQ=255)
                                                                                       gTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGt                                                                >  2:1545980/1‑150 (MQ=255)
                                                                                        tGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGtt                                                               >  2:2056109/1‑150 (MQ=255)
                                                                                        tGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGtt                                                               >  1:1636012/1‑150 (MQ=255)
                                                                                        tGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGtt                                                               >  1:3181062/1‑150 (MQ=255)
                                                                                         gAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTa                                                              >  1:3378938/1‑150 (MQ=255)
                                                                                         gAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTa                                                              >  2:3079891/1‑150 (MQ=255)
                                                                                           gCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACa                                                            >  1:4137915/1‑150 (MQ=255)
                                                                                            ccGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAg                                                           <  1:447864/150‑1 (MQ=255)
                                                                                            ccGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAg                                                           <  1:2774421/150‑1 (MQ=255)
                                                                                              gTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGa                                                         <  1:1685452/150‑1 (MQ=255)
                                                                                                tggctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATc                                                       >  2:1490140/1‑150 (MQ=255)
                                                                                                tggctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATc                                                       >  2:1895699/1‑150 (MQ=255)
                                                                                                   ctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGAGGGGAGGTCGATATCTGGCGAACCACCTACTATCACCCGATGCCGTCACACCAGGCGATTCTCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGa                                                    >  2:2645631/1‑150 (MQ=255)
                                                                                                    tggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGAc                                                   >  2:1491675/1‑150 (MQ=255)
                                                                                                     ggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGAcc                                                  >  2:4504564/1‑150 (MQ=255)
                                                                                                     ggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGAcc                                                  <  2:3544066/150‑1 (MQ=255)
                                                                                                     ggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGAcc                                                  <  1:2687335/150‑1 (MQ=255)
                                                                                                       cGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCga                                                >  1:3241845/1‑150 (MQ=255)
                                                                                                       cGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCga                                                <  2:1171565/150‑1 (MQ=255)
                                                                                                        gTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCgag                                               <  2:2080649/150‑1 (MQ=255)
                                                                                                         ttCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCgaga                                              >  1:3818985/1‑150 (MQ=255)
                                                                                                         ttCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCgaga                                              <  2:4310937/150‑1 (MQ=255)
                                                                                                          tCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGagag                                             <  2:2736201/150‑1 (MQ=255)
                                                                                                          tCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGagag                                             >  1:4152946/1‑150 (MQ=255)
                                                                                                          tCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGagag                                             <  1:4891020/150‑1 (MQ=255)
                                                                                                           cATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGc                                            <  2:3892993/150‑1 (MQ=255)
                                                                                                            aTGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCg                                           >  2:3265834/1‑150 (MQ=255)
                                                                                                               cTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                        >  1:1387242/1‑150 (MQ=255)
                                                                                                               cTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                        >  1:4646543/1‑150 (MQ=255)
                                                                                                               cTTACTACGATATTGTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                        <  2:3235462/150‑1 (MQ=255)
                                                                                                                ttACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAca                                       <  1:5046315/150‑1 (MQ=255)
                                                                                                                 tactacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcag                                      <  2:2075852/150‑1 (MQ=255)
                                                                                                                 tactacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcag                                      >  2:1078360/1‑150 (MQ=255)
                                                                                                                    tacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCCCTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcagcag                                   >  1:499632/1‑150 (MQ=255)
                                                                                                                     acGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACagcagc                                  <  2:3810019/150‑1 (MQ=255)
                                                                                                                          atTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCttttt                             >  2:3599498/1‑150 (MQ=255)
                                                                                                                                 gCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgc                      >  1:955553/1‑150 (MQ=255)
                                                                                                                                   gAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgcgc                    >  1:1082775/1‑150 (MQ=255)
                                                                                                                                   gAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgcgc                    >  2:4230331/1‑150 (MQ=255)
                                                                                                                                    aaGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCt                   <  1:2097260/150‑1 (MQ=255)
                                                                                                                                       ccGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTAcc                <  1:252906/150‑1 (MQ=255)
                                                                                                                                            tgtgAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAg           <  2:3130354/150‑1 (MQ=255)
                                                                                                                                              tgAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGAt         <  1:1856117/150‑1 (MQ=255)
                                                                                                                                                aGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGc       >  2:1785889/1‑150 (MQ=255)
                                                                                                                                                   tCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTgg    <  2:3378938/150‑1 (MQ=255)
                                                                                                                                                   tCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTgg    >  2:2555120/1‑150 (MQ=255)
                                                                                                                                                   tCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTgg    >  2:2302543/1‑150 (MQ=255)
                                                                                                                                                    cGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTGGa   >  2:104026/1‑150 (MQ=255)
                                                                                                                                                     gATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTGGaa  <  2:2031118/150‑1 (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: