Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 758,983 T→A L93Q (CTG→CAG)  sucA → 2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913758,9830TA99.4% 543.6 / NA 154L93Q (CTG→CAG) sucA2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (74/79);  total (74/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

TTCCAGCAGTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGT  >  NC_000913/758838‑759127
                                                                                                                                                 |                                                                                                                                                
ttCCAGCAGTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCa                                                                                                                                               <  5:381223/149‑1 (MQ=255)
   cagcagTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCagcagc                                                                                                                                             <  7:99222/148‑1 (MQ=255)
   cagcagTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCagcagc                                                                                                                                             <  6:492648/148‑1 (MQ=255)
   cagcagTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGGGAAGCAGGTTAAAGTCCagcagc                                                                                                                                             <  3:359531/148‑1 (MQ=255)
      gagTTACCTGGTACGGGAGTGAAACCGGATCAATTCGACTCTCTATGGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGTCCCTGACACCAATGTGAAGCAGGTTAAAGTCGTGCAGCTCAt                                                                                                                                         <  2:215756/148‑1 (MQ=255)
        gTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCAtt                                                                                                                                        >  4:321425/1‑148 (MQ=255)
        gTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCAGACCCTGCCACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCAtt                                                                                                                                        >  5:27541/1‑148 (MQ=255)
        gTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCAGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCAtt                                                                                                                                        >  5:175319/1‑148 (MQ=255)
           aCCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAGGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAAc                                                                                                                                     >  5:701662/1‑148 (MQ=255)
           aCCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCGCATTaag                                                                                                                                     <  5:511350/148‑2 (MQ=255)
            ccTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACg                                                                                                                                    >  6:489921/1‑148 (MQ=255)
             cTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCa                                                                                                                                  >  4:639742/1‑149 (MQ=255)
             cTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCa                                                                                                                                  >  3:448169/1‑149 (MQ=255)
              tGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCa                                                                                                                                  >  4:21670/1‑148 (MQ=255)
               ggTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATa                                                                                                                                >  7:544297/1‑149 (MQ=255)
               ggTACGGGAGTCAAAACGCATCAAATCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAAAGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATa                                                                                                                                >  6:128520/1‑149 (MQ=255)
                gTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCTATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATAc                                                                                                                               <  5:662313/149‑1 (MQ=255)
                     ggAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCtt                                                                                                                          <  7:503279/149‑1 (MQ=255)
                        gTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCg                                                                                                                       <  1:305680/149‑1 (MQ=255)
                        gTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCg                                                                                                                       >  3:563023/1‑149 (MQ=255)
                           ataCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTgg                                                                                                                    <  7:136841/147‑1 (MQ=255)
                           aaaCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTgg                                                                                                                    >  4:461253/1‑149 (MQ=255)
                           aaaCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACACTTCAACGATCTCCGACCCAGACACCAATGTGAAGCACGATAAAGTCCAGCAGCTCATAAACGCATACCGCTTCCGTg                                                                                                                     >  8:688244/1‑148 (MQ=255)
                              ccGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCACCTCATTAACGCATACCGCTTCCGTGGTCa                                                                                                                 >  3:578821/1‑149 (MQ=255)
                               cGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTc                                                                                                                                           <  2:423654/122‑1 (MQ=255)
                                 gATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCa                                                                                                              <  5:489921/149‑1 (MQ=255)
                                      aTTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATg                                                                                                         >  8:624581/1‑149 (MQ=255)
                                       ttCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTAAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGc                                                                                                        <  3:597753/149‑1 (MQ=255)
                                           aCTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCa                                                                                                              >  7:595201/1‑139 (MQ=255)
                                           aCTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAt                                                                                                    >  2:410017/1‑149 (MQ=255)
                                            ctctCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAtc                                                                                                   >  6:567667/1‑149 (MQ=255)
                                             tctcAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAtct                                                                                                  >  7:313667/1‑149 (MQ=255)
                                               tcAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAtctc                                                                                                 >  1:496212/1‑148 (MQ=255)
                                                 aaaCGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGa                                                                                               <  8:708944/148‑1 (MQ=255)
                                                 aaaCGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGAt                                                                                              <  7:82725/149‑1 (MQ=255)
                                                 aaaCGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGCGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGc                                                                                                        <  7:437077/139‑1 (MQ=255)
                                                  aaCGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATc                                                                                             <  5:600146/149‑1 (MQ=255)
                                                    cgcgTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCg                                                                                           <  8:544297/149‑1 (MQ=255)
                                                     gcgTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTcc                                                                                                                        >  8:193798/1‑119 (MQ=255)
                                                     gcgTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTcc                                                                                                                        <  7:193798/119‑1 (MQ=255)
                                                      cgTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCt                                                                                         >  2:285682/1‑149 (MQ=255)
                                                      cgTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCt                                                                                         <  5:311444/149‑1 (MQ=255)
                                                       gTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTGCGACCCTGACACCAATGTGAAGCAGGTTAAAGGCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTg                                                                                        <  2:522425/149‑1 (MQ=255)
                                                         gAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTGCCGTGGTCACCAGCATGCGAATCTCGCTCCGCTggg                                                                                      >  7:279520/1‑149 (MQ=255)
                                                          aaTATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGa                                                                                     <  4:448169/149‑1 (MQ=255)
                                                          aaTATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCACCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGa                                                                                     <  4:578821/149‑1 (MQ=255)
                                                              tttCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACtgt                                                                                 <  1:285682/149‑1 (MQ=255)
                                                              tttCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACtgt                                                                                 <  3:18762/149‑1 (MQ=255)
                                                              tttCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGGCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACtgt                                                                                 <  8:255825/149‑1 (MQ=255)
                                                                 ccgccgCCTGGCGAAAGACTCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTgg                                                                               <  7:355515/148‑1 (MQ=255)
                                                                 ccgccgCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTgg                                                                               <  6:107540/148‑1 (MQ=255)
                                                                  cgccgcCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGgca                                                                             <  1:410017/149‑1 (MQ=255)
                                                                  cgccgcCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGgc                                                                              <  3:639742/148‑1 (MQ=255)
                                                                       ccTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGgcagc                                                                           >  1:144804/1‑146 (MQ=255)
                                                                          ggCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCAc                                                                                                                >  1:149678/1‑106 (MQ=255)
                                                                          ggCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCAc                                                                                                                <  2:149678/106‑1 (MQ=255)
                                                                           gCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATa                                                                     >  7:505176/1‑148 (MQ=255)
                                                                             gAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAAc                                                                                                                                     >  3:633410/1‑82 (MQ=255)
                                                                              aaaGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGt                                                                 <  6:310491/149‑1 (MQ=255)
                                                                                  aCGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATa                                                                     >  8:580262/1‑141 (MQ=255)
                                                                                      ttCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATc                                                         >  2:203709/1‑149 (MQ=255)
                                                                                        cACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCt                                                        <  5:567667/148‑1 (MQ=255)
                                                                                         aCGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTg                                                       >  2:525762/1‑148 (MQ=255)
                                                                                          cGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGa                                                     <  4:563023/149‑1 (MQ=255)
                                                                                           gTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGc                                                                                                        <  7:5566/97‑1 (MQ=255)
                                                                                           gTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGc                                                                                                        >  8:5566/1‑97 (MQ=255)
                                                                                           gTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGAGGCCGATCTGGAt                                                    <  8:505176/149‑1 (MQ=255)
                                                                                              aCTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTg                                                                                                                     <  8:427449/81‑1 (MQ=255)
                                                                                              aCTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTg                                                                                                                     >  7:427449/1‑81 (MQ=255)
                                                                                                tctTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTc                                               >  7:715365/1‑149 (MQ=255)
                                                                                                  ttCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGAt                                                    >  8:507937/1‑142 (MQ=255)
                                                                                                  ttCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCtt                                             >  6:261697/1‑149 (MQ=255)
                                                                                                  ttCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCtt                                             >  6:174447/1‑149 (MQ=255)
                                                                                                   tCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGCTCCTCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCttt                                            >  6:646929/1‑149 (MQ=255)
                                                                                                    cAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCt                                              <  7:624581/146‑1 (MQ=255)
                                                                                                     aaCGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATc                                                         >  5:381861/1‑134 (MQ=255)
                                                                                                        gATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCAc                                        >  7:753716/1‑148 (MQ=255)
                                                                                                             ccGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCtt                                                                                                                          <  5:640546/61‑1 (MQ=255)
                                                                                                              cGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCt                                   <  7:415491/147‑1 (MQ=255)
                                                                                                              cGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGa                                 >  6:302560/1‑149 (MQ=255)
                                                                                                               gACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCt                                                        <  3:175682/125‑1 (MQ=255)
                                                                                                               gACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCt                                                        >  4:175682/1‑125 (MQ=255)
                                                                                                                   cTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGTAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCg                              >  8:363933/1‑147 (MQ=255)
                                                                                                                    tGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGTACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAg                           >  6:509761/1‑149 (MQ=255)
                                                                                                                        acCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGa                       >  6:722539/1‑149 (MQ=255)
                                                                                                                          cAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCg                                                 <  4:82999/121‑1 (MQ=255)
                                                                                                                          cAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCg                                                 >  3:82999/1‑121 (MQ=255)
                                                                                                                          cAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACt                     >  3:374530/1‑149 (MQ=255)
                                                                                                                          cAATGTGAAGCAGGTTAAAGTACAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACt                     >  4:344461/1‑149 (MQ=255)
                                                                                                                            aTGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTc                   <  2:354315/149‑1 (MQ=255)
                                                                                                                              gtgAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTAACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTc                   >  2:726905/1‑147 (MQ=255)
                                                                                                                               tgAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCa                 <  8:753716/148‑1 (MQ=255)
                                                                                                                               tgAAGCAGGGTAAAGTCCAGCAGCTCATTAACGCATACCGATTCCGTGGTCACCAGCATTCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAg                <  8:279520/149‑1 (MQ=255)
                                                                                                                                 aaGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAg                >  6:17107/1‑147 (MQ=255)
                                                                                                                                      ggTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAAATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCt         <  3:564728/149‑1 (MQ=255)
                                                                                                                                       gTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAg                                                                  >  7:628626/1‑91 (MQ=255)
                                                                                                                                       gTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAg                                                                  <  8:628626/91‑1 (MQ=255)
                                                                                                                                           aaGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAAc    >  3:261733/1‑149 (MQ=255)
                                                                                                                                           aaGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAAc    >  3:251478/1‑149 (MQ=255)
                                                                                                                                             gTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGt  >  4:624324/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCCAGCAGTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGT  >  NC_000913/758838‑759127

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