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

GTTTTAATTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCGGTACTAAAACTG  >  CP000730/1103664‑1103925
                                                                                                                                  |                                                                                                                                    
tttttAAGTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGAAAGAGTATTCGATACACCTTTAGCTGAAGAATGTATTGGTGGTTTAGCGATGGGTCtt                                                                                                                           <  2:635395/140‑1 (MQ=255)
gTTTTAATTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCTATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCtt                                                                                                                           <  2:450453/141‑1 (MQ=255)
gTTTTAATTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTGGCTGAATCAGGTATGGGTGGTTTAGCGATGGGTCtg                                                                                                                           >  1:293997/1‑140 (MQ=255)
 ttttAATTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGGAGTACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTGTTGGTGGTTTAGCGGTGGGTCTTg                                                                                                                          >  2:158874/1‑141 (MQ=255)
    tAAGTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAg                                                                                                                       <  2:469551/141‑1 (MQ=255)
    ataTTTTTGGTGAAGACGTGGGTGCTAACGGGGGTGTTTACCGAGATCATGAAGGACTACAAAAAG‑CCTTCGGTGAAGATCGAGTACTCGATACACCTTTCGCTGAATCAGGTATTGGTGCTTTAGCGATGGGTCTTGCAg                                                                                                                       <  1:487395/139‑1 (MQ=255)
       tttttGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTg                                                                                                                    >  2:648029/1‑141 (MQ=255)
          ttGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAg                                                                                                                 >  2:89606/1‑141 (MQ=255)
          ttGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAg                                                                                                                 <  1:429192/141‑1 (MQ=255)
          ttGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGAGGGGTCTTGCAGtgtacg                                                                                                                 >  2:226739/1‑136 (MQ=255)
           tggttgaGACGTTGGTGTTAACTGCGGTGTTTTCCGTGTTACTGAAGGAGTGAGAAATG‑GGTTTGGTGATGATAGAGTATTCGATGCACCTTTAGTTAATTCAGGTATTGGTGGTTGAGCGGTGTGTGTTGCGGTTGGAgg                                                                                                                >  2:534526/7‑141 (MQ=255)
           tGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAgg                                                                                                                >  1:629491/1‑141 (MQ=255)
           tGGTGAAGACGTTGGTGTTAACGGCGGCGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AAATTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGAGGTTTAGCGATGGGTCTTGCAGTTGaacg                                                                                                                >  2:677985/1‑139 (MQ=255)
            ggTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGa                                                                                                               <  1:412016/141‑1 (MQ=255)
            ggTGAAGACGTTGGTGTTAACGGCGGTGGTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGa                                                                                                               >  1:338948/1‑141 (MQ=255)
              gtaaGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGAtt                                                                                                             <  2:515928/139‑1 (MQ=255)
              ggAAGACGTTGGTGTTAACGGCGGTGTGTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGAtt                                                                                                             <  2:329509/140‑1 (MQ=255)
               gAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATCCACCTTTACCCGAACCAGGTATGGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATtc                                                                                                            >  1:226395/1‑141 (MQ=255)
               gAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTATCGATGGGTCTTGCAGTTGAAGGATtc                                                                                                            >  1:547652/1‑141 (MQ=255)
                 aggcGTTGGTGTTACCGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AAGTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATtccg                                                                                                          <  1:288031/138‑1 (MQ=255)
                 aGACGTTGGTGTTAACGGCGGTGTTTTAAGTGTTACTCGATGATTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATtccg                                                                                                          <  1:534549/141‑1 (MQ=255)
                              aaCGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGGTGAAGGATTCCGTCCGGTTATGGaa                                                                                             >  2:690195/1‑141 (MQ=255)
                                     gtgtTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAAt                                                                                      >  1:587587/1‑141 (MQ=255)
                                      tgtTTTCCGCGTTACTGAAGGACT‑CAAAAAGCAATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAAtt                                                                                     >  1:510140/1‑141 (MQ=255)
                                        ttttCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGGAGATAGAGTGTTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGACGGATTCCGTCCGGTTATGGGGGTGCAAtttt                                                                                   >  2:2753/1‑139 (MQ=255)
                                        ttttCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTAAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCt                                                                                   >  1:447667/1‑141 (MQ=255)
                                               tgtTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTATCGATGGGTCTTGCAGTTGAAGGCTTCCGTCGGGGTATGGAAGTACAATTCTTGGGttt                                                                            >  1:189228/1‑141 (MQ=255)
                                               tgtTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGttt                                                                            >  1:563537/1‑141 (MQ=255)
                                               tgtTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTGTAGCTGAATCAGGTATTGGTGGTTTATCGATGGGTCTTGCAGTTGAAGGATTCCATCCGGTTATGGAAGTACATTTCTTAGGttt                                                                            >  2:441097/1‑141 (MQ=255)
                                               tgtTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTAGTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGttt                                                                            >  2:348225/1‑141 (MQ=255)
                                                gtTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGAGGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACCATTCTTAGGTTTc                                                                           <  1:538251/141‑1 (MQ=255)
                                                               cAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACCATTCTTAGGTTTCGTCTTCGAAGTAttt                                                            >  2:436409/1‑141 (MQ=255)
                                                                aaaaaG‑AATTTGGTGAAGATTGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATACTTAGGTTTCGTATTCGAAGTATTTg                                                           >  1:421153/1‑141 (MQ=255)
                                                                aaaaaG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGCATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACACGTCCTAGGTCTCGTATTCGAAGTAGTTg                                                           >  2:693133/1‑141 (MQ=255)
                                                                   agg‑AATTTGTTGAAGATAGAGTATTCGCTACGCCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATg                                                        <  2:528100/139‑1 (MQ=255)
                                                                   aaG‑TATTTGGTGAAGATAGAGTATTCGATACACATTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTTCGGTTGTGGTAGTACAATTCTTAGGTTTCGTATTCGAAGTATatgatg                                                        >  2:507519/1‑141 (MQ=255)
                                                                   aaG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACGATTCTTAGGTTTCGTATTCGAAGTATTTGCTg                                                        >  2:28722/1‑141 (MQ=255)
                                                                   aaG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATg                                                        >  1:269715/1‑141 (MQ=255)
                                                                     g‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCg                                                      <  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)
                                                                                                                                  |                                                                                                                                    
GTTTTAATTTTTGGTGAAGACGTTGGTGTTAACGGCGGTGTTTTCCGTGTTACTGAAGGACTACAAAAAG‑AATTTGGTGAAGATAGAGTATTCGATACACCTTTAGCTGAATCAGGTATTGGTGGTTTAGCGATGGGTCTTGCAGTTGAAGGATTCCGTCCGGTTATGGAAGTACAATTCTTAGGTTTCGTATTCGAAGTATTTGATGCGATTGCTGGACAAATTGCACGTACTCGTTTCCGTTCAGGCGGTACTAAAACTG  >  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: ‑