Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,103,793 G→T 8.5% A67S (GCG→TCG)  phdB → pyruvate dehydrogenase (acetyl‑transferring) beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,103,7930GT 8.5% 192.4 / 6.2 70A67S (GCG→TCG) phdBpyruvate dehydrogenase (acetyl‑transferring) beta subunit
Reads supporting (aligned to +/- strand):  ref base G (38/26);  new base T (4/2);  total (42/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

GTTTTAATTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCGGTACTAAAACTG  >  CP000730/1103664‑1103925
                                                                                                                                  |                                                                                                                                    
tttttAAGTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGAAAGAGTATTCGATACACCTTTAGCTGAAGAATGTATTGGTGGTTTAGCGATGGGTCtt                                                                                                                           <  2:635395/140‑1 (MQ=255)
gTTTTAATTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCTATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCtt                                                                                                                           <  2:450453/141‑1 (MQ=255)
gTTTTAATTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTGGCTGAATCAGGTATGGGTGGTTTAGCGATGGGTCtg                                                                                                                           >  1:293997/1‑140 (MQ=255)
 ttttAATTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGGAGTACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTGTTGGTGGTTTAGCGGTGGGTCTTg                                                                                                                          >  2:158874/1‑141 (MQ=255)
    tAAGTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAg                                                                                                                       <  2:469551/141‑1 (MQ=255)
    ataTTTTTGGTGAAGACGTGGGTGCTAACGGGGGTGTTTACCGAGATCATGAAGGACTACAAAAAGCCTTCGGTGAAGATCGAGTACTCGATACACCTTTCGCTGAATCAGGTATTGGTGCTTTAGCGATGGGTCTTGCAg                                                                                                                       <  1:487395/139‑1 (MQ=255)
       tttttGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTg                                                                                                                    >  2:648029/1‑141 (MQ=255)
          ttGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAg                                                                                                                 >  2:89606/1‑141 (MQ=255)
          ttGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAg                                                                                                                 <  1:429192/141‑1 (MQ=255)
          ttGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGAGGGGTCTTGCAGtgtacg                                                                                                                 >  2:226739/1‑136 (MQ=255)
           tggttgaGACGTTGGTGTTAACTGCGGTGTTTTCCGTGTTACTGAAGGAGTGAGAAATGGGTTTGGTGATGATAGAGTATTCGATGCACCTTTAGTTAATTCAGGTATTGGTGGTTGAGCGGTGTGTGTTGCGGTTGGAgg                                                                                                                >  2:534526/7‑141 (MQ=255)
           tGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAgg                                                                                                                >  1:629491/1‑141 (MQ=255)
           tGGTGAAGACGTTGGTGTTAACGGCGGCGTTTTCCGTGTTACTGAAGGACTACAAAAAGAAATTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGAGGTTTAGCGATGGGTCTTGCAGTTGaacg                                                                                                                >  2:677985/1‑139 (MQ=255)
            ggTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGa                                                                                                               <  1:412016/141‑1 (MQ=255)
            ggTGAAGACGTTGGTGTTAACGGCGGTGGTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGa                                                                                                               >  1:338948/1‑141 (MQ=255)
              gtaaGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGAtt                                                                                                             <  2:515928/139‑1 (MQ=255)
              ggAAGACGTTGGTGTTAACGGCGGTGTGTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGAtt                                                                                                             <  2:329509/140‑1 (MQ=255)
               gAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATCCACCTTTACCCGAACCAGGTATGGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATtc                                                                                                            >  1:226395/1‑141 (MQ=255)
               gAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTATCGATGGGTCTTGCAGTTGAAGGATtc                                                                                                            >  1:547652/1‑141 (MQ=255)
                 aggcGTTGGTGTTACCGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAAGTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATtccg                                                                                                          <  1:288031/138‑1 (MQ=255)
                 aGACGTTGGTGTTAACGGCGGTGTTTTAAGTGTTACTCGATGATTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATtccg                                                                                                          <  1:534549/141‑1 (MQ=255)
                              aaCGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGGTGAAGGATTCCGTCCGGTTATGGaa                                                                                             >  2:690195/1‑141 (MQ=255)
                                     gtgtTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAAt                                                                                      >  1:587587/1‑141 (MQ=255)
                                      tgtTTTCCGCGTTACTGAAGGACTCAAAAAGCAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAAtt                                                                                     >  1:510140/1‑141 (MQ=255)
                                        ttttCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGGAGATAGAGTGTTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGACGGATTCCGTCCGGTTATGGGGGTGCAAtttt                                                                                   >  2:2753/1‑139 (MQ=255)
                                        ttttCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTAAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCt                                                                                   >  1:447667/1‑141 (MQ=255)
                                               tgtTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTATCGATGGGTCTTGCAGTTGAAGGCTTCCGTCGGGGTATGGAAGTACAATTCTTGGGttt                                                                            >  1:189228/1‑141 (MQ=255)
                                               tgtTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGttt                                                                            >  1:563537/1‑141 (MQ=255)
                                               tgtTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTGTAGCTGAATCAGGTATTGGTGGTTTATCGATGGGTCTTGCAGTTGAAGGATTCCATCCGGTTATGGAAGTACATTTCTTAGGttt                                                                            >  2:441097/1‑141 (MQ=255)
                                               tgtTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTAGTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGttt                                                                            >  2:348225/1‑141 (MQ=255)
                                                gtTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGAGGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACCATTCTTAGGTTTc                                                                           <  1:538251/141‑1 (MQ=255)
                                                               cAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACCATTCTTAGGTTTCGTCTTCGAAGTAttt                                                            >  2:436409/1‑141 (MQ=255)
                                                                aaaaaGAATTTGGTGAAGATTGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATACTTAGGTTTCGTATTCGAAGTATTTg                                                           >  1:421153/1‑141 (MQ=255)
                                                                aaaaaGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGCATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACACGTCCTAGGTCTCGTATTCGAAGTAGTTg                                                           >  2:693133/1‑141 (MQ=255)
                                                                   aggAATTTGTTGAAGATAGAGTATTCGCTACGCCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATg                                                        <  2:528100/139‑1 (MQ=255)
                                                                   aaGTATTTGGTGAAGATAGAGTATTCGATACACATTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTTCGGTTGTGGTAGTACAATTCTTAGGTTTCGTATTCGAAGTATatgatg                                                        >  2:507519/1‑141 (MQ=255)
                                                                   aaGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACGATTCTTAGGTTTCGTATTCGAAGTATTTGCTg                                                        >  2:28722/1‑141 (MQ=255)
                                                                   aaGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATg                                                        >  1:269715/1‑141 (MQ=255)
                                                                     gAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCg                                                      <  1:603758/141‑1 (MQ=255)
                                                                      aaaTTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCg                                                      <  2:170151/140‑1 (MQ=255)
                                                                       aaTTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGa                                                     <  1:404788/141‑1 (MQ=255)
                                                                       aaTTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGa                                                     <  1:479478/141‑1 (MQ=255)
                                                                          ttGGTGAAGGTAGAGTATTCGATACACCTTTAGCTGAAGCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTCGGTCTCGTATTCGAAGTATTTGATGCGATTg                                                  <  2:667541/141‑1 (MQ=255)
                                                                          ttGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTg                                                  >  1:678635/1‑141 (MQ=255)
                                                                             gTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGGTGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTg                                               >  1:11355/1‑141 (MQ=255)
                                                                             gTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCTAAGTATTTGATGCGATTGCTg                                               >  1:72737/1‑141 (MQ=255)
                                                                             gTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTTATGCGATTGCTg                                               >  2:118225/1‑141 (MQ=255)
                                                                               gAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTg                                               <  1:111636/139‑1 (MQ=255)
                                                                               gAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGa                                             <  1:293728/141‑1 (MQ=255)
                                                                               gAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGa                                             >  2:156227/1‑141 (MQ=255)
                                                                               gAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGGTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTg                                               >  2:111636/1‑139 (MQ=255)
                                                                               gAAGACAGAGGATCCGATACACATTTACTTGAAGCAGGTATTGGTGGTTTAGCGATGGGTCGTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGa                                             <  1:256411/141‑1 (MQ=255)
                                                                                  gATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACaa                                          >  2:184413/1‑141 (MQ=255)
                                                                                     agagTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTCGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAAtt                                       >  1:16347/1‑141 (MQ=255)
                                                                                              gATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACt                              >  1:520110/1‑141 (MQ=255)
                                                                                                tACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCg                            >  2:30412/1‑141 (MQ=255)
                                                                                                  cacCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGtt                          <  2:269715/141‑1 (MQ=255)
                                                                                                    ccTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTc                        >  2:524455/1‑141 (MQ=255)
                                                                                                            tGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTGGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAACTTGCACGTACTCGTTTCCGTTCAgg                >  2:587431/1‑141 (MQ=255)
                                                                                                              aagaaGGTATTGGTGGTTTAGCGATCGGTCTTGCATTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCg              <  2:99668/137‑1 (MQ=255)
                                                                                                                  aGGTATTGGTGGTTTATCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCGGTAc          >  2:165551/1‑141 (MQ=255)
                                                                                                                   ttcagctgttgtTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTAGTCGTTTCCGTTCAGGGGGTACt         >  2:180645/11‑141 (MQ=255)
                                                                                                                   ggTATTGGTGGTTTAGCGATGGGTCTTGGAGTTGAGGGATTCCGTCCGGTTATGGAAGTACAATTCTGAGGTTTCGTATTCGAAGTGGTTGAGGCGTTTGCTGAACAAATAGCACGTACTCGTTTCCGTACGGGCGCTAca         >  2:301237/1‑140 (MQ=255)
                                                                                                                    ggatgggtggtTTATCGATTGGTTTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATGCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCGGTACTa        <  2:17720/136‑1 (MQ=255)
                                                                                                                    ggaTTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTGCGTTCAGGCGGTACTa        <  1:348225/139‑1 (MQ=255)
                                                                                                                    gTATTGGTGGTTTGGCGGTGGGTTTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCGGTACTa        <  2:74298/141‑1 (MQ=255)
                                                                                                                    gTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCGGTACTa        <  2:563537/141‑1 (MQ=255)
                                                                                                                    gTATTGGGGGTTTGGGGTTGGGTTTTGCAGTGGTAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCGGTACTa        <  2:74326/141‑1 (MQ=255)
                                                                                                                       ttGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCGGTACTaaaa     >  2:9498/1‑141 (MQ=255)
                                                                                                                          gtggtTTATCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCGGTACTAAAACTg  <  1:165551/141‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                    
GTTTTAATTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAGAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCGGTACTAAAACTG  >  CP000730/1103664‑1103925

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