Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,665,772 G→T 100% E123* (GAA→TAA)  hcaD → phenylpropionate dioxygenase, ferredoxin reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,665,7720GT95.3% 353.3 / 8.2 130E123* (GAA→TAA) hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base T (55/67);  total (59/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 130 total aligned reads displayed.

TCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGA  >  CP009273/2665627‑2665916
                                                                                                                                                 |                                                                                                                                                 
tCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATaa                                                                                                                                                 <  2:2462073/148‑1 (MQ=255)
 cTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAAc                                                                                                                                                >  2:2714552/1‑148 (MQ=255)
    cATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCGGCCGTTGCGTGTTGCACTGGGATAAccca                                                                                                                                             >  1:3279500/1‑145 (MQ=255)
     aTTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgctg                                                                                                                                            <  2:3263431/148‑1 (MQ=255)
      ttCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGAGCTGGCCCCCTGCCGTTGCGTGATGCACTGGGATAACgctgc                                                                                                                                           >  1:207482/1‑148 (MQ=255)
       tCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgctgct                                                                                                                                         <  2:2710073/148‑1 (MQ=255)
            tgtAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTAcc                                                                                                                                    >  2:654717/1‑148 (MQ=255)
             gtAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAccc                                                                                                                                   >  1:2994441/1‑148 (MQ=255)
               aaCCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTa                                                                                                                                 <  2:2729034/148‑1 (MQ=255)
                 ccATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACg                                                                                                                               <  1:1813848/148‑1 (MQ=255)
                   aTCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCCCAGCCGGTGCTTGATGCACTGGGCTAACGCTGCTTTAACCTACGcc                                                                                                                             >  1:1843254/1‑148 (MQ=255)
                     cAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCCCGGCCCGTGCGTGATGCACTGGGGGAACGCTGCTGTTACCGACGccca                                                                                                                           >  1:2296058/1‑146 (MQ=255)
                       aaaCGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGc                                                                                                                         <  1:509153/148‑1 (MQ=255)
                           gCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGt                                                                                                                     <  2:743317/148‑1 (MQ=255)
                               ggTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATg                                                                                                                 <  2:2544957/148‑1 (MQ=255)
                               ggTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCGTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTGACCCTACGCCATGCCGGTGATg                                                                                                                 <  2:524406/148‑1 (MQ=255)
                               ggTCGCGACACACGAGAGTGAGTGTTAACCACCGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATg                                                                                                                 <  2:2187071/148‑1 (MQ=255)
                                   gcgACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATgccgc                                                                                                             <  2:304330/148‑1 (MQ=255)
                                   gcgACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATgccgc                                                                                                             <  1:1138368/148‑1 (MQ=255)
                                     gACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCa                                                                                                           <  2:1879430/148‑1 (MQ=255)
                                       cacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGa                                                                                                         <  1:2667508/148‑1 (MQ=255)
                                       cacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGa                                                                                                         <  1:2663854/148‑1 (MQ=255)
                                        acacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGAc                                                                                                        <  1:82330/148‑1 (MQ=255)
                                        acacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGAc                                                                                                        <  1:1451997/148‑1 (MQ=255)
                                           cGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGCGCGCGACCCCGGCCCGTGGTTGTTTGACAGGGGTTAAGCGGCGTTAACCAACGCCAAGCCGGTGATGCCGCCAgcaagc                                                                                                     >  1:317652/1‑143 (MQ=255)
                                            gagaGTTAGTGGTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCg                                                                                                    <  2:2809781/148‑1 (MQ=255)
                                             agagTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCga                                                                                                   <  2:2411037/148‑1 (MQ=255)
                                              gagTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCGCGCCCGCTGCCGTTGGTTGATTCCCTGGGTTAACGCGGGTTTTACCTAAGCCCAGGCGGGGATGGCCCCCGACTgggcg                                                                                                  >  1:1308320/1‑144 (MQ=255)
                                                 ttAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAg                                                                                               <  1:2945308/148‑1 (MQ=255)
                                                 ttAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAg                                                                                               <  1:2945262/148‑1 (MQ=255)
                                                 ttAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAg                                                                                               <  1:2716629/148‑1 (MQ=255)
                                                  tAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGt                                                                                              <  1:561333/148‑1 (MQ=255)
                                                  tAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCTGCCGGTGCTTTATTCCCTGGGGTAACGCGGCGTTAACCAAAGCCAAACCCGTGATTCCGCCAGAAAGCGAGAAGt                                                                                              >  1:1395048/1‑148 (MQ=255)
                                                   aGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGtt                                                                                             <  1:2527110/148‑1 (MQ=255)
                                                      gtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGCGCGCGACCGCGGCCGGTGGGTGTTTCACGGGGATCAAGCTGCGTTTCCCCACGCCCAGCCCGGGGTGCCGCCCGACAGCGCGAACGTCTg                                                                                          >  1:2935673/1‑148 (MQ=255)
                                                          aCCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCTCGACCGCTGCCGTTGCGTGATGCCCTGGGATAACGCTGCTTTACCCAACGCCAGGCCGGTGATGCCGCCAGACAGCGAGAAGTTCTGCAGc                                                                                      >  2:817311/1‑148 (MQ=255)
                                                          aCCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCTCGACCGCTGCCGTTGCGTGATGCACTGGGATAACGCTGCTTTACCCGACCCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGc                                                                                      >  2:1865742/1‑148 (MQ=255)
                                                            cAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCTGCCGTTGCTTGTTGCACTGGGATAACCCTGCTTTACCCTACGCCAAGCCGGGGATGCCGCCCGACTGCGAGAAGTTGTGCAGccc                                                                                    >  2:3222927/1‑148 (MQ=255)
                                                              ccGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGGTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGa                                                                                  <  1:121367/147‑1 (MQ=255)
                                                              aCGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGa                                                                                  <  1:1794910/148‑1 (MQ=255)
                                                                    aaaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGTCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTc                                                                            <  1:1646113/148‑1 (MQ=255)
                                                                    aaaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCCTTGCATGAAGCACAGGGATAAACCTGCTGTTCCCAAAGCCAAGCCGGTGATTCCGCCAGACAGCGACAAGTTCTGCAGCCCGAAccctc                                                                            >  2:1937826/1‑144 (MQ=255)
                                                                     aaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCa                                                                           <  2:2748377/148‑1 (MQ=255)
                                                                      aGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCCGTGCATGATGCACGGGGCTAACGCGGCTTTACCCAACGCCAAGCCGGTGATGCCCCCCGACAGCCAGAAGTTCAGCAGCCCCAACCGTccg                                                                          >  2:1544511/1‑146 (MQ=255)
                                                                       gCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGt                                                                         <  1:102886/148‑1 (MQ=255)
                                                                           gCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATTCACTGGGATAACGCTGCGTTTACCAACCCCAAGCCGGTGATTCCGCCAGACAGCCAGAAGTTATGCAGCCCCAAACGTCAGGCGTg                                                                     >  1:2178119/1‑148 (MQ=255)
                                                                           gCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGCTAACGCTGCGTTACCCTACGCCCAGCCGGTGATTCCGCCAGACAGCGAGAAGTTCGGCAGCCCCAACCGTCAGTCGTg                                                                     >  1:2178226/1‑148 (MQ=255)
                                                                            cACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGa                                                                    <  2:2433163/148‑1 (MQ=255)
                                                                            cACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGGTGCTTGATTCACTGGGATAACCCTGCTTTAACCTACGCCATGCCCGTGATGACGCCAGACAGCGAGAACTTCTTCAGCCCGAACGGTCAGTCGTGa                                                                    >  2:129797/1‑148 (MQ=255)
                                                                              cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGAtt                                                                  <  2:1892694/148‑1 (MQ=255)
                                                                              cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGAtt                                                                  <  2:3220011/148‑1 (MQ=255)
                                                                              cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGAtt                                                                  <  2:1579246/148‑1 (MQ=255)
                                                                                gggATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTAt                                                                <  1:1931711/148‑1 (MQ=255)
                                                                                  gATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCg                                                              <  1:2666425/148‑1 (MQ=255)
                                                                                   aTCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCgg                                                             <  1:3205487/148‑1 (MQ=255)
                                                                                   aTCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCgg                                                             <  1:3204640/148‑1 (MQ=255)
                                                                                     cAGCTTTTTATAGCAACCGGCGCGGGCGCTCGCCCCCGGCCCGTGGTTGTTTGCCGGGGAGAAACCTGCGTTTACCTACGCCCAGCCCGGGGTGCCGCCCGACAGCGCGAAGTTATGCAACCCCAACCGTCAGGCGGGATTAtatgtg                                                           >  2:564440/1‑143 (MQ=255)
                                                                                       gCTTTTTATAGCAACCGGCGCGGAAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGcc                                                         >  1:2645496/1‑148 (MQ=255)
                                                                                        cTTTTTATAGCAACCGGCGCGGCCGCTCGACCGCGGCCGTTGCTTGATGCACTGGGATAACGCTGCGTTACCCTACGCCAAGCCGGTGATGCCGCCAGACTGCGCGAAGTTCTGCAGCCCGAACGGTCAGTCGGGAGTATAGGTGCCg                                                        >  2:1847069/1‑148 (MQ=255)
                                                                                          ttttATAGCAACCGGCGCGGCCGCGCGCCCCCGGCCGGTGGTTGATGGAAATGGATAACGCAGCTTTTACCAAACGCCAGCCCGGGATGGCGCCCCACTGCGAGAACTTCAGCAGCCCGAACCCTCACTCGTGTTTTTAGGTGCCCGa                                                      >  2:3051882/1‑148 (MQ=255)
                                                                                          ttttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGa                                                      <  2:3276387/148‑1 (MQ=255)
                                                                                             tataGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACt                                                   <  1:654717/148‑1 (MQ=255)
                                                                                              ataGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTa                                                  <  2:1962421/148‑1 (MQ=255)
                                                                                              ataGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTa                                                  <  2:1123035/148‑1 (MQ=255)
                                                                                                aGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAtt                                                <  1:1892823/148‑1 (MQ=255)
                                                                                                   aaCCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGt                                             <  1:1951476/148‑1 (MQ=255)
                                                                                                   aaCCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGt                                             <  1:1846844/148‑1 (MQ=255)
                                                                                                     ccGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCAGCTTTACCCTACGCCAAGCCGGTGATTCCGCCCGACTTCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGCTTATCGGTGCCGGAACTATTGGTCt                                           >  2:2363516/1‑148 (MQ=255)
                                                                                                        gcgcGGCAGCTCGACCGCTGCCGTTGCTTGATTCACTGGGATAACGCTGCTTTAACCAACGCCAAGCCGGTGATGCCGCCCGACAGCGAGAAGTTCTGCACCCCCAACGGTCAGTCGTGATTATAGGTGCCGGAACAATTGGTATGGa                                        >  2:3180548/1‑148 (MQ=255)
                                                                                                         cgcgGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGaa                                       <  2:2729287/148‑1 (MQ=255)
                                                                                                            ggCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTg                                    >  2:3169954/1‑148 (MQ=255)
                                                                                                             gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGCACTATTGGTCTGGAACTgg                                   >  1:2223913/1‑148 (MQ=255)
                                                                                                                 cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                               <  1:3091308/148‑1 (MQ=255)
                                                                                                                 cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                               >  1:381933/1‑148 (MQ=255)
                                                                                                                 cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCAACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                               >  1:3073320/1‑148 (MQ=255)
                                                                                                                      ccGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgcg                          <  1:2579453/148‑1 (MQ=255)
                                                                                                                       cGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgcgc                         <  1:2714552/148‑1 (MQ=255)
                                                                                                                         cTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCa                       <  2:556546/148‑1 (MQ=255)
                                                                                                                         cTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCa                       <  2:556181/148‑1 (MQ=255)
                                                                                                                             cGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgca                   >  1:3123252/1‑148 (MQ=255)
                                                                                                                              gTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcag                  <  1:2186685/148‑1 (MQ=255)
                                                                                                                               ttGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGACCTGGCTGCCAGCGCCAcgcagc                 >  2:2025522/1‑148 (MQ=255)
                                                                                                                                tGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTTGTCTGGAACTGGCTGCCAGCGCCAcgcagcg                >  2:2742518/1‑148 (MQ=255)
                                                                                                                                tGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcg                >  2:2748260/1‑148 (MQ=255)
                                                                                                                                tGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcg                >  1:3071183/1‑148 (MQ=255)
                                                                                                                                  cTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgca              >  1:2196484/1‑148 (MQ=255)
                                                                                                                                  cTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgca              >  1:1490408/1‑148 (MQ=255)
                                                                                                                                   ttGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgcag             <  1:1818040/148‑1 (MQ=255)
                                                                                                                                   ttGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgcag             <  1:1816781/148‑1 (MQ=255)
                                                                                                                                   ttGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgcag             <  1:1304776/148‑1 (MQ=255)
                                                                                                                                    tGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa            <  1:2632659/148‑1 (MQ=255)
                                                                                                                                      aTGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATg          >  2:710396/1‑148 (MQ=255)
                                                                                                                                      aTGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATg          >  2:710423/1‑148 (MQ=255)
                                                                                                                                       tGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGt         <  2:3279573/148‑1 (MQ=255)
                                                                                                                                       tGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGt         <  2:3279500/148‑1 (MQ=255)
                                                                                                                                        gCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTa        >  1:2586222/1‑148 (MQ=255)
                                                                                                                                         cACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTaa       <  2:317652/148‑1 (MQ=255)
                                                                                                                                         cACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTaa       <  2:317687/148‑1 (MQ=255)
                                                                                                                                             gggATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTg   >  2:771664/1‑148 (MQ=255)
                                                                                                                                              ggATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTCTCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGa  <  1:1515707/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGA  >  CP009273/2665627‑2665916

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