Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 758,983 T→A 67.2% 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,9830TA67.2% 88.3 / 81.0 89L93Q (CTG→CAG) sucA2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (14/15);  new base A (33/27);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

TCCAGCAGTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAAC  >  NC_000913/758839‑759125
                                                                                                                                                |                                                                                                                                              
tCCAGCAGTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTATAAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAAAGAGCAGCAGGTTAAAGTCAAGCAg                                                                                                                                            >  8:186976/1‑149 (MQ=255)
  cagcagTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCt                                                                                                                                                >  4:111621/1‑143 (MQ=255)
  cagcagTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCt                                                                                                                                                <  3:111621/143‑1 (MQ=255)
  cagcagTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCt                                                                                                                                          <  8:144371/149‑1 (MQ=255)
  cagcagTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTACGACCCTGACACCAATGTGAAGCAGGTTAAAGTCAAGCAGCt                                                                                                                                          >  7:77476/1‑149 (MQ=255)
    gcagTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCa                                                                                                                                        >  6:351103/1‑149 (MQ=255)
    gcagTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCa                                                                                                                                        <  4:260384/149‑1 (MQ=255)
          aCCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACg                                                                                                                                  >  2:355241/1‑149 (MQ=255)
           ccTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACg                                                                                                                                  >  4:108711/1‑148 (MQ=255)
            cTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCa                                                                                                                                >  3:82604/1‑149 (MQ=255)
            cTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGAGCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCa                                                                                                                                <  6:223832/149‑1 (MQ=255)
               gTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATa                                                                                                                              <  2:27653/148‑1 (MQ=255)
               gTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATAc                                                                                                                             >  4:326878/1‑149 (MQ=255)
                  cGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCg                                                                                                                           <  1:355241/148‑1 (MQ=255)
                    ggAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGCCCCTGACACCAATGTGAGGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCtt                                                                                                                        >  2:193567/1‑149 (MQ=255)
                     gAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTc                                                                                                                       >  4:85151/1‑149 (MQ=255)
                        tCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCACCGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGt                                                                                                                    >  1:334080/1‑149 (MQ=255)
                        tCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGt                                                                                                                    >  1:2525/1‑149 (MQ=255)
                          aaaCCGGATCAATTCCACTCTCAAAGGCGGGAATATTTCCGCCGCCTGGCGAAAGACGATTCACGTTACGCTTCAACGATCTCCGACCCTGACACCAATGGGAAGCAGGTTACCGTCCAGCAGCTCATTCACGCATACCGCTTCCGTg                                                                                                                   <  1:275158/148‑1 (MQ=255)
                          aaaCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGCCCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTgg                                                                                                                  >  2:295093/1‑149 (MQ=255)
                          aaaCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTgg                                                                                                                  >  8:347374/1‑149 (MQ=255)
                          aaaCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTg                                                                                                                   <  2:167375/148‑1 (MQ=255)
                          aaaCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTg                                                                                                                   <  3:334244/148‑1 (MQ=255)
                           aaCCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCAGACTTGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTTAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATAACGCTTCCGTGGt                                                                                                                 <  7:236970/149‑1 (MQ=255)
                               ggATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGGCAcc                                                                                                             <  8:108962/149‑1 (MQ=255)
                                gATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCAc                                                                                                              >  8:118832/1‑147 (MQ=255)
                                gATCAATTCCACTCTCAAACGCGGGAATATTTCCGCCGCCTGGCGAAAGACGCGTCACGTTAATCTTCCACGATCTCCGCCCCTGACCCCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTACCGCATACCGCTTCCGTGGTCAc                                                                                                              >  6:328886/1‑147 (MQ=255)
                                     aTTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATg                                                                                                       >  3:233332/1‑149 (MQ=255)
                                     aTTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATg                                                                                                       >  3:63297/1‑149 (MQ=255)
                                     aTTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATg                                                                                                       >  3:356524/1‑149 (MQ=255)
                                     aTTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATg                                                                                                       >  1:199433/1‑149 (MQ=255)
                                         cACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGaa                                                                                                   >  3:138617/1‑149 (MQ=255)
                                         cACTCTCAAACGCGTGAATATATCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGa                                                                                                    <  8:145248/148‑1 (MQ=255)
                                           ctctCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAtc                                                                                                 >  4:227888/1‑149 (MQ=255)
                                           ctctCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAtc                                                                                                 <  7:118832/149‑1 (MQ=255)
                                            tctcAAACGCGTGAATATTTCAGCCGCCAGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTACGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAtct                                                                                                >  3:221587/1‑149 (MQ=255)
                                              tcAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAtctc                                                                                               >  1:115257/1‑148 (MQ=255)
                                              tcAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGAGCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAAtct                                                                                                <  4:127814/147‑1 (MQ=255)
                                                aaaCGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCGTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAGTCTCGAt                                                                                            >  8:161875/1‑149 (MQ=255)
                                                aaaCGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCATGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCg                                                                                              <  8:293106/147‑1 (MQ=255)
                                                aaaCGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGAt                                                                                            <  5:314910/149‑1 (MQ=255)
                                                  aCGCGTGAATATTTCCGCCGCCTGGCGAAAGCCGCTTCACGTTTCTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGCTcc                                                                                          >  2:177193/1‑149 (MQ=255)
                                                     cgTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCGTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAGTCTCGATCCGCt                                                                                       <  7:161875/149‑1 (MQ=255)
                                                      gTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGCCCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTg                                                                                      >  1:22362/1‑149 (MQ=255)
                                                      gTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTg                                                                                      >  8:3757/1‑149 (MQ=255)
                                                      gTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTg                                                                                      >  7:57441/1‑149 (MQ=255)
                                                      gTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATATGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTg                                                                                      <  6:181634/149‑1 (MQ=255)
                                                       tGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTgg                                                                                     <  3:326878/149‑1 (MQ=255)
                                                       tGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTgg                                                                                     <  7:3757/149‑1 (MQ=255)
                                                       tGAATATTTCCGCCGCCTGGAGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGAGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTgg                                                                                     <  8:202098/149‑1 (MQ=255)
                                                            atTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACtg                                                                                >  1:358533/1‑149 (MQ=255)
                                                             tttCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACtgt                                                                               <  3:27817/149‑1 (MQ=255)
                                                             tttCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACtg                                                                                <  3:108711/148‑1 (MQ=255)
                                                              ttCCGCCGCCTGGCGAAAGACGCTTCAAGTTACTCTTCAACGATCTCAGACCCTGACACCAATGGGACTCAGTTTAAAGTCCAGCAGCACATTAACGCATACCGCTTCCGTTGTCACCAGCATGCGAATCTCGATCCGCTGGGACtgt                                                                               <  8:77476/148‑1 (MQ=255)
                                                                 cgccgcCTTGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGCCCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGgca                                                                           <  7:186976/149‑1 (MQ=255)
                                                                 cgccgcCTGGCGAATGACGCTTCACGTTACTCTTGAACGATCTGCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCCGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGgca                                                                           <  1:155256/149‑1 (MQ=255)
                                                                 cgccgcCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTgg                                                                             >  3:277240/1‑147 (MQ=255)
                                                                   ccgccTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGgcagc                                                                         <  2:2525/149‑1 (MQ=255)
                                                                     gccTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCaa                                                                       <  2:358533/149‑1 (MQ=255)
                                                                     gccTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGAACCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCaa                                                                       <  3:227888/149‑1 (MQ=255)
                                                                      ccTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAg                                                                      <  3:339818/149‑1 (MQ=255)
                                                                           cGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATaaa                                                                 <  1:10565/149‑1 (MQ=255)
                                                                                gACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGc                                                            >  2:205941/1‑149 (MQ=255)
                                                                                 aCGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGcc                                                           >  7:343070/1‑149 (MQ=255)
                                                                                 aCGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGcc                                                           >  2:387459/1‑149 (MQ=255)
                                                                                       cACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATaaa                                                                 >  6:34920/1‑137 (MQ=255)
                                                                                       cACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATaaa                                                                 <  5:34920/137‑1 (MQ=255)
                                                                                       cACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCt                                                      <  4:356524/148‑1 (MQ=255)
                                                                                              ctctTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGt                                              <  3:85151/149‑1 (MQ=255)
                                                                                                ctTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGGCt                                            >  4:351609/1‑149 (MQ=255)
                                                                                                 ttCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCtt                                           >  5:56158/1‑149 (MQ=255)
                                                                                                 ttCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCtt                                           >  8:154609/1‑149 (MQ=255)
                                                                                                   cAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATACGTCttt                                          >  7:74792/1‑148 (MQ=255)
                                                                                                       gATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCAc                                      >  2:159961/1‑148 (MQ=255)
                                                                                                        aTCTCCGCCCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACg                                     >  4:39968/1‑148 (MQ=255)
                                                                                                         tctcCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGAt                                   >  8:159953/1‑149 (MQ=255)
                                                                                                             cGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCTCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATc                                  <  4:138617/146‑1 (MQ=255)
                                                                                                             cGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGa                               <  2:199433/149‑1 (MQ=255)
                                                                                                                  cTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGaa                          <  4:233332/149‑1 (MQ=255)
                                                                                                                   tGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGAt                                                  >  6:237423/1‑124 (MQ=255)
                                                                                                                   tGACACCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGAt                                                  <  5:237423/124‑1 (MQ=255)
                                                                                                                     acacCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGc                        <  3:351609/148‑1 (MQ=255)
                                                                                                                     acacCAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAg                         >  2:275898/1‑147 (MQ=255)
                                                                                                                     acacCAATGTGAAGCAGGTTAAACTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTGTTTCCACGATCTGACCGAAg                         <  4:123727/147‑1 (MQ=255)
                                                                                                                         cAATGTGAAGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACt                   >  8:395455/1‑149 (MQ=255)
                                                                                                                            tgtgAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCTCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTcc                <  2:115257/149‑1 (MQ=255)
                                                                                                                              tgAAGCAGGATAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATATAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCa               <  1:346761/148‑1 (MQ=255)
                                                                                                                                aaGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAg              >  8:151691/1‑147 (MQ=255)
                                                                                                                                 aGCAGGTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGATTCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGgag           <  4:277240/149‑1 (MQ=255)
                                                                                                                                      gTTAAAGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCtt      >  1:210915/1‑149 (MQ=255)
                                                                                                                                          aaGTCCAGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGAAGCAAGATAAAGTGGCCTATCTGGATCCGTCTTTCCCCGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAAc  >  8:17240/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TCCAGCAGTTACCTGGTACGGGAGTCAAACCGGATCAATTCCACTCTCAAACGCGTGAATATTTCCGCCGCCTGGCGAAAGACGCTTCACGTTACTCTTCAACGATCTCCGACCCTGACACCAATGTGAAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAAC  >  NC_000913/758839‑759125

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