Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,137,213 C→T C437Y (TGC→TAC)  gss ← fused glutathionylspermidine amidase/glutathionylspermidine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,137,2130CT98.4% 217.6 / ‑4.2 64C437Y (TGC→TAC) gssfused glutathionylspermidine amidase/glutathionylspermidine synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/35);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

GGGCGCAGCAATACGTCGATAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGCAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTGCAGCGCCTGCTCCATAAACTGCGCGTGATAG  >  NC_000913/3137078‑3137344
                                                                                                                                       |                                                                                                                                   
gggCGCAGCAATACGTCGATAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGtt                                                                                                                                 <  1:161856/140‑1 (MQ=255)
   cgcAGCAATACGTCGATAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAAc                                                                                                                              >  1:924515/1‑140 (MQ=255)
         aaTACGTCGATAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCg                                                                                                                        <  1:898454/140‑1 (MQ=255)
         aaTACGTCGATAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCg                                                                                                                        <  2:1052620/140‑1 (MQ=255)
         aaTACGTCGATAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCg                                                                                                                        <  1:331699/140‑1 (MQ=255)
            aCGTCGATAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTACCCAGTCGccc                                                                                                                     <  2:132880/140‑1 (MQ=255)
                cGATAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAg                                                                                                                           <  2:693363/130‑1 (MQ=255)
                cGATAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAg                                                                                                                           >  1:693363/1‑130 (MQ=255)
                cGATAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTcc                                                                                                                 <  2:382845/140‑1 (MQ=255)
                  aTAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTcccc                                                                                                               >  2:908684/1‑140 (MQ=255)
                  aTAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTcccc                                                                                                               >  1:999314/1‑140 (MQ=255)
                   tAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCa                                                                                                              <  2:383090/140‑1 (MQ=255)
                    aaGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCat                                                                                                             <  1:258740/140‑1 (MQ=255)
                          cACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCatcaatca                                                                                                       >  1:290776/1‑140 (MQ=255)
                            cTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTATCTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGt                                                                                                     >  1:423457/1‑140 (MQ=255)
                                 ttttGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGccc                                                                                                <  2:638799/140‑1 (MQ=255)
                                     gCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGgca                                                                                            >  1:706517/1‑140 (MQ=255)
                                      cGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTATCTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGgcag                                                                                           <  2:423457/140‑1 (MQ=255)
                                      cGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGgcag                                                                                           <  1:908684/140‑1 (MQ=255)
                                      cGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGgcag                                                                                           <  2:924515/140‑1 (MQ=255)
                                           gACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATcc                                                                                      <  1:734192/140‑1 (MQ=255)
                                           gACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATcc                                                                                      >  1:901086/1‑140 (MQ=255)
                                            aCCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATccc                                                                                     >  1:451977/1‑140 (MQ=255)
                                               ggTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATcccagc                                                                                  <  2:352751/140‑1 (MQ=255)
                                               ggTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATcccagc                                                                                  <  2:901086/140‑1 (MQ=255)
                                               ggTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATcccagc                                                                                  <  1:1088722/140‑1 (MQ=255)
                                                gTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATcccagcc                                                                                 >  1:241488/1‑140 (MQ=255)
                                                   cGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATcccagcccag                                                                              <  2:51707/140‑1 (MQ=255)
                                                       ttGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCa                                                                          >  1:382381/1‑140 (MQ=255)
                                                           cACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATcccagc                                                                                  >  1:1026319/1‑128 (MQ=255)
                                                           cACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATcccagc                                                                                  <  2:1026319/128‑1 (MQ=255)
                                                           cACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCatcca                                                                      >  1:417811/1‑140 (MQ=255)
                                                           cACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCatcca                                                                      >  2:629089/1‑140 (MQ=255)
                                                           cACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCatcca                                                                      <  2:621848/140‑1 (MQ=255)
                                                             ccGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCatccaat                                                                    <  2:706517/140‑1 (MQ=255)
                                                             ccGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCatccaat                                                                    <  2:382381/140‑1 (MQ=255)
                                                             ccGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCatccaat                                                                    <  2:373133/140‑1 (MQ=255)
                                                               gcagcaAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCatccaatcc                                                                  >  2:332716/1‑140 (MQ=255)
                                                                cagcaAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCatccaatcca                                                                 >  2:649790/1‑140 (MQ=255)
                                                                  gcaAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACg                                                               <  1:629089/140‑1 (MQ=255)
                                                                  gcaAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACg                                                               <  2:451977/140‑1 (MQ=255)
                                                                  gcaAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACg                                                               <  2:290776/140‑1 (MQ=255)
                                                                  gcaAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCACCCAATCCACg                                                               <  1:671059/140‑1 (MQ=255)
                                                                  gcaAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTACCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACg                                                               <  1:222414/140‑1 (MQ=255)
                                                                  gcaAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGCAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACg                                                               <  1:553681/140‑1 (MQ=255)
                                                                           cGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATAcg                                                      <  2:241488/140‑1 (MQ=255)
                                                                            ggTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATAcgc                                                     <  2:417811/140‑1 (MQ=255)
                                                                                        cACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGccc                                         >  1:749505/1‑140 (MQ=255)
                                                                                        cACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGccc                                         >  2:42154/1‑140 (MQ=255)
                                                                                        cACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGccc                                         >  2:630007/1‑140 (MQ=255)
                                                                                          cGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGc                                       >  1:429875/1‑140 (MQ=255)
                                                                                          cGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGc                                       >  2:1050772/1‑140 (MQ=255)
                                                                                                tGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTg                                 >  2:715935/1‑140 (MQ=255)
                                                                                                tGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTg                                 >  1:962319/1‑140 (MQ=255)
                                                                                                tGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTg                                 >  2:126124/1‑140 (MQ=255)
                                                                                                 gATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTGc                                <  2:108415/140‑1 (MQ=255)
                                                                                                      aaCGCGGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTGCAgcgc                           <  1:630007/140‑1 (MQ=255)
                                                                                                        cgcgGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTGCAGCGCCt                         <  2:749505/140‑1 (MQ=255)
                                                                                                        cgcgGTTTCCCACGCCCAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTGCAGCGCCt                         <  2:101667/140‑1 (MQ=255)
                                                                                                                      cccAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTGCAGCGCCTGCTCCATAAACTgc           <  1:649790/140‑1 (MQ=255)
                                                                                                                       ccAGGTTTTCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTGCAGCGCCTGCTCCATAAACTgcg          >  2:1016016/1‑140 (MQ=255)
                                                                                                                               tCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTGCAGCGCCTGCTCCATAAACTGCGCGTGATAg  >  1:341373/1‑140 (MQ=255)
                                                                                                                               tCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTGCAGCGCCTGCTCCATAAACTGCGCGTGATAg  >  2:540659/1‑140 (MQ=255)
                                                                                                                               tCCACACGTAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTGCAGCGCCTGCTCCATAAACTGCGCGTGATAg  >  1:1097123/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                   
GGGCGCAGCAATACGTCGATAAGACGCACTTCGTTTTGCGGATGACCGGTACGGATTGGCACCGCAGCAAACTCACGGTCGCTAACTTCACGAATCTGATCAAACGCGGTTTCCCACGCCCAGGTTTTCCACACGCAGTTAACCAGTCGCCCTTCCCCATCAATCAGTTGCCCGGCAGCATCCCAGCCCAGTTCATCCAATCCACGCAAGATACGCGTTTCAAAGCCCGCCTGGTGCAGCGCCTGCTCCATAAACTGCGCGTGATAG  >  NC_000913/3137078‑3137344

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