Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,024,049 G→C 5.6% 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 5.6% 242.0 / 5.7 72A912P (GCG→CCG) xdhDputative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding
Reads supporting (aligned to +/- strand):  ref base G (34/34);  new base C (2/2);  total (36/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 = 7.59e-01

CAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGG  >  NC_000913/3023901‑3024197
                                                                                                                                                    |                                                                                                                                                    
caGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCgggg                                                                                                                                                      <  1:242275/149‑1 (MQ=255)
   tATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGCCATTCCGCGTGATTTCCGCGCTGTCCTGGTGCCGAGCGACGATAAAGTCGGCCCGGTCgcggcg                                                                                                                                                    >  2:466233/1‑148 (MQ=255)
     tGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGaaa                                                                                                                                                 <  6:104145/149‑1 (MQ=255)
             gagaTCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAtctc                                                                                                                                         >  8:325970/1‑149 (MQ=255)
              agaTCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCg                                                                                                                                        <  7:86082/149‑1 (MQ=255)
              agaTCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCg                                                                                                                                        <  8:141631/149‑1 (MQ=255)
                    tttATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATc                                                                                                                                  <  5:45627/149‑1 (MQ=255)
                         gACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGgtg                                                                                                                              >  1:199699/1‑148 (MQ=255)
                          aCGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGgtgt                                                                                                                             >  8:104730/1‑148 (MQ=255)
                              cGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATg                                                                                                                        >  4:427736/1‑149 (MQ=255)
                              cGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATg                                                                                                                        <  2:199699/149‑1 (MQ=255)
                              cGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATg                                                                                                                        >  2:249995/1‑149 (MQ=255)
                               gAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAGCTCGGAAATCGGTGTAAATgg                                                                                                                       >  5:50896/1‑149 (MQ=255)
                                aaGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGgc                                                                                                                      <  7:104730/149‑1 (MQ=255)
                                  ggTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGgcg                                                                                                                     >  3:148084/1‑148 (MQ=255)
                                  ggTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGAGCTCGGAAATCGGTGTAAATGgcgc                                                                                                                    <  6:50896/149‑1 (MQ=255)
                                   gTCACCCGTTAACGCGTGATTTACGCAGTTACGGAGCACCGAAAATTTGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGGGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGAGAAATCGATCTCGGAAATCGGTGTAAATGGCGCa                                                                                                                   <  8:382732/149‑1 (MQ=255)
                                      aCCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCt                                                                                                                >  6:371280/1‑149 (MQ=255)
                                         agTTTACGCTTGATTTACGCTGTTACGGGGCTCCGAAAATTGGTGACATTCCGCGTTATTTGCGCGCTGTGCTGGTGCCGAGAGACTTTAAAGTCGGCCAGTTCGGGGCGAAATAGATCTCGGAAATCGGTGTAAATTGCGCAGCTCCg                                                                                                             <  8:253797/148‑1 (MQ=255)
                                           ttAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGc                                                                                                           <  7:325970/149‑1 (MQ=255)
                                           ttAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGc                                                                                                           <  4:16911/149‑1 (MQ=255)
                                                gcgTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTg                                                                                                      >  7:12191/1‑149 (MQ=255)
                                                     aTTTACGCAGTTACGGCGCACCGAAATTTGGTGACATTCCGCGGGATTTCCGCGCTGTGCTGGGGCCGAGCGACGATAACGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTAc                                                                                                  <  1:466233/148‑1 (MQ=255)
                                                     aTTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAACTGGCGCAGCTCCGGCGATTGCTAcc                                                                                                 >  5:276895/1‑149 (MQ=255)
                                                      tttACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCg                                                                                                >  2:107998/1‑149 (MQ=255)
                                                           gCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAAtt                                                                                           >  5:161787/1‑149 (MQ=255)
                                                           gCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAAtt                                                                                           >  4:308569/1‑149 (MQ=255)
                                                            cAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCAGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGact                                                                                                                                                 <  4:420284/94‑3 (MQ=255)
                                                            cAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCAGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGact                                                                                                                                                 <  4:413656/94‑3 (MQ=255)
                                                            cAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCAGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGact                                                                                                                                                 >  3:420284/1‑92 (MQ=255)
                                                            cAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCAGGTGCCGAGCGACGATAAAGTCGCCCCGAACGGGCCGact                                                                                                                                                 >  3:413656/1‑92 (MQ=255)
                                                              gTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGGGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCAc                                                                                        <  5:219523/149‑1 (MQ=255)
                                                               ttACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCAc                                                                                        >  7:177139/1‑148 (MQ=255)
                                                                        aCCGAAAATTGGTGACATTCCACGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCg                                                                              <  8:311950/149‑1 (MQ=255)
                                                                          cGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGAtt                                                                                                       >  8:326069/1‑122 (MQ=255)
                                                                          cGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGAtt                                                                                                       <  7:326069/122‑1 (MQ=255)
                                                                             aaaTTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATc                                                                         >  8:259401/1‑149 (MQ=255)
                                                                                     gACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTAcg                                                                 >  8:256582/1‑149 (MQ=255)
                                                                                      aCATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTAcgc                                                                >  6:150386/1‑149 (MQ=255)
                                                                                      aCATTCCGCGTGATTTCCGCGCTGTGCTAGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTAcgc                                                                >  8:340234/1‑149 (MQ=255)
                                                                                      aCATTCCGCGTGATTTCCGCGCTGTGCTAGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTAcgc                                                                <  7:340234/149‑1 (MQ=255)
                                                                                         ttCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGaa                                                             <  8:177139/149‑1 (MQ=255)
                                                                                           ccGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAAt                                                            <  3:308569/148‑1 (MQ=255)
                                                                                           ccGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGGCGGCCCGGTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCGCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAAt                                                            <  2:339624/148‑1 (MQ=255)
                                                                                             gcgTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATgg                                                          >  3:157154/1‑148 (MQ=255)
                                                                                             gcgTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATgg                                                          >  4:125758/1‑148 (MQ=255)
                                                                                               gTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCAt                                                       <  5:371280/149‑1 (MQ=255)
                                                                                                   tttCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTca                                                   >  1:445866/1‑149 (MQ=255)
                                                                                                   tttCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTca                                                   <  6:161787/149‑1 (MQ=255)
                                                                                                   tttCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGAGGCATGCGGCATCTGGTTACGCGAATGGCATTTca                                                   <  6:276895/149‑1 (MQ=255)
                                                                                                    ttCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGAtt                                                                                                       >  4:404131/1‑96 (MQ=255)
                                                                                                    ttCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGAtt                                                                                                       <  3:404131/96‑1 (MQ=255)
                                                                                                                     gCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATc                                                                         <  1:304005/109‑1 (MQ=255)
                                                                                                                     gCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATc                                                                         >  2:304005/1‑109 (MQ=255)
                                                                                                                     gCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCa                                 <  1:249995/149‑1 (MQ=255)
                                                                                                                           cgacgaTAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTgcgc                           >  5:431158/1‑149 (MQ=255)
                                                                                                                            gacgaTAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCt                          >  1:188451/1‑149 (MQ=255)
                                                                                                                             acgaTAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCAtt                                                      <  4:423408/120‑1 (MQ=255)
                                                                                                                             acgaTAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCAtt                                                      >  3:423408/1‑120 (MQ=255)
                                                                                                                                 tAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCGATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGaa                      >  2:171357/1‑148 (MQ=255)
                                                                                                                                  aaaGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGaaa                     >  8:142790/1‑148 (MQ=255)
                                                                                                                                   aaGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTc                                                    >  3:404888/1‑116 (MQ=255)
                                                                                                                                   aaGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTc                                                    <  4:404888/116‑1 (MQ=255)
                                                                                                                                     gTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGaaaaaa                  >  1:317946/1‑148 (MQ=255)
                                                                                                                                        ggCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAatat               <  1:462851/148‑1 (MQ=255)
                                                                                                                                         gCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATataa             <  5:379714/149‑1 (MQ=255)
                                                                                                                                           ccGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATgg                                                          <  6:240181/102‑1 (MQ=255)
                                                                                                                                           ccGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATgg                                                          >  5:240181/1‑102 (MQ=255)
                                                                                                                                           ccGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATataaa            >  6:282190/1‑148 (MQ=255)
                                                                                                                                            cGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAatat               <  4:18859/144‑1 (MQ=255)
                                                                                                                                            cGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATataagat          >  3:18859/1‑146 (MQ=255)
                                                                                                                                             gTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATataaataa         <  3:427736/149‑1 (MQ=255)
                                                                                                                                              ttCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAaataat        <  7:142790/149‑1 (MQ=255)
                                                                                                                                              ttCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAaataat        <  8:12191/149‑1 (MQ=255)
                                                                                                                                                  gggCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAata                >  8:216078/1‑137 (MQ=255)
                                                                                                                                                  gggCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAata                <  7:216078/137‑1 (MQ=255)
                                                                                                                                                  gggCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCAGCTGGTTACGCGAATGGCATTTCACCCCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAAt    >  8:110691/1‑149 (MQ=255)
                                                                                                                                                    gCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCGATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATgg  <  1:171357/149‑1 (MQ=255)
                                                                                                                                                    gCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATataa             <  5:23251/138‑1 (MQ=255)
                                                                                                                                                    gCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATataa             >  6:23251/1‑138 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGG  >  NC_000913/3023901‑3024197

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