Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,685,910 G→A A26T (GCC→ACC)  hmp → fused nitric oxide dioxygenase/dihydropteridine reductase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,685,9100GA100.0% 265.4 / NA 81A26T (GCC→ACC) hmpfused nitric oxide dioxygenase/dihydropteridine reductase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/43);  total (38/43)

CTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCGCCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAG  >  NC_000913/2685775‑2686031
                                                                                                                                       |                                                                                                                         
cTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                        >  1:200317/1‑139 (MQ=255)
  aTTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCAt                                                                                                                      <  2:200317/139‑1 (MQ=255)
            aGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAcc                                                                                                            <  2:197767/139‑1 (MQ=255)
             gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             <  1:717025/137‑1 (MQ=255)
             gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             >  2:717025/1‑137 (MQ=255)
             gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                           <  1:1108320/139‑1 (MQ=255)
             gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                           <  2:773767/139‑1 (MQ=255)
             gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                           >  2:192258/1‑139 (MQ=255)
              aTGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGt                                                                                                          <  2:1335499/139‑1 (MQ=255)
                         aTACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                               >  1:75029/1‑139 (MQ=255)
                         aTACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                               >  1:1462639/1‑139 (MQ=255)
                             aTCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                           >  1:1443748/1‑123 (MQ=255)
                             aTCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                           <  2:1443748/123‑1 (MQ=255)
                                       aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGTCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                        <  1:1041409/100‑1 (MQ=255)
                                       aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                        <  2:1069547/100‑1 (MQ=255)
                                       aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                        >  2:1041409/1‑100 (MQ=255)
                                       aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                        >  1:293912/1‑100 (MQ=255)
                                       aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                        <  2:293912/100‑1 (MQ=255)
                                       aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                        >  1:1069547/1‑100 (MQ=255)
                                       aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCa                                                                                 >  2:1400800/1‑139 (MQ=255)
                                       aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCa                                                                                 >  2:1292770/1‑139 (MQ=255)
                                       aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCa                                                                                 >  1:1308040/1‑139 (MQ=255)
                                          cAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAg                                                                              >  2:579928/1‑139 (MQ=255)
                                          cAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAg                                                                              >  2:518271/1‑139 (MQ=255)
                                           aaaaaaaGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGa                                                                             <  2:771610/139‑1 (MQ=255)
                                                   gAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTa                                                                     <  2:1462639/139‑1 (MQ=255)
                                                              gCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAGATTTTTAACATGAGTAAc                                                          <  2:771647/139‑1 (MQ=255)
                                                              gCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                          <  2:1193389/139‑1 (MQ=255)
                                                              gCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                          <  1:192258/139‑1 (MQ=255)
                                                              gCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                          <  2:75029/139‑1 (MQ=255)
                                                              gCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                          <  1:941127/139‑1 (MQ=255)
                                                                 tGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                 <  1:75044/81‑1 (MQ=255)
                                                                 tGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                 >  2:75044/1‑81 (MQ=255)
                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                               >  1:568849/1‑97 (MQ=255)
                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                               <  2:568849/97‑1 (MQ=255)
                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATa                                                                                             <  1:886169/99‑1 (MQ=255)
                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCg                                                     >  2:195331/1‑139 (MQ=255)
                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTGCTACGACCGTATGTTTACTCATa                                                                                             >  2:886169/1‑99 (MQ=255)
                                                                         aaaCCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTAc                                                                                                                >  2:1315844/1‑74 (MQ=255)
                                                                         aaaCCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTAc                                                                                                                <  1:1315844/74‑1 (MQ=255)
                                                                          aaCCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAAtt                                                                         >  1:640448/1‑112 (MQ=255)
                                                                          aaCCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAAtt                                                                         <  2:640448/112‑1 (MQ=255)
                                                                                cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             <  2:382917/70‑1 (MQ=255)
                                                                                cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             >  1:382917/1‑70 (MQ=255)
                                                                                cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGttt                                                                                                    >  1:840871/1‑79 (MQ=255)
                                                                                cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGttt                                                                                                    <  2:840871/79‑1 (MQ=255)
                                                                                cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                               <  1:508475/84‑1 (MQ=255)
                                                                                cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                               >  2:508475/1‑84 (MQ=255)
                                                                                cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCaa                                        >  2:54919/1‑139 (MQ=255)
                                                                                        tAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                        <  1:340512/51‑1 (MQ=255)
                                                                                        tAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                        >  2:340512/1‑51 (MQ=255)
                                                                                         aaaaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGcc                               <  1:518271/139‑1 (MQ=255)
                                                                                         aaaaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGcc                               <  1:1400800/139‑1 (MQ=255)
                                                                                         aaaaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGATAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGcc                               <  2:266577/139‑1 (MQ=255)
                                                                                             gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             >  1:246642/1‑57 (MQ=255)
                                                                                             gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             <  2:1364504/57‑1 (MQ=255)
                                                                                             gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             <  2:246642/57‑1 (MQ=255)
                                                                                             gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             >  1:1364504/1‑57 (MQ=255)
                                                                                               caccaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             >  2:738390/1‑55 (MQ=255)
                                                                                               caccaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             <  1:738390/55‑1 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                           >  1:810326/1‑54 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                           <  2:810326/54‑1 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                    >  2:264508/1‑93 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                    <  1:264508/93‑1 (MQ=255)
                                                                                                   aTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACg                     <  2:1138164/139‑1 (MQ=255)
                                                                                                         ttACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTg               >  1:711976/1‑139 (MQ=255)
                                                                                                         ttACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTg               <  1:195331/139‑1 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             >  2:659636/1‑39 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                             <  1:659636/39‑1 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                           >  2:413207/1‑41 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                           <  1:413207/41‑1 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTg                                    <  2:383462/112‑1 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTg                                    >  1:383462/1‑112 (MQ=255)
                                                                                                                tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACTGTATGTTTACTCATAACCCAGAGCTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTTATGCCCTGTTTAACGCTATTGCCGCCTa        <  1:54919/139‑1 (MQ=255)
                                                                                                                 ggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTAc       >  1:59206/1‑139 (MQ=255)
                                                                                                                  gAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACg      >  1:189205/1‑139 (MQ=255)
                                                                                                                      cGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAg  <  2:59206/139‑1 (MQ=255)
                                                                                                                       ggggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccg            >  1:942505/1‑128 (MQ=255)
                                                                                                                       ggggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccg            <  2:942505/128‑1 (MQ=255)
                                                                                                                         ggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgc           >  2:1255535/1‑127 (MQ=255)
                                                                                                                         ggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgc           <  1:1255535/127‑1 (MQ=255)
                                                                                                                                       |                                                                                                                         
CTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCGCCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAG  >  NC_000913/2685775‑2686031

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