Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,602,131 G→T 8.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,1310GT 8.8% 646.4 / 25.6 182D177Y (GAT→TAT) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (73/93);  new base T (7/9);  total (80/102)
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.96e-01

Only 100 of 183 total aligned reads displayed.

GTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTG  >  CP009273/1601984‑1602277
                                                                                                                                                   |                                                                                                                                                  
gTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCg                                                                                                                                                    <  1:2846268/148‑1 (MQ=255)
 tACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGa                                                                                                                                                   >  1:2149206/1‑148 (MQ=255)
  aCAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTat                                                                                                                                                  <  2:2278437/148‑1 (MQ=255)
    aGAAGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGatat                                                                                                                                                <  1:2424230/148‑1 (MQ=255)
     gATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATc                                                                                                                                               <  2:1261279/148‑1 (MQ=255)
      aTGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCt                                                                                                                                              >  2:8701/1‑148 (MQ=255)
      aTGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCt                                                                                                                                              >  1:427886/1‑148 (MQ=255)
      aTGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCt                                                                                                                                              <  1:955925/148‑1 (MQ=255)
        gCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTTAGCGAAGCCGGATGTGAGGTCGATATCTgg                                                                                                                                            >  1:2775920/1‑148 (MQ=255)
        gCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTgg                                                                                                                                            >  1:2775904/1‑148 (MQ=255)
            gATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGaa                                                                                                                                        >  1:1621421/1‑148 (MQ=255)
            gATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGaa                                                                                                                                        >  2:2722623/1‑148 (MQ=255)
                aCTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAaccac                                                                                                                                    >  2:1562912/1‑148 (MQ=255)
                aCTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAaccac                                                                                                                                    >  1:1263414/1‑148 (MQ=255)
                aCTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAaccac                                                                                                                                    <  2:760700/148‑1 (MQ=255)
                    gctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACctac                                                                                                                                <  2:2454334/148‑1 (MQ=255)
                    gctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACctac                                                                                                                                >  2:106548/1‑148 (MQ=255)
                    gctggAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACctac                                                                                                                                <  1:2437876/148‑1 (MQ=255)
                        gAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATc                                                                                                                            >  1:431273/1‑148 (MQ=255)
                             gACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAg                                                                                                                       <  1:6786/148‑1 (MQ=255)
                             gACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAg                                                                                                                       <  1:1285588/148‑1 (MQ=255)
                             gACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAg                                                                                                                       <  1:9666/148‑1 (MQ=255)
                             gACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAg                                                                                                                       <  2:1035558/148‑1 (MQ=255)
                              aCCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGa                                                                                                                      >  2:1744961/1‑148 (MQ=255)
                                ccATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATg                                                                                                                    <  1:1769615/148‑1 (MQ=255)
                                   tgtgCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCg                                                                                                                 <  2:631604/148‑1 (MQ=255)
                                      gCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTca                                                                                                              >  2:1999038/1‑148 (MQ=255)
                                         cATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCacac                                                                                                           >  2:783338/1‑148 (MQ=255)
                                         cATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCacac                                                                                                           >  2:783342/1‑148 (MQ=255)
                                         cATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCacac                                                                                                           >  2:2464023/1‑148 (MQ=255)
                                                   gTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTa                                                                                                 <  2:2056528/148‑1 (MQ=255)
                                                     tGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATc                                                                                               >  1:3313146/1‑148 (MQ=255)
                                                      gCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCg                                                                                              <  1:2447718/148‑1 (MQ=255)
                                                       ccTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGa                                                                                             <  2:1112231/148‑1 (MQ=255)
                                                        cTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGAt                                                                                            <  1:1518935/148‑1 (MQ=255)
                                                        cTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGAt                                                                                            <  1:1520168/148‑1 (MQ=255)
                                                          gggAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGAGTATCGATTg                                                                                          >  2:1229208/1‑148 (MQ=255)
                                                           ggAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGCTGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTgg                                                                                         >  2:3092098/1‑148 (MQ=255)
                                                            gAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTggg                                                                                        >  1:2979387/1‑148 (MQ=255)
                                                            gAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTggg                                                                                        <  1:537039/148‑1 (MQ=255)
                                                            gAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTggg                                                                                        >  2:1345955/1‑148 (MQ=255)
                                                              aCAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTg                                                                                      >  2:960910/1‑148 (MQ=255)
                                                                  aaCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTg                                                                                  <  1:1696407/148‑1 (MQ=255)
                                                                   aCTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGc                                                                                 <  1:1080636/148‑1 (MQ=255)
                                                                     tACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCa                                                                               >  1:2540300/1‑148 (MQ=255)
                                                                       cccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACt                                                                             >  2:1798512/1‑148 (MQ=255)
                                                                       cccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATAATCGATTGGGTGACTGCCACt                                                                             >  2:1798530/1‑148 (MQ=255)
                                                                        ccGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTg                                                                            <  1:905700/148‑1 (MQ=255)
                                                                            aTCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGAtt                                                                        <  2:431273/148‑1 (MQ=255)
                                                                             tCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTa                                                                       <  2:904211/148‑1 (MQ=255)
                                                                                cgGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGt                                                                    <  2:877211/148‑1 (MQ=255)
                                                                                    cGTGATCCGTTGTCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGATGTCTTTATCTGGCGAACCACCTACTATCACCAGTTGCCGTGACACCAGTAGATTATAGTTTGTGTGACTGCCACTGGATTACGTCCGt                                                                <  1:3159129/148‑1 (MQ=255)
                                                                                      tGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTgg                                                              >  1:3179484/1‑148 (MQ=255)
                                                                                      tGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTgg                                                              >  1:2715259/1‑148 (MQ=255)
                                                                                       gAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGt                                                             <  1:2654793/148‑1 (MQ=255)
                                                                                       gAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGt                                                             <  1:2654659/148‑1 (MQ=255)
                                                                                         gCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTa                                                           <  1:628690/148‑1 (MQ=255)
                                                                                         gCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTACGTGGTTa                                                           >  1:1463876/1‑148 (MQ=255)
                                                                                           cGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACa                                                         <  2:2085752/148‑1 (MQ=255)
                                                                                           cGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACa                                                         <  2:2323642/148‑1 (MQ=255)
                                                                                           cGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACa                                                         <  2:2323619/148‑1 (MQ=255)
                                                                                             ttGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAgg                                                       >  2:221538/1‑148 (MQ=255)
                                                                                              tggctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGa                                                      <  2:1904072/148‑1 (MQ=255)
                                                                                               ggctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGAt                                                     <  2:3152947/148‑1 (MQ=255)
                                                                                               ggctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGAt                                                     <  2:144908/148‑1 (MQ=255)
                                                                                               ggctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGAt                                                     <  2:1385730/148‑1 (MQ=255)
                                                                                                gctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATc                                                    >  2:2173314/1‑148 (MQ=255)
                                                                                                 ctggcGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCt                                                   <  2:2242267/148‑1 (MQ=255)
                                                                                                     cGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGAcc                                               <  2:1224154/148‑1 (MQ=255)
                                                                                                     cGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGAcc                                               >  2:2557280/1‑148 (MQ=255)
                                                                                                         cATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCgaga                                           >  2:566410/1‑148 (MQ=255)
                                                                                                            gCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCg                                        >  2:2855322/1‑148 (MQ=255)
                                                                                                              ttACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGaa                                      <  1:1570244/148‑1 (MQ=255)
                                                                                                               tactacGATATTTTGAGCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                     >  1:1977650/1‑148 (MQ=255)
                                                                                                               tactacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                     >  1:605999/1‑148 (MQ=255)
                                                                                                               tactacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAc                                     <  1:2795045/148‑1 (MQ=255)
                                                                                                                actacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAca                                    <  2:161794/148‑1 (MQ=255)
                                                                                                                 ctacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcag                                   <  2:2161350/148‑1 (MQ=255)
                                                                                                                  tacGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcagc                                  <  1:2169396/148‑1 (MQ=255)
                                                                                                                   acGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAAcagca                                 <  2:427886/148‑1 (MQ=255)
                                                                                                                     gATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACagcagc                               <  1:855547/148‑1 (MQ=255)
                                                                                                                     gATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACagcagc                               <  1:930003/148‑1 (MQ=255)
                                                                                                                          tttGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCttttt                          >  1:1722757/1‑148 (MQ=255)
                                                                                                                           ttGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCGCCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTc                         >  1:3145781/1‑148 (MQ=255)
                                                                                                                           ttGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTc                         <  1:2089268/148‑1 (MQ=255)
                                                                                                                           ttGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTc                         <  2:401910/148‑1 (MQ=255)
                                                                                                                           ttGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTc                         <  2:2294159/148‑1 (MQ=255)
                                                                                                                               gCGAAGCCGGATGTGAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTaa                     <  1:2023209/148‑1 (MQ=255)
                                                                                                                                 gAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgc                   >  2:3173484/1‑148 (MQ=255)
                                                                                                                                   aGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAgcgc                 >  2:2726276/1‑148 (MQ=255)
                                                                                                                                       ggATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTAcc             <  1:1376378/148‑1 (MQ=255)
                                                                                                                                        gATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCa            >  1:1955130/1‑148 (MQ=255)
                                                                                                                                        gATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCa            <  1:1345663/148‑1 (MQ=255)
                                                                                                                                           gtgAGGTCTATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCa         >  2:1444349/1‑148 (MQ=255)
                                                                                                                                            tgAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAg        >  2:1152063/1‑148 (MQ=255)
                                                                                                                                            tgAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAg        <  1:960910/148‑1 (MQ=255)
                                                                                                                                             gAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGa       <  1:1748563/148‑1 (MQ=255)
                                                                                                                                             gAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGa       <  1:106548/148‑1 (MQ=255)
                                                                                                                                               ggTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATg     <  1:1890922/148‑1 (MQ=255)
                                                                                                                                                gTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGc    >  2:2702890/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTACAGATGCCAGATAACTGGCTGGAGCCGACCCATGTGCTCATGCGCGAAGTTGCCTGGGAACAAAACTACCCGGATCGCGGGCGTGAGCCGTTGGCTGGCGTTCATGCTTACTACGATATTTTGAGCGAAGCCGGATGTGAGGTCGATATCTGGCGAACCACCTACTATCACCAGATGCCGTCACACCAGGCGATTATCGATTGGGTGACTGCCACTGGATTACGTCCGTGGTTACAGGATCTGACCGAGAGCGAACAGCAGCTTTTTCTTAAGCGCTACCATCAGATGCTG  >  CP009273/1601984‑1602277

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: ‑