Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,138,011 C→A 10.2% G342V (GGC→GTC)  gldA ← L‑1,2‑propanediol dehydrogenase/glycerol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,138,0110CA10.2% 116.1 / 7.1 53G342V (GGC→GTC) gldAL‑1,2‑propanediol dehydrogenase/glycerol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (29/17);  new base A (2/3);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

AAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACGG  >  NC_000913/4137869‑4138152
                                                                                                                                               |                                                                                                                                             
aaGCGTAGCGCATCAGGCAATTTTGCGTGCAAAGTACCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGGCCTGGTCGGCTACCAGCAGAGAGGGGGGGGCCGGGGCTG‑GCGTCgcgc                                                                                                                                               >  2:222551/1‑143 (MQ=255)
 aGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACACCCCGGCAGTTTGGAGTTGTTTATTCCCACTCTTCCCGGAAACGCTGACCGTACCGGTCCGCTACCAGCAGAGCGGCGTAAACCTACTCTG‑TCCTCCcgcc                                                                                                                                              >  1:639352/1‑143 (MQ=255)
   cGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACGCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGAGCGG‑GCGTCGcgctgc                                                                                                                                            >  2:640779/1‑143 (MQ=255)
   cGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGCTTGGAGTAGGTTATTCACACTATCGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGcgccgc                                                                                                                                            >  1:441146/1‑143 (MQ=255)
   cGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGAACGAGCCACGAGATTGCAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGAGCGG‑GCGGCGcgcccc                                                                                                                                            >  2:607139/1‑141 (MQ=255)
    gTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCgccgcc                                                                                                                                           >  1:213950/1‑143 (MQ=255)
     tAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGCCCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCa                                                                                                                                          <  1:853253/143‑1 (MQ=255)
      aGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGGCGCGACGCCAg                                                                                                                                         <  2:274585/143‑1 (MQ=255)
         gcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGGCGCGACGCCAgatc                                                                                                                                      <  2:581112/143‑4 (MQ=255)
         gcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGGCGCGACGCCAgatc                                                                                                                                      <  1:232835/143‑4 (MQ=255)
          cATCAGGCAATTTTGCGTTCAAGCTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCAt                                                                                                                                     <  2:905365/143‑1 (MQ=255)
               ggCAATTTTGTGTGCTGACTCCCGGACGAGCCGGGAGTTTGGAGTATGTTATTCCCTCTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGCGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCAtgttgt                                                                                                                                <  2:123523/143‑1 (MQ=255)
               ggCAATTTTGCGTTCAAACTCCCGGACAAGTCGGGAGTTTGTAGTAGGATATTCCCACACTTGCCGGCAACGCTGACCGTACTGGTCGGCTACCCGCAGAGCGGCGTAACCCTGATCTG‑GCGTCGCGCCCCCCGGccggctga                                                                                                                                >  2:171671/1‑136 (MQ=255)
               ggCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAGACGCTGACCGTACTGGTCGGCTACAAGCAGAGCGGCGTAAACTCGATCTGCGCG‑CGCGCCGCCAGGCAtgtggc                                                                                                                                >  2:142760/1‑140 (MQ=255)
                gCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGCGGTTGGAGTAGGTTATTCCCACTCTTGCAGGAAAAGCTGACCGCAATGGGCGGCTACCACGCGAGAGGCGGAAACCTCATCTG‑GGGTCGCGCACCCAGGAGTGTtgtg                                                                                                                               >  2:869358/1‑143 (MQ=255)
                gCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTtgtg                                                                                                                               <  1:785591/143‑1 (MQ=255)
                  aaTTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTCAACCTGATCTG‑GCGTCGCGGCGCCAGGCACCTTGTGaa                                                                                                                             >  2:64009/1‑143 (MQ=255)
                  aaTTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGGTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGGATGTTGTGaa                                                                                                                             >  1:234563/1‑143 (MQ=255)
                   aTTTTGCGTTCAATCTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCGCTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGGGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAAt                                                                                                                            <  2:913157/143‑1 (MQ=255)
                    ttttGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCCGCAGAGCGGCGTAAACCTGCTCTG‑GCGTCGCGCCGCCCGGCATGTcgtgcacg                                                                                                                           >  2:517732/1‑139 (MQ=255)
                         cGTTCGAGCTCCCGGACAAGACGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTc                                                                                                                      <  2:360321/143‑1 (MQ=255)
                          gTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCa                                                                                                                     >  1:221777/1‑143 (MQ=255)
                               aaCTCCGGGAGAGGCCGGGCGCTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTc                                                                                                                <  1:762020/143‑1 (MQ=255)
                                  tCCCGGACAAGCCGGGCGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGc                                                                                                             >  1:651075/1‑143 (MQ=255)
                                    ccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGGGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGcac                                                                                                           >  1:242913/1‑143 (MQ=255)
                                       gACAAGCCGGGAGTATGGAGTAGGTTATTCCCACTCTTGCCCGAAACGCTGACCGTACTGGCCGGCTACCAGCAGAGCGGCGTAAACCTGATCTC‑GCGTCGCGCCGCCAAGCATGCTGTGACTGGTGGCACCTGCTGCACATg                                                                                                        >  2:423560/1‑143 (MQ=255)
                                                    tttGGAGTAGGTTATTCCCAGTCTTGCACGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCcac                                                                                           <  1:7673/143‑1 (MQ=255)
                                                    tttGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTAACCTTCTGCACATGCCGCTTCTGCcac                                                                                           >  2:100664/1‑143 (MQ=255)
                                                     ttGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGCAAACCTGACCTG‑GCGTCTCGCCGCCGGGAGTGCTCCCAATCGTGTCACCCTCCGCACACCCCGCGCCTGccccc                                                                                          >  2:647325/1‑140 (MQ=255)
                                                       ggAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAAt                                                                                        >  1:732299/1‑143 (MQ=255)
                                                         aGTAGGTTATTCCCACTCTTGCAGGAGTCGCGTTCCGTGCTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGATTCACCTTCTGCGCATGCCGCTTCTGCCACAATTc                                                                                      >  2:113109/1‑143 (MQ=255)
                                                         aGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTCTCACCTTCTGCACATGCCGCCTCTGCCACAACTc                                                                                      >  2:118334/1‑143 (MQ=255)
                                                           tAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAAGCTGGTCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGGTTCGCCTTCTGCACATGCCGCTTCTGCCACAATTCGc                                                                                    >  1:502109/1‑143 (MQ=255)
                                                                 gttCCCGCTCTGGCAGGAAACGCTGTCCGTACTGGTCGGCTACCAGCAGAGAGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGGATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTc                                                                              <  2:179528/142‑1 (MQ=255)
                                                                    cccACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGcc                                                                           <  1:47345/143‑1 (MQ=255)
                                                                     ccACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATGCGCATTTTCGCCg                                                                          >  1:484826/1‑143 (MQ=255)
                                                                      cACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGACTGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCgg                                                                         >  1:798066/1‑143 (MQ=255)
                                                                       aCTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGCTCTG‑GCGTCGCGCCGCCCGGCATGTTCTGAATGCTTTCACCTTCTGCACATGCCACTTCTGCCACAATTCGCATTTTCAccagc                                                                        >  1:354295/1‑140 (MQ=255)
                                                                                 gAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGCCATcttctt                                                              >  2:269682/1‑143 (MQ=255)
                                                                                   aaCGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTa                                                            >  1:467796/1‑143 (MQ=255)
                                                                                   aaCGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCATAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTa                                                            >  1:39016/1‑143 (MQ=255)
                                                                                      gCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCCCAATTCGCATTTTCGCCGGGACCTCTTCCTTActa                                                         >  2:760665/1‑143 (MQ=255)
                                                                                        tGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATc                                                       <  2:484826/143‑1 (MQ=255)
                                                                                             gTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACACCTTCTTTAATATCCCGtt                                                  >  1:750056/1‑143 (MQ=255)
                                                                                                    tCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCgag                                           >  1:723685/1‑143 (MQ=255)
                                                                                                      ggCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGTACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACACCTTCTTTAATATCCAGTTGAGCGagag                                         >  1:162144/1‑143 (MQ=255)
                                                                                                         gcgccgcagAGCAGCGGCAACCTGATCTA‑GCGTTGTGTCGCCAGTCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTa                                      <  2:888099/138‑1 (MQ=255)
                                                                                                            cagcagAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTg                                   <  2:675698/143‑1 (MQ=255)
                                                                                                              gcagAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGCGTTATTGGc                                 >  1:87291/1‑143 (MQ=255)
                                                                                                              gcagAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGc                                 <  2:147261/143‑1 (MQ=255)
                                                                                                              gcagAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGc                                 <  1:278708/143‑1 (MQ=255)
                                                                                                              gcagAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTCCTTTACTATCCAGCTGAGCGAGAGTTATTgga                                 >  2:685079/1‑142 (MQ=255)
                                                                                                                  agCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGTCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAAc                             <  1:618263/143‑1 (MQ=255)
                                                                                                                   gcggcgTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTGCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAGTCTCCAGTTGAGCGAGAGTTATTGGCAAAcc                            >  2:40480/1‑143 (MQ=255)
                                                                                                                       cgTAAACCTGATCTG‑GCGTCGGGCGGGCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTAcc                        <  2:221777/143‑1 (MQ=255)
                                                                                                                                  ccTG‑GCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAg             >  1:779556/2‑143 (MQ=255)
                                                                                                                                   cTG‑GCGGCGCGACGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAgg            >  1:839798/1‑143 (MQ=255)
                                                                                                                                       gccgcgcgCCACCATGCATGTTGTGTTTGGTTTGACCTGCTGCACATGCCGCTTCTGCCACAATTCGCAGTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCa         <  2:267042/139‑1 (MQ=255)
                                                                                                                                         gTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCCAGGGGCAGc       >  1:335280/1‑143 (MQ=255)
                                                                                                                                            gcgcCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACAGCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTAc    <  1:88133/143‑1 (MQ=255)
                                                                                                                                             cgccgcCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGACCATTATGCGCAAACCTACCGCATGGCTAAGGGCAGCTACg   >  1:865670/1‑143 (MQ=255)
                                                                                                                                             cgccgcCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAAAAAGGGGTCGAGGCAGTCCTACCGCACGGCTAAGGGCAGCTaag   >  1:381834/1‑141 (MQ=255)
                                                                                                                                              gccgccTGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTGTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACgg  <  2:558528/143‑1 (MQ=255)
                                                                                                                                              gccgccAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCTAACCTACCGCGTGGCTACGGGCTGCTGCgt  >  1:304508/1‑142 (MQ=255)
                                                                                                                                              gccgccAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACgg  <  2:213950/143‑1 (MQ=255)
                                                                                                                                              gccgccAGCCATGTTGTGGATGGTGTCCCTTTCGGCGCAGGCCGTTTCTGCCACAGTTGGCATTTTCGCCGGTACATCTTCTTTACTATCCAGGTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACgg  <  2:892183/143‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
AAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTG‑GCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTGCCACAATTCGCATTTTCGCCGGGACATCTTCTTTAATATCCAGTTGAGCGAGAGTTATTGGCAAACCTACCGCATGGCTAAGGGCAGCTACGG  >  NC_000913/4137869‑4138152

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