Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,024,049 G→C 7.0% A912P (GCG→CCG)  xdhD → putative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,024,0490GC 7.0% 246.7 / 8.0 71A912P (GCG→CCG) xdhDputative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base C (2/3);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

ATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAA  >  NC_000913/3023917‑3024180
                                                                                                                                    |                                                                                                                                   
aTCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCgg                                                                                                                      >  3:364091/1‑148 (MQ=255)
 tCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCg                                                                                                                                   <  3:29654/134‑1 (MQ=255)
    tttATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAAt                                                                                                                  >  1:84682/1‑148 (MQ=255)
       aTGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCgg                                                                                                               <  4:364091/148‑1 (MQ=255)
        tGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCg                                                                                                                >  1:197829/1‑146 (MQ=255)
        tGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCg                                                                                                                <  2:197829/146‑1 (MQ=255)
          aCGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGa                                                                                                                                  <  3:1110640/126‑1 (MQ=255)
          aCGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGa                                                                                                                                  >  4:1110640/1‑126 (MQ=255)
              cGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGgtg                                                                                                             <  4:326912/143‑1 (MQ=255)
              cGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGgtg                                                                                                             >  3:326912/1‑143 (MQ=255)
                aaGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATgg                                                                                                      >  3:994681/1‑148 (MQ=255)
                     cACGCGTTAACGGGTGATTTACGCCGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtctc                                                                                                                        >  1:1015935/1‑125 (MQ=255)
                     cACGCGTTAACGGGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtctc                                                                                                                        <  2:1015935/125‑1 (MQ=255)
                     cACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCg                                                                                                                             <  2:8274/120‑1 (MQ=255)
                     cACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCg                                                                                                                             >  1:8274/1‑120 (MQ=255)
                       cccGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCgg                                                                                                                      >  4:166550/1‑125 (MQ=255)
                       cccGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCgg                                                                                                                      <  3:166550/125‑1 (MQ=255)
                       cccGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAAt                                                                                                                  <  3:312867/129‑1 (MQ=255)
                       cccGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAAt                                                                                                                  >  4:312867/1‑129 (MQ=255)
                           ttAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCgg                                                                                           >  4:543087/1‑148 (MQ=255)
                              aCGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGa                                                                                        <  2:84682/148‑1 (MQ=255)
                               cgcgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGAt                                                                                       <  4:994681/148‑1 (MQ=255)
                                 cgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaaa                                                                                                                                >  2:620603/1‑105 (MQ=255)
                                 cgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaaa                                                                                                                                <  1:620603/105‑1 (MQ=255)
                                 cgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTg                                                                                     >  4:839423/1‑148 (MQ=255)
                                 cgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGa                                                                                                                                  <  3:80054/103‑1 (MQ=255)
                                 cgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGa                                                                                                                                  >  4:80054/1‑103 (MQ=255)
                                   tGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGc                                                                                    >  3:1120779/1‑147 (MQ=255)
                                   tGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGc                                                                                    <  4:1120779/147‑1 (MQ=255)
                                      tttACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATg                                                                                                       >  4:418478/1‑125 (MQ=255)
                                      tttACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATg                                                                                                       <  3:418478/125‑1 (MQ=255)
                                               ttACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCAc                                                                       >  4:83105/1‑148 (MQ=255)
                                                tACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTAc                                                                                 <  3:152941/137‑1 (MQ=255)
                                                tACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTAc                                                                                 >  4:152941/1‑137 (MQ=255)
                                                       cACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGatgcatg                                                               >  3:908940/1‑148 (MQ=255)
                                                            aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCgg                                                                                                               >  4:78226/1‑95 (MQ=255)
                                                            aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCgg                                                                                                               >  2:359205/1‑95 (MQ=255)
                                                            aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCgg                                                                                                               <  1:359205/95‑1 (MQ=255)
                                                            aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCgg                                                                                                               <  3:78226/95‑1 (MQ=255)
                                                            aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTAcc                                                                                >  3:22369/1‑126 (MQ=255)
                                                            aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTAcc                                                                                <  4:22369/126‑1 (MQ=255)
                                                            aaaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCa                                                          <  3:543087/148‑1 (MQ=255)
                                                                ttGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATCCATGCGGCATCTg                                                      <  3:83105/148‑1 (MQ=255)
                                                                 tGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTgg                                                     <  3:839423/148‑1 (MQ=255)
                                                                  ggTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGt                                                    <  4:908940/148‑1 (MQ=255)
                                                                         ttCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACg                                                                      >  3:351292/1‑123 (MQ=255)
                                                                         ttCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACg                                                                      <  4:351292/123‑1 (MQ=255)
                                                                               gTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtc                                                                                                                          >  3:1090877/1‑65 (MQ=255)
                                                                               gTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtc                                                                                                                          <  4:1090877/65‑1 (MQ=255)
                                                                               gTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTg                                                                                     <  1:1002227/102‑1 (MQ=255)
                                                                               gTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTg                                                                                     >  2:1002227/1‑102 (MQ=255)
                                                                                                 tGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAAt                     >  2:16565/1‑148 (MQ=255)
                                                                                                    tGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTgg                                                     <  2:187459/113‑1 (MQ=255)
                                                                                                    tGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTgg                                                     >  1:187459/1‑113 (MQ=255)
                                                                                                     gCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGatgcatg                                                               >  3:70721/1‑102 (MQ=255)
                                                                                                     gCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTCCTACCGCAATTCACGatgcatg                                                               <  4:70721/102‑1 (MQ=255)
                                                                                                         aGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTg             <  1:16565/148‑1 (MQ=255)
                                                                                                         aGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTg             >  4:590988/1‑148 (MQ=255)
                                                                                                           cgacgaTAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTgcg           <  3:590988/148‑1 (MQ=255)
                                                                                                               gaTAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATGTCACACCGGAGAAAATACTCACTgcgc          >  3:866495/1‑145 (MQ=255)
                                                                                                               gaTAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATGTCACACCGGAGAAAATACTCACTgcgc          <  4:866495/145‑1 (MQ=255)
                                                                                                                  aaaGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTg        <  2:1043576/144‑1 (MQ=255)
                                                                                                                  aaaGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTg        >  1:1043576/1‑144 (MQ=255)
                                                                                                                         gCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGaaaaa  <  2:295580/143‑1 (MQ=255)
                                                                                                                         gCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGaaaaa  >  1:295580/1‑143 (MQ=255)
                                                                                                                          cccGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGaaa                       >  4:591750/1‑121 (MQ=255)
                                                                                                                          cccGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGaaa                       <  3:591750/121‑1 (MQ=255)
                                                                                                                           ccGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACg                                                                      <  2:45254/73‑1 (MQ=255)
                                                                                                                           ccGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACg                                                                      >  1:460087/1‑73 (MQ=255)
                                                                                                                           ccGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACg                                                                      <  2:460087/73‑1 (MQ=255)
                                                                                                                           ccGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACg                                                                      >  1:45254/1‑73 (MQ=255)
                                                                                                                                  gggCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGatgcatg                                                               <  3:336165/73‑1 (MQ=255)
                                                                                                                                  gggCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGatgcatg                                                               >  4:336165/1‑73 (MQ=255)
                                                                                                                                    |                                                                                                                                   
ATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAA  >  NC_000913/3023917‑3024180

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