Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 53,379 T→C 10.5% E13G (GAG→GGG)  pdxA ← 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091353,3790TC10.5% 232.1 / 15.2 76E13G (GAG→GGG) pdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (34/34);  new base C (4/4);  total (38/38)
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.94e-01

GGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTC  >  NC_000913/53235‑53524
                                                                                                                                                |                                                                                                                                                 
ggggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGc                                                                                                                                                >  3:521784/1‑148 (MQ=255)
   gAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCg                                                                                                                                              <  2:589112/147‑1 (MQ=255)
   gAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCg                                                                                                                                              >  1:589112/1‑147 (MQ=255)
    aaTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCggg                                                                                                                                            >  2:705782/1‑148 (MQ=255)
     aTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGa                                                                                                                                           <  4:1044339/148‑1 (MQ=255)
            cGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                                                                        >  2:504521/1‑144 (MQ=255)
            cGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                                                                        <  1:504521/144‑1 (MQ=255)
                    gAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGt                                                                                                                            >  1:280837/1‑148 (MQ=255)
                             aCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTa                                                                                                                   >  1:991695/1‑148 (MQ=255)
                                 aGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCa                                                                                                               >  3:28191/1‑148 (MQ=255)
                                     tCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCa                                                                                                               <  2:376648/144‑1 (MQ=255)
                                     tCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCa                                                                                                               >  1:376648/1‑144 (MQ=255)
                                       gCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGtt                                                                                                         >  1:1010231/1‑148 (MQ=255)
                                         tGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTg                                                                                                        <  3:622158/147‑1 (MQ=255)
                                         tGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTg                                                                                                        >  4:622158/1‑147 (MQ=255)
                                          gCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCt                                                                                                      <  2:1010231/148‑1 (MQ=255)
                                           cccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATc                                                                                                                                     <  4:269118/116‑1 (MQ=255)
                                           cccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATc                                                                                                                                     >  3:269118/1‑116 (MQ=255)
                                              ggTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAgg                                                                                                  <  4:521784/148‑1 (MQ=255)
                                                  ggTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAttt                                                                                              <  4:28191/148‑1 (MQ=255)
                                                     aaGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCACTCCCGGCGAGCCCGCCGGGatt                                                                                                                                         <  3:571379/102‑3 (MQ=255)
                                                     aaGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCACTCCCGGCGAGCCCGCCGGGatt                                                                                                                                         >  4:571379/1‑100 (MQ=255)
                                                     aaGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGa                                                                                                 <  2:236154/142‑1 (MQ=255)
                                                     aaGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGa                                                                                                 >  1:236154/1‑142 (MQ=255)
                                                     aaGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTaa                                                                                           <  1:705782/148‑1 (MQ=255)
                                                         agagTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTa                                                                                       <  2:280837/148‑1 (MQ=255)
                                                         agagTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTa                                                                                       >  1:18470/1‑148 (MQ=255)
                                                              ggCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATca                                                                                                                                    <  3:210770/98‑1 (MQ=255)
                                                              ggCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATca                                                                                                                                    >  4:210770/1‑98 (MQ=255)
                                                                cATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCacaa                                                                                                                                 >  3:1069269/1‑99 (MQ=255)
                                                                cATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCacaa                                                                                                                                 <  4:1069269/99‑1 (MQ=255)
                                                                   cGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTggg                                                                                                                        >  4:922012/1‑105 (MQ=255)
                                                                   cGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTggg                                                                                                                        <  3:922012/105‑1 (MQ=255)
                                                                     gCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTa                                                                                            <  2:456810/131‑1 (MQ=255)
                                                                     gCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTa                                                                                            >  1:456810/1‑131 (MQ=255)
                                                                      cacaAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTa                                                                                                                   >  2:940347/1‑107 (MQ=255)
                                                                      cacaAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTa                                                                                                                   <  1:940347/107‑1 (MQ=255)
                                                                      cacaAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCggg                                                                                                                                            <  4:899886/82‑1 (MQ=255)
                                                                      cacaAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCggg                                                                                                                                            >  3:899886/1‑82 (MQ=255)
                                                                          aacaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                                                                        >  1:736204/1‑82 (MQ=255)
                                                                          aacaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                                                                        <  2:736204/82‑1 (MQ=255)
                                                                           acaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGcc                                                                                                                                               >  2:410660/1‑74 (MQ=255)
                                                                           acaacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGcc                                                                                                                                               <  1:410660/74‑1 (MQ=255)
                                                                            caacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                                                                    <  2:18470/148‑1 (MQ=255)
                                                                            caacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                                                                    <  2:991695/148‑1 (MQ=255)
                                                                             aacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTaa                                                                                                                  >  1:770254/1‑101 (MQ=255)
                                                                             aacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTaa                                                                                                                  <  2:770254/101‑1 (MQ=255)
                                                                                cAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAtt                                                                                                             <  2:104204/103‑1 (MQ=255)
                                                                                cAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAtt                                                                                                             >  1:104204/1‑103 (MQ=255)
                                                                                 aGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTgg                                                                               >  3:160927/1‑132 (MQ=255)
                                                                                 aGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTgg                                                                               <  4:160927/132‑1 (MQ=255)
                                                                                               cTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGctt                                                 >  1:1071788/1‑148 (MQ=255)
                                                                                                       gtgCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACgtt                                                                         >  1:192116/1‑116 (MQ=255)
                                                                                                       gtgCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACgtt                                                                         <  2:192116/116‑1 (MQ=255)
                                                                                                                 gACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGaa                                       >  4:65555/1‑140 (MQ=255)
                                                                                                                 gACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGaa                                       <  3:65555/140‑1 (MQ=255)
                                                                                                                  aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAtt                                                                                               <  1:920219/83‑1 (MQ=255)
                                                                                                                  aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAtt                                                                                               >  2:920219/1‑83 (MQ=255)
                                                                                                                  aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTaa                                                                                           >  2:153421/1‑87 (MQ=255)
                                                                                                                  aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTaa                                                                                           <  1:153421/87‑1 (MQ=255)
                                                                                                                     aCTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcat                           >  3:116708/1‑148 (MQ=255)
                                                                                                                      cTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcatc                          <  2:1071788/148‑1 (MQ=255)
                                                                                                                      cTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTtcatc                          <  4:116708/148‑1 (MQ=255)
                                                                                                                        aaGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTg                                                                                                        >  3:366695/1‑68 (MQ=255)
                                                                                                                        aaGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTg                                                                                                        <  4:366695/68‑1 (MQ=255)
                                                                                                                                   aaTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCg                                           >  1:798148/1‑118 (MQ=255)
                                                                                                                                   aaTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCg                                           <  2:798148/118‑1 (MQ=255)
                                                                                                                                     tCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAtt                                                                                                             >  2:389783/1‑50 (MQ=37)
                                                                                                                                     tCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAtt                                                                                                             <  1:389783/50‑1 (MQ=37)
                                                                                                                                      cccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                                                                    >  2:842748/1‑90 (MQ=255)
                                                                                                                                      cccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                                                                    <  1:842748/90‑1 (MQ=255)
                                                                                                                                      cccGGCGAGCCCGCCGGGATTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCAc          >  1:146580/1‑148 (MQ=255)
                                                                                                                                       ccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGgcgc                                                                                  <  2:937099/75‑1 (MQ=255)
                                                                                                                                       ccGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGgcgc                                                                                  >  1:937099/1‑75 (MQ=255)
                                                                                                                                          gcggGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATc      >  2:774906/1‑148 (MQ=255)
                                                                                                                                              gCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTc  <  1:774906/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGCTTCTTCCGAGAACTTACGGTTCATCAGCATGCGGTATGCACGATCTTTC  >  NC_000913/53235‑53524

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