Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 138,784 C→T 9.5% intergenic (+151/+51) 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,7840CT 9.5% 193.2 / 9.5 63intergenic (+151/+51)cueO/gcdmulticopper oxidase (laccase)/glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base T (4/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

AGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGC  >  NC_000913/138636‑138931
                                                                                                                                                    |                                                                                                                                                   
aGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCAc                                                                                                                                                     >  1:264467/1‑149 (MQ=255)
        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                                                                                                                                >  3:539927/1‑149 (MQ=255)
                     ccGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                                                                                                                <  4:393805/149‑1 (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                                                                                               <  4:539927/149‑1 (MQ=255)
                                                      cGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTAc                                                                                               <  4:435144/149‑1 (MQ=255)
                                                      cGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTAc                                                                                               <  1:390426/149‑1 (MQ=255)
                                                       gCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGTATCCGGCATGAACAACGCGCACGTTGTCAACAATCTGAg                                                                                                         >  2:513665/1‑138 (MQ=255)
                                                       gCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGa                                                                                                          >  1:600121/1‑137 (MQ=255)
                                                       gCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGa                                                                                                          <  2:600121/137‑1 (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                                                                      <  4:587781/149‑1 (MQ=255)
                                                                               ttATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAg                                                                      <  1:503122/149‑1 (MQ=255)
                                                                                tATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGc                                                                     <  4:51119/149‑1 (MQ=255)
                                                                                  tCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca                                                                   <  2:437681/149‑1 (MQ=255)
                                                                                  tCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca                                                                   >  2:397124/1‑149 (MQ=255)
                                                                                  tCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca                                                                   >  4:8686/1‑149 (MQ=255)
                                                                                  tCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca                                                                   >  2:161335/1‑149 (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                            >  2:490300/1‑149 (MQ=255)
                                                                                                                         aaTTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGAc                            >  3:391219/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                                 <  2:608058/121‑1 (MQ=39)
                                                                                                                                                ttCATGCCGGATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAAcc                                 >  3:93900/1‑121 (MQ=39)
                                                                                                                                                ttCATGCCGGATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAAcc                                 >  1:608058/1‑121 (MQ=39)
                                                                                                                                                ttCATGCCGGATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAAcc                                 <  4:93900/121‑1 (MQ=39)
                                                                                                                                                   aCGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGc  <  2:212974/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGC  >  NC_000913/138636‑138931

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