Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,012,953 C→A 45.1% R247S (CGT→AGT)  ubiE → bifunctional 2‑octaprenyl‑6‑methoxy‑1,4‑benzoquinone methylase/ S‑adenosylmethionine:2‑DMK methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,012,9530CA45.1% 155.5 / 218.4 201R247S (CGT→AGT) ubiEbifunctional 2‑octaprenyl‑6‑methoxy‑1,4‑benzoquinone methylase/ S‑adenosylmethionine:2‑DMK methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (51/58);  new base A (41/47);  total (92/109)
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.97e-01

Only 100 of 203 total aligned reads displayed.

GTCGCGAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGCGC  >  CP009273/4012806‑4013097
                                                                                                                                                   |                                                                                                                                                
gTCGCGAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATc                                                                                                                                                  >  1:1399705/1‑148 (MQ=255)
       acgacgCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTta                                                                                                                                           >  2:2268287/1‑148 (MQ=255)
       acgacgCCGACAGCTACCGTTATCGGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACGACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTta                                                                                                                                           <  1:3078986/148‑1 (MQ=255)
        cgacgCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTtat                                                                                                                                          >  2:2816166/1‑148 (MQ=255)
          acgCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATaa                                                                                                                                        <  2:678933/148‑1 (MQ=255)
             ccGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGtt                                                                                                                                     <  2:2065499/148‑1 (MQ=255)
             ccGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGtt                                                                                                                                     <  1:1638178/148‑1 (MQ=255)
              cGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTc                                                                                                                                    <  1:2572589/148‑1 (MQ=255)
              cGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTc                                                                                                                                    <  1:2145223/148‑1 (MQ=255)
              cGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTc                                                                                                                                    <  1:885736/148‑1 (MQ=255)
                aCAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTg                                                                                                                                  >  1:422383/1‑148 (MQ=255)
                 cAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGa                                                                                                                                 <  1:1038983/148‑1 (MQ=255)
                 cAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGa                                                                                                                                 <  1:1340094/148‑1 (MQ=255)
                     tACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAgg                                                                                                                             <  2:2924495/148‑1 (MQ=255)
                     tACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAgg                                                                                                                             <  2:2896198/148‑1 (MQ=255)
                       ccGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGgag                                                                                                                           >  1:3082017/1‑148 (MQ=255)
                          ttATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGAcc                                                                                                                        >  2:2464507/1‑148 (MQ=255)
                              cTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGaa                                                                                                                    <  1:199555/148‑1 (MQ=255)
                               tGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGaaa                                                                                                                   <  1:1128172/148‑1 (MQ=255)
                                 gCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATg                                                                                                                 >  2:747688/1‑148 (MQ=255)
                                  cAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGc                                                                                                                >  1:1674557/1‑148 (MQ=255)
                                     aaTCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCtt                                                                                                             <  1:3551539/148‑1 (MQ=255)
                                         catccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTaa                                                                                                         <  2:1819308/148‑1 (MQ=255)
                                          atccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTaaa                                                                                                        >  2:402640/1‑148 (MQ=255)
                                               tATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCttt                                                                                                   <  1:3420536/148‑1 (MQ=255)
                                                aTGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGGAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAATTGCGTTTTAAACCTTTa                                                                                                  >  1:2239440/1‑148 (MQ=255)
                                                   cATCCCGATCATGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTg                                                                                               <  2:1890501/148‑1 (MQ=255)
                                                   cATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTg                                                                                               <  1:2288084/148‑1 (MQ=255)
                                                   cATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTg                                                                                               <  1:493844/148‑1 (MQ=255)
                                                     tCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGAc                                                                                             >  1:1783370/1‑148 (MQ=255)
                                                      cccGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACg                                                                                            <  2:3231019/148‑1 (MQ=255)
                                                      cccGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACg                                                                                            >  2:48392/1‑148 (MQ=255)
                                                      cccGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACg                                                                                            >  1:1468741/1‑148 (MQ=255)
                                                       ccGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACgg                                                                                           <  1:908695/148‑1 (MQ=255)
                                                           tCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAgg                                                                                       >  1:1386766/1‑148 (MQ=255)
                                                            cAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGa                                                                                      <  2:1929382/148‑1 (MQ=255)
                                                            cAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGa                                                                                      <  2:1929027/148‑1 (MQ=255)
                                                            cAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGa                                                                                      <  1:836600/148‑1 (MQ=255)
                                                            cAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGa                                                                                      <  2:1472452/148‑1 (MQ=255)
                                                            cAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGa                                                                                      <  2:2842174/148‑1 (MQ=255)
                                                            cAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGa                                                                                      <  2:2844291/148‑1 (MQ=255)
                                                            cAGGAGACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGa                                                                                      <  2:2844236/148‑1 (MQ=255)
                                                             aGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGaa                                                                                     <  1:515829/148‑1 (MQ=255)
                                                               gATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAAtt                                                                                   >  1:633607/1‑148 (MQ=255)
                                                               gATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCGTTAGTGACGGCAGGAAtt                                                                                   >  2:2264125/1‑148 (MQ=255)
                                                                aTACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGa                                                                                      >  1:2089601/1‑144 (MQ=255)
                                                                aTACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTg                                                                                  >  2:619530/1‑148 (MQ=255)
                                                                 tACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGa                                                                                 <  1:78458/148‑1 (MQ=255)
                                                                  aCCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGaa                                                                                >  2:920736/1‑148 (MQ=255)
                                                                      tGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTc                                                                            <  1:2968562/148‑1 (MQ=255)
                                                                        aaaGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTg                                                                          <  1:1457215/148‑1 (MQ=255)
                                                                         aaGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGc                                                                         >  2:1261988/1‑148 (MQ=255)
                                                                          aGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCt                                                                        <  1:1187285/148‑1 (MQ=255)
                                                                          aGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCt                                                                        >  1:3311523/1‑148 (MQ=255)
                                                                          aGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCt                                                                        >  1:2350155/1‑148 (MQ=255)
                                                                           gCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTc                                                                       >  2:1718834/1‑148 (MQ=255)
                                                                                gatgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAacac                                                                  <  1:3195336/148‑1 (MQ=255)
                                                                                 atgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACAcc                                                                 >  1:2869173/1‑148 (MQ=255)
                                                                                  tgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCt                                                                >  1:2213985/1‑148 (MQ=255)
                                                                                   gCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCtt                                                               <  2:1647285/148‑1 (MQ=255)
                                                                                      ggATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCt                                                            >  1:232781/1‑148 (MQ=255)
                                                                                       gATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTg                                                           <  1:591531/148‑1 (MQ=255)
                                                                                          gCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTAt                                                        >  2:2101795/1‑148 (MQ=255)
                                                                                            cGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCg                                                      <  2:1143649/148‑1 (MQ=255)
                                                                                            cGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCg                                                      <  2:1966533/148‑1 (MQ=255)
                                                                                              gATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCt                                                    >  1:170516/1‑148 (MQ=255)
                                                                                               aTTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTc                                                   >  1:476096/1‑148 (MQ=255)
                                                                                                ttCGAAAGTGTCGACTACTACAATATGACGGAAGGGGTTGTGGCGCTGCATATTGGTTATAAGATAAGACAGGATACAGGAAATTCCTTTTAAACCTTTAGTGAAGTCAGGAATTTAAAGTATGATCAAAAACTTCATGTATCGCtaa                                                  >  1:1223591/1‑146 (MQ=255)
                                                                                                 tCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCAc                                                 >  2:821073/1‑148 (MQ=255)
                                                                                                 tCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCAc                                                 >  1:2423795/1‑148 (MQ=255)
                                                                                                 tCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCAc                                                 <  2:842819/148‑1 (MQ=255)
                                                                                                               tactacAATCTGACGGCGGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACTCCTTCCGGTATCGCTCACCCGCGCTGCCAACg                                   >  2:1828258/1‑148 (MQ=255)
                                                                                                                actacAATCTGCCGGCAGGGGTTGTGGCGCTGCATAGTGGTTATACGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACgg                                  <  1:3540616/148‑1 (MQ=255)
                                                                                                                  tacAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGcc                                >  2:258969/1‑148 (MQ=255)
                                                                                                                  tacAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGcc                                >  2:2616149/1‑148 (MQ=255)
                                                                                                                    cAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGTAAATGCATTTTAAACCTTTAGTGACGGCAGTAATTTAATGTCTGCTCAACACCTTCCTTTATCGCTCACCCGCGCTGAAAACGGCCCg                              <  2:49460/148‑1 (MQ=255)
                                                                                                                     aaTCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGc                             >  1:182472/1‑148 (MQ=255)
                                                                                                                           aCGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTATGACAGGAGACCGGAAATGCCTTTTAAACCGTTAGTGACGGCAGGAATAGAAAGTATTCTCAACACCTTCCATTATCGCTCACCCGCGCTGAAAACGCCCCGCTCGCGt                       >  1:3453321/1‑148 (MQ=255)
                                                                                                                              gctggggTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTAAAACCTTTAGTGACGGCAGTAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctg                    <  1:3512846/145‑1 (MQ=255)
                                                                                                                              gCAGGGGTTGTGGCGCTGAATTGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCGTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctg                    <  2:319605/148‑1 (MQ=255)
                                                                                                                               cAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctgc                   >  2:279921/1‑148 (MQ=255)
                                                                                                                               cAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctgc                   >  1:590532/1‑148 (MQ=255)
                                                                                                                                 ggggGTGTGGCGCTGCATAGTGTTTATAAGTTCTTACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctgctg                 <  2:1919611/148‑1 (MQ=255)
                                                                                                                                  gggTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTgg                >  2:396118/1‑148 (MQ=255)
                                                                                                                                  gggTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTgg                >  1:2936788/1‑148 (MQ=255)
                                                                                                                                  gggTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTgg                >  1:1157184/1‑148 (MQ=255)
                                                                                                                                  gggTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTgg                >  2:3111598/1‑148 (MQ=255)
                                                                                                                                      tgtgGCGGTCCTTTGTTGTTATAATTTCGGACAGGAGACCGGAATTGCCTTTTAAACCTTTCTTCACGGAAGGAATTTAATGTCTGCTCAACCCCTTCCTGTATCGCTCACCCGCGGTGAAAACGGCCCGCTCGCGTCTGCTGGGTaa            <  1:3396381/148‑1 (MQ=255)
                                                                                                                                          gcgcTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGCCGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTa        <  1:2816166/148‑1 (MQ=255)
                                                                                                                                          gcgcTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTa        >  1:657790/1‑148 (MQ=255)
                                                                                                                                          gcgcTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTa        <  2:2350155/148‑1 (MQ=255)
                                                                                                                                          gcgcTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTa        <  2:558122/148‑1 (MQ=255)
                                                                                                                                           cgcTGCTTCGTGGTTATCTGTTGTCACAGGAGCCCGGATTTGACTTTTAAATCTTTCGTGTCGTCAGGATTTTAAAGTCTGCTCAACCCCTTCGTTTATCGGTCCCCGGCGCTGAACACGGCCCGCTCGCGTCTGCTGGGTAAAGTAt       <  1:2442878/148‑1 (MQ=255)
                                                                                                                                           cgcTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCGGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTAt       <  2:2273968/148‑1 (MQ=255)
                                                                                                                                           cgcTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTAt       <  1:2268287/148‑1 (MQ=255)
                                                                                                                                            gcTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTAAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACCCCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTAtt      <  1:747441/148‑1 (MQ=255)
                                                                                                                                             cTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTg     >  1:1071583/1‑148 (MQ=255)
                                                                                                                                               gCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTgcg   >  1:710680/1‑148 (MQ=255)
                                                                                                                                               gCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTgcg   >  1:3157266/1‑148 (MQ=255)
                                                                                                                                                cATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCATTAGTGACGGCAGGAAATGAAATTCTGATCAACACCTTCCTGTATCGCTCACCAGCGCTGAAAACTTCCCGCTCGCGTCTGCTGGGTAAAGTATTgcgc  >  2:1689719/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTCGCGAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGCGC  >  CP009273/4012806‑4013097

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: