Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,354,537 C→G 6.3% intergenic (‑487/‑188) yhcC ← / → gltB putative Fe‑S oxidoreductase, Radical SAM superfamily protein/glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,354,5370CG 6.3% 224.6 / 4.6 63intergenic (‑487/‑188)yhcC/gltBputative Fe‑S oxidoreductase, Radical SAM superfamily protein/glutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base G (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGCCACGCTCTTTCTGTGCCCGTGCGCAGCGGTTCG  >  NC_000913/3354395‑3354685
                                                                                                                                              |                                                                                                                                                    
aaGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCtc                                                                                                                                                <  2:156489/149‑1 (MQ=255)
aaGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCtc                                                                                                                                                >  2:105819/1‑149 (MQ=255)
 aGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCtca                                                                                                                                               <  2:101097/149‑1 (MQ=255)
    tAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGGCTGCtcatcc                                                                                                                                            <  4:281887/149‑2 (MQ=255)
         tCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATAttt                                                                                                                                       >  3:522147/1‑149 (MQ=255)
          ccATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATATCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTa                                                                                                                                      <  3:512408/149‑1 (MQ=255)
           cATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTAt                                                                                                                                     <  3:295064/149‑1 (MQ=255)
             ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCata                                                                                                                                          <  4:333443/142‑1 (MQ=255)
             ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCata                                                                                                                                          >  3:333443/1‑142 (MQ=255)
             ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGc                                                                                                                                   >  4:131119/1‑149 (MQ=255)
                                        ttttGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGGCTGCtcatcc                                                                                                                                            <  3:176315/113‑2 (MQ=255)
                                        ttttGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGGCTGCtcatcc                                                                                                                                            >  4:176315/1‑112 (MQ=255)
                                              gcTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAAcc                                                                                                  >  3:70415/1‑149 (MQ=255)
                                                tAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTAt                                                                                                       >  2:159846/1‑142 (MQ=255)
                                                tAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTAt                                                                                                       <  1:159846/142‑1 (MQ=255)
                                                  aaGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGc                                                                                                                                                  >  3:555987/1‑97 (MQ=255)
                                                  aaGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGc                                                                                                                                                  <  4:555987/97‑1 (MQ=255)
                                                    gCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCg                                                                                            >  2:529396/1‑149 (MQ=255)
                                                        aTTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGaaaa                                                                                        >  4:270817/1‑149 (MQ=255)
                                                        aTTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGaaaa                                                                                        <  2:600404/149‑1 (MQ=255)
                                                        aTTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGaaaa                                                                                        >  2:9790/1‑149 (MQ=255)
                                                           tCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGa                                                                                     <  1:302407/149‑1 (MQ=255)
                                                           tCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGa                                                                                     <  4:403637/149‑1 (MQ=255)
                                                                aCCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGacaaca                                                                                >  2:618572/1‑149 (MQ=255)
                                                                 ccAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGAcaacaa                                                                               >  4:618350/1‑149 (MQ=255)
                                                                 ccAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGAcaacaa                                                                               <  1:268087/149‑1 (MQ=255)
                                                                  cAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAg                                                                              >  2:496136/1‑149 (MQ=255)
                                                                   aaTAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAgg                                                                             >  1:111856/1‑149 (MQ=255)
                                                                     tAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAgggg                                                                           <  4:222683/149‑1 (MQ=255)
                                                                     tAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAgggg                                                                           <  3:618350/149‑1 (MQ=255)
                                                                       aGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGc                                                                         <  4:288442/149‑1 (MQ=255)
                                                                       aGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTACTGGATGAGCAGCCTTCTCATCCTATTTATGCAGTAATTGAGCTTCACTCTTCACCGTATTAACCGATGCGAAAAGGACAACAATGGGGc                                                                         <  4:200287/149‑1 (MQ=255)
                                                                           tGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGc                                                                                             >  2:619296/1‑125 (MQ=255)
                                                                           tGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGc                                                                                             <  1:619296/125‑1 (MQ=255)
                                                                           tGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAAt                                                                     <  3:270817/149‑1 (MQ=255)
                                                                             ccATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAggg                                                                            >  4:592532/1‑140 (MQ=255)
                                                                             ccATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAggg                                                                            <  3:592532/140‑1 (MQ=255)
                                                                               aTTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGa                                                                 <  4:96487/149‑1 (MQ=255)
                                                                               aTTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGa                                                                 <  4:522147/149‑1 (MQ=255)
                                                                                   gACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGgcg                                                             >  1:383305/1‑149 (MQ=255)
                                                                                    aCCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGgcgc                                                            >  2:389359/1‑149 (MQ=255)
                                                                                      cTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGgcgcgc                                                          <  1:529396/149‑1 (MQ=255)
                                                                                      cTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGgcgcgc                                                          >  3:397785/1‑149 (MQ=255)
                                                                                       tGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTa                                                                                                        >  2:61218/1‑102 (MQ=255)
                                                                                       tGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTa                                                                                                        <  1:61218/102‑1 (MQ=255)
                                                                                           tCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATg                                                     >  1:317819/1‑149 (MQ=255)
                                                                                               cTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGacac                                                 <  1:447993/149‑1 (MQ=255)
                                                                                                tttCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAggggg                                                                          >  1:229258/1‑123 (MQ=255)
                                                                                                tttCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAggggg                                                                          <  2:229258/123‑1 (MQ=255)
                                                                                                 ttCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGc                                               <  1:105819/149‑1 (MQ=255)
                                                                                                    ccGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCaaa                                            >  1:162853/1‑149 (MQ=255)
                                                                                                    ccGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCaaa                                            >  4:506356/1‑149 (MQ=255)
                                                                                                            ttGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGATATTCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATTCGAGGCGCGCGTATGACACGCAAACCCCGTCg                                    >  4:35497/1‑149 (MQ=255)
                                                                                                            ttGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGc                                                                   >  1:600428/1‑118 (MQ=255)
                                                                                                            ttGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGc                                                                   <  2:600428/118‑1 (MQ=255)
                                                                                                            ttGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCg                                    >  1:422858/1‑149 (MQ=255)
                                                                                                              ggAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGcc                                  >  1:170108/1‑149 (MQ=255)
                                                                                                                   tCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGCCACGct                             >  2:103940/1‑149 (MQ=255)
                                                                                                                       cTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGCCACGCTCttt                         <  1:618572/149‑1 (MQ=255)
                                                                                                                           aaGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGCCACGCTCTTTCtgt                     <  1:496136/149‑1 (MQ=255)
                                                                                                                               tttCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGCCACGCTCTTTCTGTGccc                 >  1:51330/1‑149 (MQ=255)
                                                                                                                                 tCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGCCACGCTCTTTCTGTGCCCGt               <  2:422858/149‑1 (MQ=255)
                                                                                                                                        gAGCAGGCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGCCACGCTCTTTCTGTGCCCGTGCGCAGc        >  4:533885/1‑149 (MQ=255)
                                                                                                                                          gcagcCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGTGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGCCACGCTCTTTCTGTGCCCGTGCGCAGCgg      >  4:191868/1‑149 (MQ=255)
                                                                                                                                             gcCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGCCACGCTCTTTCTGTGCCCGTGCGCAGCGGTTc   >  1:336168/1‑149 (MQ=255)
                                                                                                                                              ccTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGCCACGCTCTTTCTGTGCCCGTGCGCAGCGGTTCg  >  3:58054/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
AAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGGCGAATGCGAGGCGCGCGTATGACACGCAAACCCCGTCGCCACGCTCTTTCTGTGCCCGTGCGCAGCGGTTCG  >  NC_000913/3354395‑3354685

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