Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,961,917 G→C G414A (GGC→GCC)  rep → DNA helicase and single‑stranded DNA‑dependent ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,961,9170GC100.0% 267.7 / NA 73G414A (GGC→GCC) repDNA helicase and single‑stranded DNA‑dependent ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (37/36);  total (37/36)

TACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGACCCGCTT  >  NC_000913/3961781‑3962043
                                                                                                                                        |                                                                                                                              
tACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGccc                                                                                                                              >  1:241257/1‑139 (MQ=255)
        atCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACg                                                                                                                      <  2:1297624/139‑1 (MQ=255)
                gtACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGaaaaa                                                                                                              <  2:725370/139‑1 (MQ=255)
                 tACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGaaaa                                                                                                               >  1:511781/1‑137 (MQ=255)
                 tACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGaaaa                                                                                                               <  2:511781/137‑1 (MQ=255)
                 tACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAg                                                                                                             >  1:100780/1‑139 (MQ=255)
                 tACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAg                                                                                                             >  1:910190/1‑139 (MQ=255)
                   cgtcgtTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAg                                                                                                             >  2:196400/1‑137 (MQ=255)
                   cgtcgtTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAg                                                                                                             <  1:196400/137‑1 (MQ=255)
                   cgtcgtTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCt                                                                                                           >  1:864238/1‑139 (MQ=255)
                     tcgtTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTgg                                                                                                         <  2:241257/139‑1 (MQ=255)
                       gtTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGt                                                                                                       >  2:1266401/1‑139 (MQ=255)
                       gtTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGt                                                                                                       >  1:530740/1‑139 (MQ=255)
                                   tCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATg                                                                                           >  2:1159097/1‑139 (MQ=255)
                                      tGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg                                                                                        <  2:910190/139‑1 (MQ=255)
                                      tGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg                                                                                        <  1:294466/139‑1 (MQ=255)
                                        aaaTCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcgcg                                                                                      <  1:1107385/139‑1 (MQ=255)
                                           tCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCaa                                                                                   <  1:1266401/139‑1 (MQ=255)
                                              aGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATaa                                                                                >  1:137211/1‑139 (MQ=255)
                                                   ttGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCa                                                                           <  2:530740/139‑1 (MQ=255)
                                                           ttATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcgcg                                                                                      <  2:990845/120‑1 (MQ=255)
                                                           ttATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcgcg                                                                                      >  1:990845/1‑120 (MQ=255)
                                                             aTCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGc                                                                 <  1:1250545/139‑1 (MQ=255)
                                                              tCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGcc                                                                >  2:615176/1‑139 (MQ=255)
                                                               cTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCa                                                               <  1:1221829/139‑1 (MQ=255)
                                                                   gcgTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGaaaa                                                                                                               >  2:966527/1‑87 (MQ=255)
                                                                   gcgTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGaaaa                                                                                                               <  1:966527/87‑1 (MQ=255)
                                                                   gcgTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCtt                                                           >  1:187629/1‑139 (MQ=255)
                                                                   gcgTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCtt                                                           >  1:1298959/1‑139 (MQ=255)
                                                                        cTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCt                                                            >  1:76012/1‑133 (MQ=255)
                                                                        cTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCt                                                            <  2:76012/133‑1 (MQ=255)
                                                                               aCCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGc                                                                 <  2:952099/121‑1 (MQ=255)
                                                                               aCCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGc                                                                 >  1:952099/1‑121 (MQ=255)
                                                                                cccGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCa                                                                           <  1:187733/110‑1 (MQ=255)
                                                                                cccGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGAGCGATGACGCGCAATAAAAGCa                                                                           >  2:187733/1‑110 (MQ=255)
                                                                                     aCGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCa                                         >  1:712453/1‑139 (MQ=255)
                                                                                          gACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGc                                    >  1:252490/1‑139 (MQ=255)
                                                                                               cgcATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGaaaa                                                                                                               >  2:437655/1‑59 (MQ=255)
                                                                                               cgcATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGaaaa                                                                                                               <  1:437655/59‑1 (MQ=255)
                                                                                               cgcATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGc                               >  2:511797/1‑139 (MQ=255)
                                                                                               cgcATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGc                               >  1:554333/1‑139 (MQ=255)
                                                                                                  aTTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGa                            <  2:1298959/139‑1 (MQ=255)
                                                                                                  aTTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGa                            <  2:1304491/139‑1 (MQ=255)
                                                                                                   tttCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTg                                                                                                                   <  1:911608/51‑1 (MQ=255)
                                                                                                   tttCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTg                                                                                                                   >  2:911608/1‑51 (MQ=255)
                                                                                                   tttCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGAc                           <  2:712453/139‑1 (MQ=255)
                                                                                                    ttCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACg                          <  2:554333/139‑1 (MQ=255)
                                                                                                    ttCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACg                          <  1:511797/139‑1 (MQ=255)
                                                                                                          gTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTa                    >  2:1164171/1‑139 (MQ=255)
                                                                                                              cGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGc                                                                 <  1:1006679/90‑1 (MQ=255)
                                                                                                              cGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGc                                                                 >  2:1006679/1‑90 (MQ=255)
                                                                                                              cGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGaa                >  2:738461/1‑139 (MQ=255)
                                                                                                              cGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACAGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGaa                >  2:877352/1‑139 (MQ=255)
                                                                                                                ttAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGc              >  1:1142909/1‑139 (MQ=255)
                                                                                                                  aaCACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCa                                                                                    >  1:411385/1‑67 (MQ=255)
                                                                                                                  aaCACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCa                                                                                    <  2:411385/67‑1 (MQ=255)
                                                                                                                   acacGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg                                                                                        >  1:761859/1‑62 (MQ=255)
                                                                                                                   acacGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg                                                                                        <  2:761859/62‑1 (MQ=255)
                                                                                                                     acGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcgcg                                                                                      >  1:1265892/1‑62 (MQ=255)
                                                                                                                     acGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcgcg                                                                                      <  2:1265892/62‑1 (MQ=255)
                                                                                                                      cGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGAc        >  2:917202/1‑139 (MQ=255)
                                                                                                                       gCCGAAGCAAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGAcc       <  2:998052/139‑1 (MQ=255)
                                                                                                                        ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg                                                                                        >  2:342812/1‑57 (MQ=255)
                                                                                                                        ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg                                                                                        <  1:342812/57‑1 (MQ=255)
                                                                                                                        ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGc                                                                 >  1:523281/1‑80 (MQ=255)
                                                                                                                        ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGc                                                                 <  2:523281/80‑1 (MQ=255)
                                                                                                                        ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCtt                                                           >  1:315297/1‑86 (MQ=255)
                                                                                                                        ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCtt                                                           <  2:315297/86‑1 (MQ=255)
                                                                                                                        ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGAccc      <  1:738461/139‑1 (MQ=255)
                                                                                                                        ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCACAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGAccc      >  2:127725/1‑139 (MQ=255)
                                                                                                                           aaGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGACCCGCt   <  1:1159097/139‑1 (MQ=255)
                                                                                                                            aGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCACAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGACCCGCtt  <  1:127725/139‑1 (MQ=255)
                                                                                                                                  agaTTGCCCCGGCTACGCTGCAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCa                                         <  1:807858/94‑1 (MQ=255)
                                                                                                                                  agaTTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCa                                         >  2:807858/1‑94 (MQ=255)
                                                                                                                                        |                                                                                                                              
TACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGACCCGCTT  >  NC_000913/3961781‑3962043

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