Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,569,309 G→A pseudogene (312/921 nt) yjiP → pseudogene, transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,569,3090GA98.7% 235.9 / NA 74pseudogene (312/921 nt)yjiPpseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/38);  minor base C (1/0);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAC  >  NC_000913/4569175‑4569439
                                                                                                                                      |                                                                                                                                  
cTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGcc                                                                                                                                <  1:968962/139‑1 (MQ=255)
cTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGcc                                                                                                                                <  2:387347/139‑1 (MQ=255)
        cGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATagag                                                                                                                        <  1:1184804/139‑1 (MQ=255)
          ctatTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATCTCTATGTGGTGATTGAACATCAGAGCCGCGAGGACGTTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATATACCGCCATATAGAGCa                                                                                                                      >  1:1041339/2‑139 (MQ=255)
           taatttATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCAt                                                                                                                     >  1:645300/4‑139 (MQ=255)
             ttttATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGa                                                                                                                   >  2:1369349/1‑139 (MQ=255)
                  tGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAAc                                                                                                              >  1:1692668/1‑139 (MQ=255)
                   ggTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACg                                                                                                             <  2:645300/139‑1 (MQ=255)
                     tCGGTAAAGACCCGCGAAGGCGATGGGTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGACTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCAtata                                                                                                                           <  1:1237702/123‑1 (MQ=255)
                     tCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCAtata                                                                                                                           >  2:1237702/1‑123 (MQ=255)
                       ggTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAg                                                                                                         >  1:1343499/1‑139 (MQ=255)
                        gTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAacg                                                                                                                                  >  2:564867/1‑113 (MQ=255)
                          aaaGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                      <  2:844089/139‑1 (MQ=255)
                          aaaGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                      <  2:725543/139‑1 (MQ=255)
                          aaaGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                      <  1:1369349/139‑1 (MQ=255)
                           aaGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGc                                                                                                     <  1:805629/139‑1 (MQ=255)
                           aaGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGc                                                                                                     <  2:1628599/139‑1 (MQ=255)
                           aaGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGc                                                                                                     <  2:549964/139‑1 (MQ=255)
                                   cgAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGc                                                                                                     <  1:868138/131‑1 (MQ=255)
                                   cgAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGc                                                                                                     >  2:868138/1‑131 (MQ=255)
                                   cgAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTg                                                                                             >  1:448226/1‑139 (MQ=255)
                                     aaGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGt                                                                                           <  2:1692668/139‑1 (MQ=255)
                                     aaGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGt                                                                                           <  2:1343499/139‑1 (MQ=255)
                                              gCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGAt                                                                                  <  2:448226/139‑1 (MQ=255)
                                                  tatCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCata                                                                                                                             >  1:645460/1‑92 (MQ=255)
                                                  tatCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCata                                                                                                                             <  2:645460/92‑1 (MQ=255)
                                                        tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCata                                                                                                                             >  2:1254018/1‑86 (MQ=255)
                                                        tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCata                                                                                                                             <  1:1254018/86‑1 (MQ=255)
                                                        tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGcc                                                                                                           >  1:50334/1‑104 (MQ=255)
                                                        tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGcc                                                                                                           <  2:50334/104‑1 (MQ=255)
                                                        tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                        >  2:953152/1‑139 (MQ=255)
                                                        tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                        <  2:1551187/139‑1 (MQ=255)
                                                        tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                        >  2:1466124/1‑139 (MQ=255)
                                                        tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                        >  1:1088574/1‑139 (MQ=255)
                                                            gtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                        >  1:983122/1‑135 (MQ=255)
                                                            gtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                        <  2:983122/135‑1 (MQ=255)
                                                            gtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACg                                                                    >  2:921485/1‑139 (MQ=255)
                                                              gATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCg                                                                                                <  2:825821/109‑1 (MQ=255)
                                                              gATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCg                                                                                                >  1:825821/1‑109 (MQ=255)
                                                               aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAacg                                                                                                                                  >  2:230991/1‑74 (MQ=255)
                                                               aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAacg                                                                                                                                  <  1:230991/74‑1 (MQ=255)
                                                               aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                      >  1:1557878/1‑118 (MQ=255)
                                                               aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                      <  2:1557878/118‑1 (MQ=255)
                                                               aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTa                                                                 >  2:1134597/1‑139 (MQ=255)
                                                               aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTa                                                                 >  2:564870/1‑139 (MQ=255)
                                                                  gAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACg                                                                                                             >  2:1143815/1‑92 (MQ=255)
                                                                  gAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCATGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACg                                                                                                             <  1:1143815/92‑1 (MQ=255)
                                                                       tCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGc                                                                                                                       <  1:157453/77‑1 (MQ=255)
                                                                       tCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGc                                                                                                                       >  2:157453/1‑77 (MQ=255)
                                                                          gagCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                      <  1:1466124/139‑1 (MQ=255)
                                                                           agCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCg                                                                                    >  2:185392/1‑108 (MQ=255)
                                                                           agCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCg                                                                                    <  1:185392/108‑1 (MQ=255)
                                                                                cgAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTAcccc                                                <  1:921485/139‑1 (MQ=255)
                                                                                  aGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGGTGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                      >  2:1001021/1‑83 (MQ=255)
                                                                                  aGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                      <  1:1001021/83‑1 (MQ=255)
                                                                                             atGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGCTAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTgct                                   >  2:359036/1‑139 (MQ=255)
                                                                                               ggCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                      >  2:1372340/1‑118 (MQ=255)
                                                                                               ggCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                      <  1:1372340/118‑1 (MQ=255)
                                                                                                     tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTAcc                                                                                                 >  2:367116/1‑69 (MQ=255)
                                                                                                     tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTAcc                                                                                                 <  1:367116/69‑1 (MQ=255)
                                                                                                     tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGAc                           >  2:548441/1‑139 (MQ=255)
                                                                                                             tGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGc                   >  2:1136684/1‑139 (MQ=255)
                                                                                                                gcTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGccga                <  1:548441/139‑1 (MQ=255)
                                                                                                                 cTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCAc                                                                     <  2:123175/85‑1 (MQ=255)
                                                                                                                 cTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCAc                                                                     >  1:123175/1‑85 (MQ=255)
                                                                                                                      ccATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGccgacccgac          <  2:1088574/139‑1 (MQ=255)
                                                                                                                      ccATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGccgacccgac          <  1:1134597/139‑1 (MQ=255)
                                                                                                                      ccATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGccgacccgac          <  1:564870/139‑1 (MQ=255)
                                                                                                                         tGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGt                                                                  >  2:620553/1‑80 (MQ=255)
                                                                                                                         tGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGt                                                                  <  1:620553/80‑1 (MQ=255)
                                                                                                                              gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGt                                                                                           <  2:471486/50‑1 (MQ=255)
                                                                                                                              gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGt                                                                                           >  1:471486/1‑50 (MQ=255)
                                                                                                                              gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCg                                                             >  1:905930/1‑80 (MQ=255)
                                                                                                                              gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCg                                                             <  2:905930/80‑1 (MQ=255)
                                                                                                                              gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc  >  1:694044/1‑139 (MQ=255)
                                                                                                                              gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc  >  1:674271/1‑139 (MQ=255)
                                                                                                                              gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc  >  1:1067185/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAC  >  NC_000913/4569175‑4569439

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