Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,602,131 G→T 20.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,1310GT20.5% 682.0 / 128.1 249D177Y (GAT→TAT) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (107/91);  new base T (22/29);  total (129/120)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

CAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTGGA  >  CP009273/1601982‑1602279
                                                                                                                                                     |                                                                                                                                                    
 aGTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGa                                                                                                                                                     >  2:1345449/1‑150 (MQ=255)
    aCAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGatat                                                                                                                                                  >  2:3465517/1‑150 (MQ=255)
     cAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATc                                                                                                                                                 <  1:721047/150‑1 (MQ=255)
     cAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATc                                                                                                                                                 >  1:622162/1‑150 (MQ=255)
       gATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTg                                                                                                                                               <  2:1092587/150‑1 (MQ=255)
          gCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCg                                                                                                                                            >  2:1068675/1‑150 (MQ=255)
           ccAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGa                                                                                                                                           >  1:1058722/1‑150 (MQ=255)
             aGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAac                                                                                                                                         >  2:2735581/1‑150 (MQ=255)
             aGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCCTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAac                                                                                                                                         <  2:1452576/150‑1 (MQ=255)
               aTAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAacca                                                                                                                                       <  2:1858824/150‑1 (MQ=255)
               aTAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAacca                                                                                                                                       >  1:2846391/1‑150 (MQ=255)
                 aaCTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAaccacc                                                                                                                                     <  1:2285494/150‑1 (MQ=255)
                       ctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTAt                                                                                                                               <  1:873766/150‑1 (MQ=255)
                       ctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTAt                                                                                                                               >  2:47990/1‑150 (MQ=255)
                        tggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATc                                                                                                                              <  2:1658529/150‑1 (MQ=255)
                          gAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCAc                                                                                                                            <  1:3636776/150‑1 (MQ=255)
                          gAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCAc                                                                                                                            <  1:3943139/150‑1 (MQ=255)
                           aGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCAcc                                                                                                                           >  1:1727628/1‑150 (MQ=255)
                                  ccATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGcc                                                                                                                    <  2:3175086/150‑1 (MQ=255)
                                   cATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCg                                                                                                                   >  2:4278462/1‑150 (MQ=255)
                                     tgtgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTc                                                                                                                 >  2:1040300/1‑150 (MQ=255)
                                     tgtgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTc                                                                                                                 >  1:3741407/1‑150 (MQ=255)
                                      gtgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTca                                                                                                                <  1:4050013/150‑1 (MQ=255)
                                            aTGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAg                                                                                                          <  2:102074/150‑1 (MQ=255)
                                             tGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAgg                                                                                                         <  2:2955651/150‑1 (MQ=255)
                                                 cgAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGAt                                                                                                     <  2:3442541/150‑1 (MQ=255)
                                                  gAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGAtt                                                                                                    <  2:359628/150‑1 (MQ=255)
                                                  gAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGagt                                                                                                    >  1:97378/1‑148 (MQ=255)
                                                  gAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGAtt                                                                                                    >  2:3830133/1‑150 (MQ=255)
                                                  gAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGAtt                                                                                                    >  2:293699/1‑150 (MQ=255)
                                                     gTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATTCTTAATACGATATTTTGAGCGAAGCCGGATGTGATGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATc                                                                                                 >  1:1447317/1‑150 (MQ=255)
                                                      ttGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCg                                                                                                <  2:961805/150‑1 (MQ=255)
                                                              gAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCCGATGCCGTCACACCAGGCGAGTATCGATTGGGTg                                                                                        >  1:1569774/1‑150 (MQ=255)
                                                                 cAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACt                                                                                     >  1:4314801/1‑150 (MQ=255)
                                                                 cAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACt                                                                                     <  1:3989772/150‑1 (MQ=255)
                                                                  aaaaCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTg                                                                                    <  2:724993/150‑1 (MQ=255)
                                                                   aaaCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCGTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGc                                                                                   >  1:3536284/1‑150 (MQ=255)
                                                                     aCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCa                                                                                 <  1:2765300/150‑1 (MQ=255)
                                                                     aCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCa                                                                                 <  1:2189399/150‑1 (MQ=255)
                                                                       tACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACt                                                                               <  1:1753349/150‑1 (MQ=255)
                                                                        aCCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCAcag                                                                              >  1:1622295/1‑148 (MQ=255)
                                                                        aCCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTg                                                                              <  1:3020609/150‑1 (MQ=255)
                                                                        aCCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTg                                                                              <  1:615650/150‑1 (MQ=255)
                                                                         cccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTgg                                                                             >  2:3865589/1‑150 (MQ=255)
                                                                          ccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGa                                                                            <  2:3307513/150‑1 (MQ=255)
                                                                             gATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTa                                                                         <  1:184802/150‑1 (MQ=255)
                                                                              aTCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTAc                                                                        <  1:3873697/150‑1 (MQ=255)
                                                                                cgcgGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGt                                                                      <  1:2056145/150‑1 (MQ=255)
                                                                                cgcgGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGt                                                                      <  1:3545655/150‑1 (MQ=255)
                                                                                   gggCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCg                                                                   <  1:4002182/150‑1 (MQ=255)
                                                                                    ggCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGt                                                                  >  2:3858592/1‑150 (MQ=255)
                                                                                     gCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGACGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTg                                                                 <  2:3544779/150‑1 (MQ=255)
                                                                                      cGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTgg                                                                <  1:2245253/150‑1 (MQ=255)
                                                                                       gTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGt                                                               >  1:2340790/1‑150 (MQ=255)
                                                                                       gTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGt                                                               <  1:2897056/150‑1 (MQ=255)
                                                                                        tGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAATCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACAGGATTACGTCCGTGGtt                                                              >  1:731758/1‑150 (MQ=255)
                                                                                        tGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACTAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGtt                                                              <  2:3009654/150‑1 (MQ=255)
                                                                                        tGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGtt                                                              <  1:170310/150‑1 (MQ=255)
                                                                                         gAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTa                                                             >  2:4057884/1‑150 (MQ=255)
                                                                                           gCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACa                                                           <  2:1157136/150‑1 (MQ=255)
                                                                                              gTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGa                                                        <  1:1086693/150‑1 (MQ=255)
                                                                                              gTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGa                                                        <  2:2114795/150‑1 (MQ=255)
                                                                                                 ggctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCt                                                     <  1:1746088/150‑1 (MQ=255)
                                                                                                   ctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGa                                                   >  2:463386/1‑150 (MQ=255)
                                                                                                   ctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGa                                                   <  1:423272/150‑1 (MQ=255)
                                                                                                    tggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGAc                                                  <  1:2771238/150‑1 (MQ=255)
                                                                                                        gTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCgag                                              >  2:2745690/1‑150 (MQ=255)
                                                                                                          tCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGagag                                            >  2:394048/1‑150 (MQ=255)
                                                                                                            aTGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCg                                          >  1:2579660/1‑150 (MQ=255)
                                                                                                            aTGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCg                                          <  1:1104717/150‑1 (MQ=255)
                                                                                                             tGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGa                                         <  2:2508520/150‑1 (MQ=255)
                                                                                                              gCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGaa                                        <  2:4276664/150‑1 (MQ=255)
                                                                                                               cTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                       >  2:226821/1‑150 (MQ=255)
                                                                                                               cTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                       >  2:816747/1‑150 (MQ=255)
                                                                                                               cTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                       <  2:1456232/150‑1 (MQ=255)
                                                                                                                  actacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcagc                                    >  1:2323396/1‑150 (MQ=255)
                                                                                                                      cGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCt                                <  2:1940134/150‑1 (MQ=255)
                                                                                                                      cGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCt                                >  2:3222341/1‑150 (MQ=255)
                                                                                                                       gATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCtt                               <  1:1971424/150‑1 (MQ=255)
                                                                                                                         tatTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCtttt                             >  2:990239/1‑150 (MQ=255)
                                                                                                                           ttttGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTc                           >  1:3218375/1‑150 (MQ=255)
                                                                                                                            tttGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCt                          >  1:3724630/1‑150 (MQ=255)
                                                                                                                                aGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAg                      <  2:375347/150‑1 (MQ=255)
                                                                                                                                aGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAg                      <  1:67099/150‑1 (MQ=255)
                                                                                                                                 gCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgc                     <  1:1349020/150‑1 (MQ=255)
                                                                                                                                  cGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgcg                    >  1:1045912/1‑150 (MQ=255)
                                                                                                                                   gAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgcgc                   >  2:2940184/1‑150 (MQ=255)
                                                                                                                                     aGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTa                 <  2:2871702/150‑1 (MQ=255)
                                                                                                                                     aGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTa                 <  2:333421/150‑1 (MQ=255)
                                                                                                                                     aGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTa                 >  2:1749230/1‑150 (MQ=255)
                                                                                                                                     aGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTa                 >  2:4158656/1‑150 (MQ=255)
                                                                                                                                      gCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTAc                <  1:4082542/150‑1 (MQ=255)
                                                                                                                                       ccGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTAcc               <  2:2844513/150‑1 (MQ=255)
                                                                                                                                        cGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCa              >  1:3568760/1‑150 (MQ=255)
                                                                                                                                         ggATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCAt             >  1:1955657/1‑150 (MQ=255)
                                                                                                                                         ggATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCAt             >  1:862809/1‑150 (MQ=255)
                                                                                                                                          gATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATc            >  1:1509919/1‑150 (MQ=255)
                                                                                                                                          gATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATc            <  1:2576221/150‑1 (MQ=255)
                                                                                                                                                  gTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTg    >  2:1978478/1‑150 (MQ=255)
                                                                                                                                                   tCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTgg   >  2:920069/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTGGA  >  CP009273/1601982‑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: