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% 321.9 / NA 98A26T (GCC→ACC) hmpfused nitric oxide dioxygenase/dihydropteridine reductase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (48/50);  total (48/50)

CTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCGCCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTA  >  NC_000913/2685775‑2686033
                                                                                                                                       |                                                                                                                           
cTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          >  1:667799/1‑139 (MQ=255)
       ttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                   <  2:124882/139‑1 (MQ=255)
       ttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                   <  1:65138/139‑1 (MQ=255)
        tataAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTAc                                                                                                                  <  2:1414864/139‑1 (MQ=255)
           aaGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                               <  1:1407608/139‑1 (MQ=255)
                   tttGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGtt                                                                                                       <  1:540375/139‑1 (MQ=255)
                           aCATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATa                                                                                               >  2:507354/1‑139 (MQ=255)
                             aTCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          >  2:552245/1‑110 (MQ=255)
                             aTCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          <  1:552245/110‑1 (MQ=255)
                             aTCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCt                                                                                                                    <  2:1531715/116‑1 (MQ=255)
                             aTCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCt                                                                                                                    >  1:1531715/1‑116 (MQ=255)
                              tCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAAcc                                                                                            >  1:1620227/1‑139 (MQ=255)
                                 aTTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAg                                                                                         >  1:385138/1‑139 (MQ=255)
                                       aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCa                                                                                   >  2:769680/1‑139 (MQ=255)
                                           aaaaaaaGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGa                                                                               <  2:667799/139‑1 (MQ=255)
                                           aaaaaaaGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGa                                                                               <  1:507354/139‑1 (MQ=255)
                                                aaGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAAttt                                                                          >  1:1022021/1‑139 (MQ=255)
                                                  ggAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAAttttt                                                                        <  2:314771/139‑1 (MQ=255)
                                                   gAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAGTTTTTa                                                                       >  2:1543033/1‑139 (MQ=255)
                                                    aaGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAGTTTTTaa                                                                      <  1:1543033/139‑1 (MQ=255)
                                                    aaGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                      <  2:1620227/139‑1 (MQ=255)
                                                       aCCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACAt                                                                   >  2:1681337/1‑139 (MQ=255)
                                                           tatGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          >  2:1488157/1‑80 (MQ=255)
                                                           tatGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          <  1:1488157/80‑1 (MQ=255)
                                                           tatGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGt                                                               >  2:1140087/1‑139 (MQ=255)
                                                              gCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                            <  1:1140087/139‑1 (MQ=255)
                                                                  gACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          <  1:270590/73‑1 (MQ=255)
                                                                  gACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          >  2:270590/1‑73 (MQ=255)
                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCGTTTCTACGACCGTATGTTTACTCATATCCCAGAACTCAAAGAAATTTTTaa                                                                      >  1:1052047/1‑124 (MQ=255)
                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                   <  1:1716751/79‑1 (MQ=255)
                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                   >  2:1716751/1‑79 (MQ=255)
                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTAc                                                                                                                  <  1:305538/80‑1 (MQ=255)
                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTAc                                                                                                                  >  2:305538/1‑80 (MQ=255)
                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGATGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                      <  2:1052047/124‑1 (MQ=255)
                                                                      cgcAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAATGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTaa                                                    <  2:1022021/137‑1 (MQ=255)
                                                                       tCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          <  1:329547/68‑1 (MQ=255)
                                                                       tCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          >  2:329547/1‑68 (MQ=255)
                                                                               tCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                 >  2:239331/1‑85 (MQ=255)
                                                                               tCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                 <  1:239331/85‑1 (MQ=255)
                                                                                     cAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          >  2:254473/1‑54 (MQ=255)
                                                                                     cAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          <  1:254473/54‑1 (MQ=255)
                                                                                       gTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAg                                   >  1:127218/1‑139 (MQ=255)
                                                                                          aaaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                            >  2:703645/1‑111 (MQ=255)
                                                                                          aaaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                            <  1:703645/111‑1 (MQ=255)
                                                                                             gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCa                                                                                   >  1:369744/1‑85 (MQ=255)
                                                                                             gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCa                                                                                   <  2:369744/85‑1 (MQ=255)
                                                                                             gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGt                             >  1:1450785/1‑139 (MQ=255)
                                                                                             gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGt                             >  2:829226/1‑139 (MQ=255)
                                                                                             gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGt                             >  1:984179/1‑139 (MQ=255)
                                                                                               caccaTCCCTTTACTGGTGGAAACGGGACCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                      >  1:199552/1‑96 (MQ=255)
                                                                                               caccaTCCCTTTACTGGTGGAAACGGGACCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                      <  2:199552/96‑1 (MQ=255)
                                                                                                accaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                               >  1:647477/1‑54 (MQ=255)
                                                                                                accaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                               <  2:647477/54‑1 (MQ=255)
                                                                                                accaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTa                          >  1:232296/1‑139 (MQ=255)
                                                                                                 ccaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                             <  2:390212/55‑1 (MQ=255)
                                                                                                 ccaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                             >  1:390212/1‑55 (MQ=255)
                                                                                                 ccaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTaa                         <  2:385138/139‑1 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          <  2:230207/41‑1 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                          >  1:230207/1‑41 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGt                                                                                                        <  1:414941/59‑1 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGt                                                                                                        >  2:414941/1‑59 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAAcc                                                                                            >  2:250242/1‑71 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAAcc                                                                                            <  1:250242/71‑1 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                      <  1:1216551/93‑1 (MQ=255)
                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                      >  2:1216551/1‑93 (MQ=255)
                                                                                                   aTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACg                       <  1:769680/139‑1 (MQ=255)
                                                                                                           aCTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAg                                                                                <  2:1356796/74‑1 (MQ=255)
                                                                                                           aCTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAg                                                                                >  1:1356796/1‑74 (MQ=255)
                                                                                                             tggtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgc             <  1:1681337/139‑1 (MQ=255)
                                                                                                              ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                   <  2:684451/36‑1 (MQ=255)
                                                                                                              ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                   >  1:684451/1‑36 (MQ=255)
                                                                                                              ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgcc            <  2:203307/139‑1 (MQ=255)
                                                                                                              ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgcc            <  1:829226/139‑1 (MQ=255)
                                                                                                              ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgcc            <  2:1450785/139‑1 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                               <  2:1244943/39‑1 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                               >  1:1244943/1‑39 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGt                                                                                                        >  2:1385277/1‑46 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGt                                                                                                        <  1:1385277/46‑1 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                 >  2:649046/1‑53 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                 <  1:649046/53‑1 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                      <  2:85365/80‑1 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                      >  1:85365/1‑80 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCt           >  2:1675973/1‑139 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCt           >  1:265764/1‑139 (MQ=255)
                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCt           >  2:325359/1‑139 (MQ=255)
                                                                                                                tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTa                                                                                                     <  1:1374806/48‑1 (MQ=255)
                                                                                                                tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTa                                                                                                     >  2:1374806/1‑48 (MQ=255)
                                                                                                                tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                      >  1:236462/1‑79 (MQ=255)
                                                                                                                tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                      <  2:236462/79‑1 (MQ=255)
                                                                                                                tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCCATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTa          >  2:1693470/1‑139 (MQ=255)
                                                                                                                 ggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTAc         <  1:1675973/139‑1 (MQ=255)
                                                                                                                  gAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTaa                         <  2:203716/122‑1 (MQ=255)
                                                                                                                  gAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTaa                         >  1:203716/1‑122 (MQ=255)
                                                                                                                    aaCGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGttt                                                                                                      <  1:1625631/43‑1 (MQ=255)
                                                                                                                    aaCGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGttt                                                                                                      >  2:1625631/1‑43 (MQ=255)
                                                                                                                    aaCGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGtt                                                                                                       >  1:153127/1‑42 (MQ=255)
                                                                                                                    aaCGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGtt                                                                                                       <  2:153127/42‑1 (MQ=255)
                                                                                                                        gggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTa  <  2:232296/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                           
CTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCGCCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTA  >  NC_000913/2685775‑2686033

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: ‑