Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 138,789 C→G 9.6% intergenic (+156/+46) cueO → / ← gcd multicopper oxidase (laccase)/glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913138,7890CG 9.6% 208.7 / 12.3 63intergenic (+156/+46)cueO/gcdmulticopper oxidase (laccase)/glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base G (4/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

TTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTT  >  NC_000913/138644‑138933
                                                                                                                                                 |                                                                                                                                                
ttCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATc                                                                                                                                               <  4:594599/149‑1 (MQ=255)
   ggAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                            >  1:437681/1‑149 (MQ=255)
      aaTATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                         >  2:390426/1‑149 (MQ=255)
         atGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaa                                                                                                                                      <  2:69462/149‑1 (MQ=255)
             ccGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                                                                                                                  <  4:393805/149‑1 (MQ=255)
             ccGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                                                                                                                  >  3:539927/1‑149 (MQ=255)
              cGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg                                                                                                                                 >  3:572278/1‑149 (MQ=255)
                    tctcCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACg                                                                                                                           >  3:51119/1‑149 (MQ=255)
                      tcCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGtt                                                                                                                         >  1:610286/1‑149 (MQ=255)
                        cGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGt                                                                                                                       <  4:272645/149‑1 (MQ=255)
                          ggCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCa                                                                                                                     <  2:352087/149‑1 (MQ=255)
                           gCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAg                                                                                                                    <  2:336922/149‑1 (MQ=255)
                             aTTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCa                                                                                                                  >  2:70005/1‑149 (MQ=255)
                                   cTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGa                                                                                                            <  1:149464/149‑1 (MQ=255)
                                   cTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGa                                                                                                            <  4:512603/149‑1 (MQ=255)
                                    tgatgaCAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAg                                                                                                           >  2:503122/1‑149 (MQ=255)
                                        gaCAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAg                                                                                                       >  4:275433/1‑149 (MQ=255)
                                        gaCAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAg                                                                                                       >  2:555911/1‑149 (MQ=255)
                                             aCGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTa                                                                                                  <  2:318147/149‑1 (MQ=255)
                                              cGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTAc                                                                                                 <  1:390426/149‑1 (MQ=255)
                                              cGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTAc                                                                                                 <  4:435144/149‑1 (MQ=255)
                                              cGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTAc                                                                                                 <  4:539927/149‑1 (MQ=255)
                                               gCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGTATCCGGCATGAACAACGCGCACGTTGTCAACAATCTGAg                                                                                                           >  2:513665/1‑138 (MQ=255)
                                               gCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGa                                                                                                            <  2:600121/137‑1 (MQ=255)
                                               gCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGa                                                                                                            >  1:600121/1‑137 (MQ=255)
                                               gCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACt                                                                                                <  1:555911/149‑1 (MQ=255)
                                               gCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAACAATCTGAg                                                                                                           <  1:513665/138‑1 (MQ=255)
                                                 aaaTCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTc                                                                                              <  4:166791/149‑1 (MQ=255)
                                                      tGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatc                                                                                         <  1:4262/149‑1 (MQ=255)
                                                                       ttATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTa                                                                                                  >  4:219286/1‑123 (MQ=255)
                                                                       ttATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTa                                                                                                  <  3:219286/123‑1 (MQ=255)
                                                                       ttATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAg                                                                        <  1:503122/149‑1 (MQ=255)
                                                                       ttATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAg                                                                        <  4:587781/149‑1 (MQ=255)
                                                                        tATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGc                                                                       <  4:51119/149‑1 (MQ=255)
                                                                          tCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca                                                                     >  2:161335/1‑149 (MQ=255)
                                                                          tCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca                                                                     >  4:8686/1‑149 (MQ=255)
                                                                          tCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca                                                                     >  2:397124/1‑149 (MQ=255)
                                                                          tCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca                                                                     <  2:437681/149‑1 (MQ=255)
                                                                              gCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAat                                                                 >  2:613829/1‑149 (MQ=255)
                                                                                   cGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGt                                                            >  1:212974/1‑149 (MQ=255)
                                                                                       aGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGcc                                                        >  2:195336/1‑149 (MQ=255)
                                                                                       aGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGcc                                                        <  1:195336/149‑1 (MQ=255)
                                                                                        gCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGccc                                                       >  2:571652/1‑149 (MQ=255)
                                                                                                aTCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcaca                                                                                           >  3:146261/1‑105 (MQ=255)
                                                                                                aTCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcaca                                                                                           <  4:146261/105‑1 (MQ=255)
                                                                                                   cATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTAcc                                            >  1:106300/1‑149 (MQ=255)
                                                                                                    aTTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCa                                           >  4:457669/1‑149 (MQ=255)
                                                                                                              cacaATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGt                                 <  1:613829/149‑1 (MQ=255)
                                                                                                               acaATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTg                                >  2:73032/1‑149 (MQ=255)
                                                                                                                 aaTTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGAc                              >  3:391219/1‑149 (MQ=255)
                                                                                                                 aaTTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGAc                              >  2:490300/1‑149 (MQ=255)
                                                                                                                   ttttGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCg                            >  4:233729/1‑149 (MQ=255)
                                                                                                                             tGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGgaga                  >  2:221730/1‑149 (MQ=255)
                                                                                                                               gATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATc                <  1:70005/149‑1 (MQ=255)
                                                                                                                                   cgttgTTCATGCCGGATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATcacca            >  4:358623/5‑149 (MQ=255)
                                                                                                                                   aGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATcacca            <  2:610286/149‑1 (MQ=255)
                                                                                                                                    gttgTTCATGCCGGATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATcaccac           >  1:29776/4‑149 (MQ=39)
                                                                                                                                    ggCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATcaccac           >  2:10952/1‑149 (MQ=255)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAAcc                                   <  4:93900/121‑1 (MQ=39)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAAcc                                   >  3:93900/1‑121 (MQ=39)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAAcc                                   >  1:608058/1‑121 (MQ=39)
                                                                                                                                        ttCATGCCGGATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAAcc                                   <  2:608058/121‑1 (MQ=39)
                                                                                                                                          catgccggATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGgag                   <  4:108064/132‑1 (MQ=255)
                                                                                                                                          catgccggATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGgag                   >  3:108064/4‑135 (MQ=255)
                                                                                                                                           aCGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGc    <  2:212974/149‑1 (MQ=255)
                                                                                                                                             gccgcATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATc                                                    <  1:306331/99‑1 (MQ=255)
                                                                                                                                             gccgcATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATc                                                    >  2:306331/1‑99 (MQ=255)
                                                                                                                                             gccgcATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCtt  >  4:283149/1‑149 (MQ=255)
                                                                                                                                             gccgcATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCtt  >  3:432058/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTT  >  NC_000913/138644‑138933

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