Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 1,782,894 C→A A68E (GCA→GAA)  aroH → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,782,8940CA90.5% 210.0 / 15.4 74A68E (GCA→GAA) aroH3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base A (36/31);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGAT  >  CP009273/1782766‑1783024
                                                                                                                                |                                                                                                                                  
tCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGaaa                                                                                                                                  <  1:308740/131‑1 (MQ=255)
 cccGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTAc                                                                                                                          >  1:3746090/1‑138 (MQ=255)
     tAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAAAACTGAATGGGGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCa                                                                                                                      >  1:838031/1‑138 (MQ=255)
       aCGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCAccc                                                                                                                   >  1:2790349/1‑139 (MQ=255)
       aCGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCAccc                                                                                                                   >  1:2838901/1‑139 (MQ=255)
        cGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCg                                                                                                                  >  1:3752072/1‑139 (MQ=255)
          ccTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGc                                                                                                                        >  1:2481616/1‑131 (MQ=255)
                      ccccATGTCCCCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCACCCGTCTGCAGTCGCTgc                                                                                                    <  1:2098114/138‑1 (MQ=255)
                      ccccATACCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTgc                                                                                                    <  1:3646643/138‑1 (MQ=255)
                              ccccGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAaccag                                                                                            <  1:3421522/137‑1 (MQ=255)
                                cctacccccGCCGCAGAATTGAAAAAATACTGAATGGCGAAGATAAGCGCCTGTTTGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAaccagta                                                                                          <  1:2356609/133‑1 (MQ=255)
                                  gACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAaccagtac                                                                                         >  1:1084038/1‑138 (MQ=255)
                                   aCTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAaccagtacca                                                                                       >  1:1584672/1‑139 (MQ=255)
                                     tCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGGCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAaccagtaccag                                                                                      >  1:3525339/1‑138 (MQ=255)
                                      cccGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTc                                                                                    >  1:391595/1‑139 (MQ=255)
                                                aaTTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCCTTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTg                                                                                                     <  1:4137689/112‑1 (MQ=255)
                                                 aTTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGaaa                                                                         <  1:488820/139‑1 (MQ=255)
                                                 aTTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGaa                                                                          <  1:2564318/138‑1 (MQ=255)
                                                     aaaaaaTACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGaaa                                                                                                                                  <  1:3672190/78‑1 (MQ=255)
                                                           tACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCg                                                                                                                  <  1:3044563/88‑1 (MQ=255)
                                                               gAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGc                                                                                                                        >  1:1501318/1‑78 (MQ=255)
                                                               gAATGGTGAAGATAAGCGACTGTGGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCa                                                             >  1:1093375/1‑137 (MQ=255)
                                                                 aTGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTa                                                         <  1:1646676/139‑1 (MQ=255)
                                                                    gTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAac                                                                                               <  1:2505215/98‑1 (MQ=255)
                                                                    gTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATg                                                                 <  1:351730/128‑1 (MQ=255)
                                                                     tGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTc                                                                                                         <  1:1759680/87‑1 (MQ=255)
                                                                     tGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCg                                                                                                        >  1:872529/1‑88 (MQ=255)
                                                                       aaGATAAGCGCCTGTTGGTCAGTATTGGCCCTTGCTCGATCCTCGATCTCAACGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTgcg                                                                                                   <  1:4342631/91‑1 (MQ=255)
                                                                            aaGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGaaaaa                                               >  1:1698663/1‑138 (MQ=255)
                                                                             aGCGACTGTTGGTCATTATTGGCCCCTACTCGATCCTCGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATg                                                                 <  1:1545768/119‑1 (MQ=255)
                                                                              gCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAaccagtaccagt                                                                                     >  1:332384/1‑98 (MQ=255)
                                                                              gCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTAtt                                                       >  1:2524522/1‑128 (MQ=255)
                                                                               cGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAac                                                                                               >  1:3930007/1‑87 (MQ=255)
                                                                               cGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCAc                                           >  1:102288/1‑139 (MQ=255)
                                                                                gaatGTTGGTCAATATTGGCCCCTGATCGAGCCACGATCGCCCCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGCCCCGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCa                                            <  1:3010096/134‑1 (MQ=255)
                                                                                gACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCAc                                           <  1:3513448/138‑1 (MQ=255)
                                                                                gACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACg                                          <  1:3921462/139‑1 (MQ=255)
                                                                                 aCTGTTGGGCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAaccagtac                                                                                         >  1:2173145/1‑91 (MQ=255)
                                                                                   tGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGa                                         >  1:1131234/1‑137 (MQ=255)
                                                                                   tGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGa                                         >  1:31445/1‑137 (MQ=255)
                                                                                   tGTTGGTCATTAATGGCCCCTGCTCGATCCACGATCTCACCGCCGCAATGGAGTACGCCACCCGTCTGCAGTc                                                                                                         <  1:3774902/73‑1 (MQ=255)
                                                                                     ttGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAacc                                                                                              <  1:846021/82‑1 (MQ=255)
                                                                                       ggTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACtgtt                                   <  1:2296872/139‑1 (MQ=255)
                                                                                       ggTCATTATTGGCCCCTACTCGATCCACGATCTCACCGATGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCTGTACCAGTCACGGCTGGAAAGCGTAATGCGCACCTAtt                                                       >  1:3716546/1‑119 (MQ=255)
                                                                                            agacaggCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCgg                              <  1:2570310/134‑1 (MQ=255)
                                                                                             tattGGCCCCTGCTCGATCCACGATCTCACCGCTGaaa                                                                                                                                  <  1:3187315/38‑1 (MQ=255)
                                                                                             tattGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAac                                                                                               >  1:1604205/1‑73 (MQ=255)
                                                                                             tattGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGt                                                                                                          <  1:1098869/62‑1 (MQ=255)
                                                                                             tattGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAac                                                                                               <  1:210096/73‑1 (MQ=255)
                                                                                               ttGGCCCCTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTc                                                                                    >  1:3018059/1‑82 (MQ=255)
                                                                                                  ggcccTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGaa                        <  1:3988306/137‑1 (MQ=255)
                                                                                                   gcccTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGaaa                       >  1:1818068/2‑139 (MQ=255)
                                                                                                   ccccTGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGaa                        <  1:2530449/138‑1 (MQ=255)
                                                                                                       tGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGc                                                                                                                        >  1:1015502/1‑38 (MQ=255)
                                                                                                       tGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCaa                                                                                                >  1:3270550/1‑62 (MQ=255)
                                                                                                       tGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAaccagtaccagt                                                                                     >  1:2383114/1‑73 (MQ=255)
                                                                                                       tGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAgg                     >  1:3521019/1‑137 (MQ=255)
                                                                                                       tGCTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAgg                     >  1:2409876/1‑137 (MQ=255)
                                                                                                         cTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTa                                                                    <  1:2500327/88‑1 (MQ=255)
                                                                                                         cTCGATCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTa                 >  1:357666/1‑139 (MQ=255)
                                                                                                             aTCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGa                    >  1:397597/1‑132 (MQ=255)
                                                                                                              tCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAac                                                                                               <  1:2031878/56‑1 (MQ=255)
                                                                                                              tCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTCCCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGa                         <  1:1874051/126‑1 (MQ=255)
                                                                                                              tCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCg                                                                      <  1:512301/81‑1 (MQ=255)
                                                                                                              tCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATgcg                                                               <  1:3786025/88‑1 (MQ=255)
                                                                                                              tCCACGATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAAtctc            <  1:2892633/139‑1 (MQ=255)
                                                                                                                   gATCTCACCGCTGCAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCt                                                          >  1:327484/1‑88 (MQ=255)
                                                                                                                   gATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAac                                                                                               <  1:2189808/51‑1 (MQ=255)
                                                                                                                   gATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATc       <  1:3198115/139‑1 (MQ=255)
                                                                                                                   gATCTCACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATc       <  1:617940/139‑1 (MQ=255)
                                                                                                                        cACCGCTGAAATGGAGTACGCTACCCGTCCGCAGTCGCTGCGCAACCAGTACCAGTCTCGTCTGGAAATCGTAATGAGCACCTATTTTGAAAAACCAGGAACTGTTGTCGGCGGGAAATGACTAAt               >  1:4044177/1‑126 (MQ=255)
                                                                                                                        cACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAaccagtaccagt                                                                                     >  1:4201641/1‑56 (MQ=255)
                                                                                                                        cACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCAc                                                            >  1:2684886/1‑81 (MQ=255)
                                                                                                                        cACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTAtttt                                                     >  1:1615359/1‑88 (MQ=255)
                                                                                                                        cACCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGAt  >  1:1026317/1‑139 (MQ=255)
                                                                                                                         aCCGCTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCa                                                                                                 <  1:3928317/43‑1 (MQ=255)
                                                                                                                             cTGAAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAaccagtaccagt                                                                                     >  1:16460/1‑51 (MQ=255)
                                                                                                                                |                                                                                                                                  
TCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCACCCGTCTGCAGTCGCTGCGCAACCAGTACCAGTCACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGAT  >  CP009273/1782766‑1783024

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